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THE VARIATION OF MESOTRON INTENSITY ! WITH 
ALTITUDE AND LATITUDE, TOGETHER WITH 
ALLIED PHENOMENA, AND THE BEAR- 

ING OF THESE MATTERS ON THE 
NATURE OF THE PRIMARY 
PARTICLES. 


BY 
W. F. G. SWANN, 


Director, Bartol Research Foundation of The Franklin Institute, Swarthmore, Pa. 
PART I. 


For convenience of publication, and limitation of space, 
the present paper is divided into two parts, the second of 
which will be published in the following issue. Part I will 
contain a statement of the problem, an outline of the pro- 
cedure, and a summary of the conclusions. For convenience 
in correlating the two parts, an outline of the individual sec- 
tions of Part II will be given in the present Part I. 

In several former communications? the writer has had 
occasion to invoke the assumption of an incoming heavy 


1 Several a the omnis established in this paper Minted ae given by the 
writer in condensed form in former publications. Part II of the present paper is 
devoted in part to the detailed proof of some of the formule formerly cited. 

2W. F. G. Swann, Phys. Rev., 56, 209 (1939); Rev. Mod. Phys., 11, 242 
(1939) (see in particular pp. 251-254). Phys. Rev., 58, 200 (1940). JOUR. 
FRANK. INST., 226, 757 (1938). (See, in particular, pp. 780 and 792.) Phys. 
Rev., 63, 210 (1943). 

(Note—The Franklin Institute is not responsible for the statements and opinions advanced 


by contributors in the JouRNAL.) 
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particle which, on collision with the atoms of the air, breaks 
up into smaller particles, now regarded as mesotrons, which 
mesotrons travel through the atmosphere and finally dis- 
appear, either by completing their range through loss by 
ionization or by more sudden disappearance resulting from 
the finite extent of the mesotron life. 

The mesotrons give birth to electrons either by knock-on 
processes closely related to ionization or during their dis- 
integration, in which each mesotron that dies gives birth to 
an electron and a neutrino. The writer has, moreover, given 
reasons for supposing that the aforesaid processes are the only 
processes which are responsible for the appearance of electrons 
in the atmosphere, and that the direct contribution of outer 
space to the electron intensity is negligible. 

2. If we consider a spectral distribution of energies for the 
mesotrons produced at any altitude, the mesotrons of different 
energies will die at different rates on account of the fact that 
the mean life 7 is given by 

7 = 7,/(1 — B?)"?, (1) 
where 7, is the mean life for a stationary mesotron, and 
8 = v/c. Here vis the mesotron velocity and c is the velocity 
of light. Thus the spectral distribution will change with 
depth of penetration into the atmosphere; and the matter is 
further complicated by the fact that on account of energy 
loss by ionization, even the mean life of a mesotron will 
change as it progresses through the atmosphere. 

3. We shall take all of the foregoing considerations into 
account ‘ and in Part II, Section I, we shall consider the 


8 The term “mean life,’’ as applied to a mean life which is changing, presents 
an apparent confusion of terms, which confusion is only temporary, however. 
It is always possible to write for the number of mesotrons 7, in a beam of meso- 
trons of energy e, 


dn, 1 


dt Te 


Be, 


and this equation may be taken to define the “instantaneous mean life” for the 


energy e. 

4 However, we shall consider the ionization loss per unit of path at constant 
pressure to be independent of the energy of the mesotrons. Following the 
development of the formula deduced in Part II, Section J, for the variation of 
intensity with altitude, Dr. M. E. Rose carried out the corresponding calcula- 
tions using the Bethe-Bloch expression for the energy loss. Dr. Rose’s paper, 
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problem of the vertical mesotron intensity at an altitude h 
above sea-level § resulting from production of mesotrons with 
energy distribution defined by the function F(E)dE for the 
range dE, the mesotrons being produced at an altitude where 
the density is p.D., where D, is the sea-level density.* It will 
be shown that the aforesaid vertical intensity is obtained by 
integrating with respect to E the expression 


E + G. — Gi a | 


4 {e-toleroly +1+Go) eh/cro(y+1+4 Go) F(E)dE} (2) 


(An), = 


where m, is the rest mass of the mesotron, x, is the altitude of 
production, A is the constant defined by the atmospheric 
density formula D, = D,exp[— 4], @ is the ionization 
energy loss per centimeter of path at sea-level, and 


G, = ap./d\m,c?; Gy, = ap,/d\m.c? and y = E/m.c?. 


In integrating this equation over the range of E involved, 
we must observe that for any depth x below the altitude of 
production we are concerned only with values of E (at x = 0) 
greater than E,, the loss of energy over the range 0 tox. In 
other words, our lower limit of EF is 


— 
E, = E, = cs (e* — 1) = — (p. — po) 


unless the earth’s magnetic field, through its action on the 
primaries, has imposed a greater lower limit. Thus, if Esp is 
the limit for entry of the primaries, and if the primaries split 
into N mesotrons under the condition that the mesotrons are 
born at rest in the system of axes in which the primary moved, 
then the lower limit of the kinetic energy for the mesotrons will 
be Emup/N. Strictly speaking, it will be even lower than this 
by (1/N) times the energy lost by the primary by ionization in 


which is published in the present issue of the JOURNAL OF THE FRANKLIN INSTI- 
TUTE, covers in general a field of application different from that involved in the 
present paper. 

5 This altitude is supposed expressed in actual distance and not in water 
equivalent distance. 

6 The analytical equivalent of this formula was quoted by the writer, without 


proof, in Phys. Rev., 60, 470 (1941). 
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descending to the level of production. This correction would 
be negligible for altitudes of production sufficiently high to 
have significance in the present theory, particularly on account 
of the factor (1/N) in the correction. 

In Part II, Section I, we shall also deduce the formula for 
the spectral distribution at sea-level in terms of the spectral 
distribution at the level of production. 

4. In Part II, Section II, we shall extend the formula of 
Section I to the case where there is an assigned rate of pro- 
duction of mesotrons at all altitudes. 

5. In Part II, Section III, we shall revert to the formula 
deduced in Part II, Section I, with a view to seeing how far 
the facts may be invoked as a test of the hypothesis of 
mesotron production by the disintegration of heavy primaries. 
Our calculations will involve an assumption of the spectral 
distribution formula for the mesotrons, and we shall obtain 
this distribution by supposing that the mesotrons have a 
spectral distribution defined by an F(£) proportional to 7/E*. 
Such a conclusion would result from the assumption that the 
true primaries are particles heavier than mesotrons,’ with an 
inverse cube spectral distribution, and that on entering the 
atmosphere they split up into mesotrons.in such fashion that 
the mesotrons born from any primary are born at rest in the 
system of axes in which the primary was moving. Indeed, this 
last assumption results, firstly, in the mesotrons perpetuating 
everywhere the direction of motion of the primary particles, 
and secondly, in the expression K’dE/E* applicable to the 
primaries being replaced by the expression K’'dE/N?F? for the 
mesotrons, where NV is the number of mesotrons produced 
from a primary, and where the lower limit of the spectral dis- 
tribution for the mesotrons is 7/N times the lower limit for 
the primaries. The value of 8(= v/c) for a mesotron pro- 
duced as above is obviously the same as that for the parent 
primary. 


7 In addition to the citations from the present writer’s former publications, 
the postulation of protons as primary particles has also been made by A. H. 
Compton and H. Bethe, Nature, 134, 734 (1934), by A. H. Compton, Rev. Sci. 
Insir., 7, 71 (1936), by T. H. Johnson, Phys. Rev., 45, 569 (1934), and by M. 
Schein, W. P. Jesse and E. O. Wollan, Phys. Rev., 59, 615 (1941). In most of 
these cases the protons were regarded as responsible for only a part of the incoming 
radiation, electrons being invoked in addition. 
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6. As regards experimental data for the purposes of test, 
we have the observations of Schein, Jesse, and Wollan * on 
the variation of vertical mesotron intensity with altitude. 
Then we have a certain amount of data on the sea-level lati- 
tude effect. We shall endeavor to utilize this data to fix the 
unknown quantities, altitude of mesotron production and 
lower limit of energy and so of 8 in the mesotron spectral 
energy distribution. On the basis of the hypothesis of meso- 
tron formation cited above, this value of 8 is also the lower 
limit of 8 for the primaries. Now the lower limit of the 
momentum of the primaries for vertical entry is known from 
the magnetic field data for any latitude, and so the mass of 
the primary becomes determined. It turns out that, within 
the limits of accuracy of the experimental data, the primary 
may be taken as a proton, or at least that the primaries re- 
sponsible for the mesotrons observed from sea-level up to 
semi-high altitudes (Pike’s Peak, for example) are protons. 

7. In Part II, Section III, we shall also make a comparison 
of the experimental data for the spectral energy distribution 
at sea-level and at magnetic latitude 50°, with the results pre- 
dicted by the theoretical formula deduced in Section I, using 
for that purpose the altitude of production assigned by the 
considerations cited above. It may be remarked that the 
lower limit for 8 does not become involved in these calculations 
since mesotrons with that lower limit could not reach down to 
sea-level. 

8. In Part IT, Section [V, we consider a matter in a category 
to some extent different from those discussed in Sections [-II]. 
It concerns the angular distribution of the electrons emitted 
by the mesotrons at death. In the system of axes which 
travels with a moving mesotron, the electrons are emitted, 
on the average, symmetrically in all directions. When viewed 
from stationary axes, however, the emission favors the direc- 
tion of propagation of the mesotron just prior to death. In 
Part II, Section IV, the general problem of electron emission 
at different angles is discussed and, in particular, expressions 
are derived for the horizontal and vertical electron intensities 
resulting from the death of vertically directed mesotrons. 


8M. Schein, W. P. Jesse, and E. O. Wollan, Phys. Rev., 59, 615 (1941). 
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Now at high altitudes the horizontal intensity (presumably 
an electron intensity) is so great * as to imply a considerable 
symmetry in electron emission at mesotron death when 
measurements are made from fixed axes. Indeed, it was the 
large horizontal intensity which caused the writer 2 to propose 
the view that the electrons in the upper atmosphere come 
mainly from the death of mesotrons, and that only a negligible 
fraction can be attributed to incoming primary electrons. 
However, it turns out that the lower limit of the mesotron 
velocities necessary to account for the sea-level latitude ratio 
and the intensity-altitude data is too high to permit sufficient 
symmetry in electron emission at mesotron death to provide 
for the observed horizontal intensity at high altitudes in the 
case of proton born mesotrons. In order to obtain a suffi- 
ciently low mesotron velocity to provide for the observed sym- 
metry, it is necessary to provide for a correspondingly low 
velocity in the incoming primaries which give rise to the 
mesotrons, and such a low velocity demands a mass larger 
than the proton mass for these primaries in order that they 
may enter through the earth’s magnetic field. It is shown 
that the conditions are satisfactorily provided for if singly 
charged helium atoms are invoked for the primaries. It is 
inferred, therefore, that in addition to the protons already 
invoked to explain, through their progeny, phenomena at low 
altitudes and at semi-high altitudes (Pike’s Peak, for example) 


* It has come to the writer’s attention that some doubt has been expressed 
(not in writing, but in conversations) as to the validity of these results, the criti- 
cism being to the effect that the observations may have resulted from shower 
rays passing simultaneously through the coincidence counters. It is, therefore, 
worth while to point out that, in addition to using triple coincidence counter 
telescopes, the individual units of the triplet sets were protected from showers 
such as the above by a wall of anti-coincidence counters. This was done in the 
case of half of the ten telescope assemblies, the other half being unprotected 
against showers for comparison. A full description of the details is given in 
published papers on this matter; and in the case of the flight of Explorer II, 
the intensity-altitude curves are given for the protected and the unprotected sets 
for various zenith angles ranging from zero to 90°. See W. F. G. Swann and 
G. L. Locher, Nat. Geog. Soc. Contr. Tech. Pap., Stratosphere Series, No. 1, 
pp. 7-14 (1935); also, W. F. G. Swann, G. L. Locher, and W. E. Danforth, 
Nat. Geog. Soc. Contr. Tech. Pap., Stratosphere Series, No. 2, pp. 16-25 (1936); 
also, W. F. G. Swann, Jour. FRANK. INsT., 222, 1 (1936) (see, in particular, 
PP- 45-55). 
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we may have to invoke singly charged helium primaries in 
order to account for phenomena observable at high altitudes 
where the horizontal intensity becomes important. The 
slowly moving mesotrons born of the helium atoms die, of 
course, on account of their short mean lives, before they have 
penetrated far into the atmosphere. 

Of course, it must be realized that all cf the considerations 
involved in all but Section II, applicable, as they are, to 
only a single altitude of mesotron production, are limited in 
approximation by that circumstance. They tend to validity 
only in places reasonably far below the region where produc- 
tion is appreciable, and there only to the extent that we can 
invoke a kind of “‘center of gravity’’ of production as an 
adequate representative of the facts. Naturally they do not 
contain the story of the maximum of mesotron intensity to 
be expected at some altitude in the region of production, since 
that phenomenon involves in its representation the law of 
disappearance of primaries with penetration into the atmos- 
phere. 


(To be continued in August issue.) 
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Crystals Used to Measure High Accelerations. (Westinghouse 
Engineer, Vol. 3, No. 1.) The use of quartz spans the entire history 
of mankind. Among the few remaining artifacts of earliest man 
are arrowheads of flint, a hard form of quartz. The frequencies of 
radio broadcast stations are controlled by quartz crystals. And 
now an ultra-modern use—the measurement of acceleration—has 
been made of their fascinating physical and electrical properties. 
Quartz has a mysterious ability. When an alternating voltage is 
applied between the faces of a quartz crystal, properly cut, the crystal 
vibrates mechanically at the frequency of the applied voltage. 
This is the property utilized in radio frequency control. But this 
action can be worked in reverse. Vibrate the crystal mechanically, 
and an alternating voltage of the same frequency appears across the 
faces. This fact has been put to work at the Westinghouse Research 
Laboratories as a new kind of accelerometer. This device is in- 
tended to measure high rates of accelerations. It consists of a 
steel case bolted solidly to the object whose acceleration is being 
measured. On the inside, stacked in line.in the direction of motion, 
are a quartz crystal, a thin layer of copper, another quartz crystal 
and a small mass all held against the case by a soft spring. When 
the acceleration occurs, this mass, because of its inertia, exerts a 
force on the quartz crystals. Because of their peculiar ability the 
crystals translate this force into a charge proportional to the 
acceleration. This charge is taken off at the copper leaf, led to an 
amplifying system, and thenee to a cathode ray oscillograph, where 
it is photographed on a moving film. The film speed is such that 
the events of one thousandth of an inch are spread over an inch of 
film. The quartz accelerometer faithfully records accelerations of 
high order occurring at frequencies up to 10,000 per second, which 
is well below its own natural frequency. The accelerometer itself 
is small enough to be carried in a coat pocket, being about 1{ inches 
square by 3 inches long. Its special merit is ruggedness. While 
the output of quartz is not as high as some other crystals, it is able 
to withstand without injury the high forces associated with high 
accelerations. 
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ON THE ABSORPTION OF THE HARD COMPONENT OF 
THE COSMIC RADIATION. 


BY 
M. E. ROSE, 


Bartol Research Foundation of The Franklin Institute, Swarthmore, Pa.* 


ABSTRACT. 


A phenomenological theory of absorption of mesotrons in 
the atmosphere and in dense materials at and below sea-level 
is given. The calculations are based on the following as- 
sumptions: (1) absorption due entirely to energy loss caused 
by electric interactions of charged particles and to spon- 
taneous decay, (2) the localization of mesotron production in 
an equivalent layer 10 per cent. below the top of the homo- 
geneous atmosphere, and (3) neglect of scattering. From 
the absorption of vertically incident mesotrons in thick layers 
of rock the energy distribution in the formation layer is 
deduced. The energy distribution at sea-level is then de- 
termined, without approximations in the analysis, and com- 
pared with observations. It is shown that fair agreement is 
obtained for any reasonable value of the proper lifetime but 
that a value of about 2 X 107° sec. allows the most favorable 
comparison of the calculated results and the data. The 
integral energy distribution and the angular distribution at 
sea-level are obtained. For the latter it is assumed, as is 
plausible, that the initial angular distribution is isotropic and 
that the primary energy distribution is independent of the 
original direction of motion. The agreement between calcu- 
lated and observed angular distributions is very satisfactory. 
The absorption in lead and in light materials for vertically 
incident mesotrons at sea-level is considered. The need for 
a more precise definition of equivalent thicknesses is pointed 
out and the stopping powers of lead and light substances are 
compared. The absorption in lead is calculated and com- 
parison with the data indicates a slightly larger absorption 


* Now at Princeton University. 


10 M. E. Rose. i. Pe f. 


than is found theoretically. However, the absorption at 
large lead thicknesses (over 100 cm.) and the qualitative 
behavior of the absorption coefficient with increasing depth 
below sea-level are in accord with experiment. The considera- 
tion of absorption above sea-level indicates that there should 
be appreciable mesotron production at elevations comparable 
to that of Mt. Evans (4,300 meters). This is borne out by 
recent observations. 


I. INTRODUCTION. 


As a result of a considerable number of investigations 
carried out in recent years the nature of the hard, or pene- 
trating, component of the cosmic radiation is now fairly well 
established. These radiations are known to consist of charged 
particles,! mesotrons, with a mass of roughly 200 times the 
mass of the electron. Moreover, it is now quite certain that 
the mesotrons are unstable and disintegrate spontaneously in 
their passage through the atmosphere. From this fact alone 
it follows that these particles are of secondary origin and are 
produced in the atmosphere—mainly in the stratosphere. 
Concerning the nature of the primary particle and the 
mechanism of the production process very little can be stated 
with certainty, except that the process is probably a nuclear 
one.’ 

The process of production constitutes one of the principle 
problems in the analysis of the hard component. A second 
problem is concerned with the subsequent behavior of the 
mesotrons once they have been produced. It is clear that 
these two problems are independent except insofar as the 
production process determines the initial energy and angular 
distributions of the mesotrons, and even this information is 
deducible from experiments on the absorption of the hard 


component. 


1 For the considerations presented here the sign of the charge is immaterial. 
Actually there is a slight excess of positively charged particles. In addition the 
existence of neutral mesotrons has been postulated but not as yet confirmed. 
Our considerations are also independent of the question of the existence of such 


uncharged particles. 
2W. F. G. Swann, Phys. Rev., 59, 770, 836 (1941); T. H. Johnson, Rev. 


Mod. Phys., 11, 208 (1939). 
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Since the mesotrons are charged particles it is to be 
expected that their behavior in passing through matter is to 
be described by the well-known collision processes to which 
such particles are subject; viz., ionization and _ radiative 
collisions. These processes, in addition to the aforementioned 
spontaneous decay should determine the manner in which the 
particles are absorbed. While it is true that a mechanism of 
absorption involving nuclear collisions occurs, as is evidenced 
by the appearance of ‘‘stars’’ on photographic plates,* these 
effects are rare and do not play an essential role in the cosmic- 
ray phenomena to be discussed here. 

Once the mechanisms of absorption and, in particular, 
that of energy loss have been specified (cf. Section ITI, 
below), it becomes possible to determine the energy distribu- 
tion (and angular distribution) as well as the total number of 
particles arriving at any given depth below the point of 
formation of the particles, provided the energy (and angular) 
distribution as well as the spatial distribution at the moment 
of formation are known. In the absence of information 
concerning the mechanism of production, the latter informa- 
tion could be obtained only from direct experimental data. 
Unfortunately, direct information as to the distribution of the 
mesotrons at birth is not now available. Therefore one may 
attempt to solve the converse problem: to obtain information 
as to the place of origin and the distribution in energy and 
angle of the produced mesotrons from the existing data on 
absorption in the atmosphere and below sea-level. Such a 
program could be carried through if the experimental data 
were sufficiently complete. While this is not quite the case, 
the information that is available is sufficient to permit what 
may be termed a first-order solution of the problem. From 
the results obtained it will be possible to draw several conclu- 
sions of interest. At the same time it is the purpose of the 
present paper to indicate in what directions the measurements 
should be extended in order that more complete information 
concerning the process of mesotron production may be 
obtained. The extent to which the theory can be carried at 
present is indicated in the next section. 


3W. M. Powell, Phys. Rev., 61, 670 (1942). See also Section III of the 
article by M. M. Shapiro, Rev. Mod. Phys., 13, 58 (1941). 
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II. PRIMARY MESOTRON DISTRIBUTION. 


The experimental data thus far obtained are sufficient to 
determine the energy distribution of the mesotrons at their 
place of origin.t The measurements which provide this 
information are the following: 


(a) The energy distribution of vertically incident meso- 
trons at sea-level. 

(6) The absorption of mesotrons in rock (or water) up to 
depths below sea-level equivalent in mass to as much 
as I,100 meters of water. 

(c) The angular distribution of mesotrons of all energies 
at sea-level.‘ 


The manner in which these measurements are to be used is 
considered in detail in Sections IV and V. An indication of 
the consistency of the results may then be obtained from the 
comparison of theoretical results with observations in the 
case of further experiments. The angular distribution at 
sea-level is considered in Section VI. In Section VII the 
results previously obtained are applied to a consideration of 
the absorption in lead and other materials for the case of 
vertically incident mesotrons at sea-level. In Section VIII 
the absorption in the atmosphere is considered. 

It may be noted that no conclusions regarding the place 
of formation, or the spatial distribution of mesotrons at 
formation, are advanced. In fact, since the intensity of the 
hard component in the atmosphere or underground depends 
not only on the distance traversed by the rays but also on 
the shape of the primary energy distribution, it may appear 
surprising that one can obtain information concerning the 
one and not the other. By way of explanation it may be 
pointed out that the original energy distribution of the 
particles is determined primarily from the absorption in very 
thick layers of rock (see (b) above). Because of the large 
thicknesses involved and the extremely large energies required 
to penetrate such thicknesses it does not matter much where 
the rays are formed, since a completely negligible energy loss 


‘The angular distribution may also be determined. Alternatively, one can 
make a plausible assumption for the initial distribution in angle and compare 
the calculated and observed distributions at sea-level to check the theory. 
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and absorption takes place in the atmosphere. Hence the 
energy distribution can be obtained from the differential 
absorption without assumptions as to the place of formation. 
It is to be emphasized however, that only the distribution for 
primary energies greater than about 10!° e.v. is fixed in this 
way because particles of lower energy are not transmitted 
through the very thick layers of dense material. The energy 
distribution for lower energies is obtained from the absorption 
in thinner layers of material and from the energy distribution 
as observed at sea-level. However, this still leaves the 
primary distribution undetermined for energies below 2 < 10° 
e.v., since slower particles, in general, are absorbed before 
reaching sea-level. 

. While no detailed information is at hand concerning the 
spatial distribution of the mesotrons at the moment of 
formation, it appears that production takes place at all 
elevations throughout the atmospheric layer. However, there 
is evidence for production with appreciable intensity only 
above 4,000 meters or thereabouts.*> Moreover, the intensity 
of production increases so rapidly with elevation that for 
most purposes it is sufficient to replace the entire producing 
region by an equivalent layer. That is, one assumes that all 
mesotrons are produced in a certain layer of fixed depth 
below the top of the homogeneous atmosphere.’ It is 
customary to take the position of this equivalent layer at 10 
per cent. below the top of the homogeneous atmosphere: at 
an elevation of 9.3 meters water equivalent above sea-level. 
As a justification for this procedure it may be stated that the 
energy distribution and/or the total intensity at sea-level 
and below is insensitive to the choice of position of the 
equivalent layer as long as this layer is placed near the top of 
the homogeneous atmosphere. It is obvious that shifting 
the position of the layer from which the particles start will 
have virtually no effect on the absorption due to energy loss 
since this depends on the amount of matter traversed. A 
somewhat greater effect on absorption due to decay would be 


5M. Schein, W. P. Jesse and E. O. Wollan, Phys. Rev., 59, 615 (1941). 

6 That is, the atmospheric layer considered as compressed in such a way 
that the pressure is uniform and equal to the normal sea-level value of 76 cm. 
Hg. The height of the homogeneous atmosphere is thus 8.35 X 10% meters. 
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expected (dependent on distance traversed) but, as will be 
seen in the following, the effect of decay does not appreciably 
change the shape of the energy distribution nor the shape of 
the absorption curve in dense materials. However, when one 
attempts to extend the calculations to absorption in the 
atmosphere and at the same time retaining the assumption of 
an equivalent producing layer, no agreement between theory 
and observed results can be obtained (see Sec. VIII). This 
is to be expected since, for elevations where mesotron pro- 
duction is of some importance, the results are much more 
sensitive to the assumed position of the equivalent producing 
layer. 
Ill. ENERGY LOSS AND DISINTEGRATION. 

For the absorption of mesotrons in the atmosphere or for 
the penetration through dense materials of low atomic number 
the energy loss of the particles is almost entirely determined 
by the ionizing collisions. Energy loss due to radiative 
collisions become of equal importance only for very high 
energies; namely, for energies of the order 2 X Io” e.v. in 
light materials such as rock, and for about 3 X 10" e.v. in 
lead. In addition, as will be shown below, the intensity of 
such mesotrons, whose energy is so great as to experience 
appreciable radiative losses, depends only slightly on the total 
energy loss. Therefore, radiation energy losses will con- 
tribute only a small correction term to the intensity. 

The precise form for the rate of energy loss depends on 
the spin of the mesotron and it is not yet possible to assign 
a value of the spin with any degree of assurance. However, 
in the case of energy loss due to ionization the collision 
probability, and therefore the rate of energy loss, becomes 
independent of the spin when the energy given to the second- 
ary electron is small compared to the mesotron energy. 
This will be the case in the vast majority of events since 
large energy transfers result from very close collisions (scat- 
tering through large angles), and such collisions are exceed- 
ingly rare due to the long range nature of the Coulomb forces 
responsible for the scattering.?. In the case of radiative 


7 Cf., for example, B. Rossi and K. Greisen, Rev. Mod. Phys., 13, 240 (1941). 
It is also necessary that the energy transfer be small compared to 2 X 10" e.v. 
which will be fulfilled in practically all cases of interest. 
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losses the question of spin dependence does not concern us in 
view of the remarks made above. 

In the following it will be convenient to measure distances 
in terms of mass traversed. Thus, if ¢ represents a distance 
in g/cm.? or cm. of water, and x represents the same distance 
in conventional units (cm.) we have 


dt = pdx, (1) 


where p is the density of the material under consideration. 
It is also convenient to introduce the energy of the mesotron 
in units of the rest energy. Thus if & represents the energy 
in c.g.s. units we introduce . 


e = E/uc?, (2) 


where pw is the rest mass of the mesotron. For uw = 200m 
(i.e., w = 1.8 X 107%5 g.) the rest energy of the mesotron is 
just 10° e.v. Thus the energy £ is given by 10%¢ e.v.§ 

The rate of energy loss by ionization is now given by ° 


=) 2mry?>Z m 2mc?B?W , | 
B = a = : ae —-— 2 2 A 
Ce) (3 Amn wel Raa — 6 2) 


Here ro = e?/mc? is the classical electronic radius, e and m 
being the electronic charge and mass; Z and A are the atomic 
number and mass number respectively; mm is the mass of the 
hydrogen atom, 1.66 X 10-*4 g., so that Amy is the atomic 
weight. J(Z) is the average ionization energy of the atom of 
atomic number Z. Following Bloch '° we set 


I(Z) = 13.52 e.v. (3a) 


8 Recent determinations of the mass would indicate that u = 180m is a 
more reliable value. Inasmuch as this work was completed before the better 
mass determinations appeared the value » = 200 m had been adopted. However, 
the error thereby induced is, in general, unimportant and is certainly smaller 
than the experimental errors involved in absorption measurements, etc. See 
also Sec. IV. 

°Cf., for example, J. A. Wheeler and R. Ladenburg, Phys. Rev., 60, 754 
(1941). We neglect the so-called density effect investigated by Fermi, Phys. 
Rev., 56, 1242 (1940); 57, 485 (1940). See also O. Halpern and H. Hall, Phys. 
Rev., 57, 459 (1940), and reference 7, p. 247. 

10F. Bloch, Zeits. f. Phys., 81, 363 (1933). 
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As usual 8 = v/c where v is the mesotron velocity, and we have 
e(e + 2) 
(e+ 3)" 


W,, is the maximum energy that may be transferred to the 
secondary electron 


? = (3b) 


W 2myp*c2(e + 2) (3c) 
5 Eee tarrnenntnciniaaemninnen: @, 3c 
m? + uw? + 2mp(e + 1) 
The rate of energy loss due to radiative collisions is 
given by !! 


Bio = - (4) 


= - rit_(* ) (e + » [xe + 4log4z va] (4) 


137Amn \ wu m 
where 


and x(£) is a function given in Fig. 2 of reference 11. The 
total rate of energy loss is then 


B(e) = Bile) + B,(e). (5) 


It may be noted that for a given velocity the ionization loss 
dE /dt is independent of the mesotron mass for all but the 
very high energies, viz.: for « <u/m or for energies small 
compared to 2 X 10!" e.v., W,, and therefore dE/dt depends 
only on the velocity. £B,(e) is then inversely proportional to 
the mesotron mass for a given velocity. This fact has the 
important consequence that under conditions in which only 
ionization loss is responsible for absorption, the absorption 
coefficient depends on the mass in a very simple manner 
(Sec. IV). 

The range of a mesotron with initial energy €» is given in 
terms of the stopping function B by 
de 


"a ae (Sa) 


R(€o) 


11 Cf., for example, H. J. Bhabha, Proc. Roy. Soc., A164, 257 (1938). Eq. 
(4) applies to particles of spin 3. 
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From the remarks made above it follows that for all but very 
high energies (€ > 200) the range is proportional to the 
mass uw. Curves giving the rate of energy loss due to ion- 
ization for air, Fe and Pb up to e = 100 are given in Fig. 2 
of reference 7. Range-momentum curves are given in Fig. 3 
of the same reference. 

In addition to the attenuation due to energy loss the 
intensity of the mesotron component of the cosmic radiation 
decreases as a result of decay." If 7 is the lifetime of the 
mesotron measured in a frame of reference moving with the 
particle, the fractional decrease of intensity due to decay in 
traversing a distance dx is 


= af (6) 


In terms of the lifetime +> measured in the laboratory frame 
of reference this becomes 


di. ( zy" dx 6 

Tl tt dor cpro’ — 
since tT = T9(1 — v?/c?)~'/*. Here p is the momentum of the 
mesotron in units of wc. For the consideration of absorption 
in the atmosphere it will be desirable to express this attenua- 
tion due to decay in terms of lengths measured in g./cm.? 
(cm. of water). Using the barometric formula 


t = Le z/ Zs 


where ¢, is the atmospheric thickness (10.3 meters water 
equivalent) and x, = ¢t,/p,, where p, is the sea-level density of 
air at 15° C. and 76cm. Hg. Then (6a) becomes 


dI_t, dt _ adt (7) 
I tpcrop tp 

The numerical value of a depends on 7, the value of the 
lifetime of the mesotron at rest. The experimental determi- 


12 A partial summary of the ample evidence for decay is given by B. Rossi, 
Rev. Mod. Phys., 11, 296 (1939). More recently the lifetime has been measured 
by F. Rasetti, Phys. Rev., 60, 198 (1941), and G. Bernardini, B. N. Cacciapuotti, 
E. Pancini, O. Piccioni and G. C. Wick, Phys. Rev., 60, 910 (1941), which may 
be consulted for further references. 


18 M. E. Rose. [J. F: 1. 


nations of ro give values ranging from 1.3 to about 5 X 107° 
sec. The weight of the experimental evidence, however, 
points to a value lying between 2 and 3 X I10-® sec. The 
procedure adopted here is to choose that value of t» which 
gives the best fit between the observed and calculated energy 
distribution for vertically incident mesotrons at sea-level. 


IV. SEA-LEVEL ENERGY DISTRIBUTION FOR VERTICAL INCIDENCE.': 


On the assumption that the ordinary processes of energy 
loss due to electromagnetic interactions and the decay suffice 
to describe the behavior of mesotrons, the energy distribution 
of the particles at any depth below the formation layer may 
be determined with the following auxiliary assumptions: 

(1) The specification of the primary energy distribution. 
Actually this is determined from the absorption of the 
mesotrons in thick layers of rock !* and, properly speaking, 
the calculation of this absorption curve involves a knowledge 
of the sea-level distribution, which in turn depends on the 
primary distribution in energy. However, since those meso- 
trons with sufficient energy to penetrate thick rock layers do 
not undergo appreciable decay (Eqs. 6 and 7), it is permissible 
to replace the atmosphere by a layer of equivalent stopping 
power. The absorption is then entirely due to energy loss 
and it becomes a simple matter to determine the energy 
distribution which gives rise to the observed absorption 
(Sec. V). For that part of the primary spectrum which is of 
importance for the sea-level energy distribution we have for 
the number of particles formed per sec. per cm.” with energy 
between ¢« and e + de with the direction of motion confined 
to the solid angle dw at the zenith 


F(e)dedw = Coe~%deda, (8) 
where C) is a constant to be adjusted empirically and y = 2.90. 


This result does not apply for very low energies, i.e., for 
¢ < 20 which is the minimum energy necessary to traverse 


13 The results presented in this section were also obtained by D. R. Hartree 
(cf. reference 12). Hartree’s calculations are unpublished, and were repeated 
here for the purpose of an independent check. The details are presented for 
the purpose of illustrating the procedure used in this and in other parts of the 


work. 
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the atmosphere. Moreover, the distribution (8) does not 
apply for very large energies: « > 800 (see Eq. 17). For 
such energies the primary intensity and the intensity at 
sea-level are too small to observe, and in fact, beyond the 
range of curvature measurements. Thus (8) applies to such 
energies as are of pertinence for the sea-level distribution. 

(2) The foregoing implies that the mesotrons are observed 
at sea-level only after having traversed a large part of the 
atmosphere. This is true for most of the mesotrons detected 
at low elevations and, in accordance with the discussion of 
Section II, we assume that all the observed mesotrons origi- 
nate in a thin layer near the top of the homogeneous atmos- 
phere. We thereby overestimate the decay due to particles 
produced at lower elevations and underestimate the decay 
effect for particles originating above the assumed position of 
the formation layer. The effect of the assumption of a 
production layer in the upper atmosphere on ionization loss 
is not serious since most of the energy loss must take place at 
very low elevation where the atmospheric density is greatest. 
The height of the production layer will be taken at 10 per cent. 
below the top of the homogeneous atmosphere, that is, the 
layer at which the atmospheric pressure is 76 mm. Hg. 

(3) Finally, we assume that scattering may be neglected 
so that the mesotrons observed in vertical incidence, for 
example, have described a vertical path from the point of 
origin to the point of observation. This assumption should 
be valid (on account of the strong forward scattering) for 
all but the very low energies (e = 1), and these low energy 
particles do not contribute greatly to the total intensity. 
An estimate of the effect of scattering may be made and it is 
concluded that the maximum error caused by its neglect is 
about one percent. 

With these assumptions we may proceed to the calculation 
of the energy distribution at sea-level. The considerations 
of this section will be confined to mesotrons incident vertically. 
The case of oblique incidence is considered in Section VI. 
The number of mesotrons incident per cm.? per sec. in the 
unit solid angle at the zenith with energy between ¢ and 
e + de at a point ¢ cm. water (g./cm.?) below the top of the 
homogeneous atmosphere will be denoted by F(e, t)de. With 


20 M. E. Roser. [j. F. 1. 


the origin of ¢ at the top of the homogeneous atmosphere we 
have ¢t = t, = 1033 g./cm.? at sea-level and ¢ = ty) = 103.3 
g./cm.* at the formation layer. 

Consider first the case in which decay is neglected and the 
absorption is entirely due to stopping by energy loss. By 
virtue of the assumption of straight line paths (see (3) above), 
the energy loss from the point of production is known. 
Denoting the energy loss in units of uc? by Ae we have 


aa ke de (9) 
a ee ‘ 


Thus, for the atmospheric layer between the point of produc- 
tion and sea-level, or for any other layer of absorbing material, 
we may define a pair of corresponding energies: the energy « 
at the bottom of the layer (point of observation) and the 
energy « = e+ Ae at the top of the layer. The loss Ae, 
and therefore e’, is a slowly varying function of e. 

It is now clear that the number of mesotrons whose energy 
lies in any interval (e, € + de) at sea-level is equal to the 
number of mesotrons whose energy lies in the corresponding 
interval (e’, e’ + de’) at the position of the production layer. 
Hence the distribution at sea-level F(e, t,)de is obtained by 
shifting the distribution in the production layer Fo(e’)de’ 
towards smaller energies by the amount Ae, and taking into 
account the change in the energy interval de’ — de: 


F(e, t,)\de = Fy(e’)de’, (10) 
or, from (9), 
de’ aa 
- 9 Hs _ ; ‘) =- . F of ‘ s 
F(e, t,) i (e’) Ble) €’) 
Then from (8) we have 
B A ; 
F(e, t,) = Co aD ) (e + Ae)~7, (11) 


where B is given by (3), (4) and (5). 

The effect of decay is now to introduce an attenuation 
factor which is readily obtained from (7). Taking into 
consideration the energy dependence of the probability for 
decay and the fact that the energy is changing during the 
traversal through the atmosphere, we find for the energy 


July, 1943-]} HaArbD COMPONENT OF CosMIc RADIATION. 21 


distribution at sea-level: 


: _B A ; 
F(e,t.) = Co — é) (e + Ae)~e-=, (12) 
where 
“Ro dR b, 
st =J. tp(R)’ ae psCTo ceca 


In (12a) we have introduced the residual range R as the 
variable of integration. The limits of integration Ry and R, 
are the ranges in air (of normal density) for the corresponding 
energies «’ = €, and e, respectively (cf. 5a). Obviously, 
we have 


It is evident that the result (12), aside from the trivial 
normalizing constant Cy, completely defines the energy distri- 
bution at sea-level for the vertical rays. However, it may be 
readily appreciated that specific results are obtainable only 
by a rather tedious numerical procedure. The decay integral 
(12a), which depends on the momentum-range relationship— 
Eq. 5a, cannot be evaluated analytically without simplifying 
assumptions; moreover, a separate numerical integration 
must be carried out for each value of «. <A great part of the 
labor involved may be eliminated if one adopts the simplifica- 
tion introduced by Euler and Heisenberg." Here one assumes 
a constant rate of energy loss and this allows an exact expres- 
sion for the momentum-range relationship and an analytic 
evaluation of the decay integral. The assumption of a 
constant energy loss, however, while fairly suitable for 
intermediate energies is a poor approximation for small 
energies (e 5 1). As is evident from (3) or from the figure 
referred to above,’ the rate of energy loss increases very 
rapidly for energies just below the rest energy. Hence the 
Euler-Heisenberg approximation cannot be expected to pro- 
vide a valid description for such phenomena as decay, which 
depend primarily on slow particles. While it is true that 
iccurate measurements of the intensity of slow particles, such 
as could be provided by the cloud chamber or proportional 


4H. Euler and W. Heisenberg, Ergeb. d. Exakt. Naturwiss., 17, 41 (1938). 
See especially pp. 39-44. 
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counter, are very difficult and, in fact, are not yet at hand, 
the importance of the decay phenomenon makes it worth 
while to consider a more accurate solution of the problem. 
In addition, it is of interest to compare the observed results 
with a theoretical analysis based on as few assumptions as 
possible in order to obtain a better indication of the im- 
portance of the assumptions concerning the primary distri- 
bution of mesotrons which forms the basis of the present 
calculations. Therefore, the simplification of constant energy 
loss will be avoided and the results are to be obtained by a 
rigorous calculation. 

In carrying out the calculations it is, of course, necessary 
to specify the numerical value of the constants which enter. 
Of these the mass and lifetime, u and ro, are uncertain with 
the experimental uncertainty in the case of u (about 20 per 
cent.) being considerably smaller than that of 7». In the 
procedure adopted here the mass was taken as pn = 200m!‘ 
and the lifetime tr» was adjusted so as to obtain as good 
agreement as possible with the observed energy distribution 
at sea-level. It is therefore pertinent to consider in what 
manner the results may be changed if another value of the 
mass were adopted. For this purpose in the energy distri- 
bution (12) we consider separately the factors directly due to 
energy loss and the exponential factor arising from the decay. 
In the factors arising from the energy loss the mass enters in 
the Ae which is essentially inversely proportional to the mass 
(cf. remarks preceding Eq. 5a). Since the energy loss Ae is a 
very slowly varying function of e, a decrease in mass u would 
mean, cet par. a distribution function varying less rapidly 
with energy. This will be more noticeable at energies for 
which Ae = €; i.e., € 5 20. In addition, the maximum of 
the energy distribution, see Fig. 1 below, arising from the 
rapid decrease of the energy loss factors in (12) and the 
increase in the decay exponential factor with e, would shift 
toward slightly higher energy, provided the decay factor were 
unaltered. Since the effect of a small change in mass (from 
200m to 180m for example) would have a rather unimportant 
effect on the energy loss factors we may conclude that the 
results given below will apply for other masses in the neighbor- 
hood of the chosen value if, as before, the decay factor remains 
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unchanged. The decay factor, as may be seen from (12a), 
is essentially proportional to u/7». Hence if the values given 
in the following for the lifetime are changed in proportion to 
the change of mass, the energy distribution will be relatively 
insensitive to the uncertain mass. 
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Fic. 1. The energy distribution at sea-level for vertically incident mesotrons 

as a function of momentum (in units wc) plotted on a logarithmic scale (lower 

scale). The upper abscissa scale gives the corresponding vaiues of e«, the energy 

in units wc?. The circles and crosses give the experimental points. The calculated 

distribution for ro = 2.1 X 1078 sec. and 2.9 X 107* sec. are given by the full 

and dashed curves, respectively. The curve marked H was calculated by 

Hartree ® with ro = 2.7 X 1076 sec. and an initial distribution proportional to 
p?-’. All curves normalized to a maximum value of unity. 


The method adopted in carrying out the numerical 
integration in (12) is described in the appendix. The results 
obtained for the energy distribution are given in Fig. I. 
Both the experimental points and the calculated curves have 
been normalized to a maximum value of unity. The full 


1 For ro = 2.1 X 10°* sec., see below, the normalization constant Cp has 
the value 7.4 X 104 with this normalization. With this choice of the normaliza- 
tion constant the total intensity in the zenith is 32.8 mesotrons per cm.? per sec. 
per unit solid angle. The transition to other normalizations is obvious and needs 
no comment. 
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curve (to = 2.1 X 107° sec.) is seen to agree quite well with 
the experimental points of Blackett !* and of Jones,'” except 
at low energies. It is to be remembered that the data for 
small energies are very uncertain and subject to large sta- 
tistical errors due to the comparatively low intensity of slow 
particles. We may exclude the possibility that the dis- 
crepancy between observations and the calculated form of 
the distribution is due to an appreciable number of mesotrons 
which originate below the assumed position of the formation 
layer. For such mesotrons there would be relatively less 
attenuation due to decay and, since the decay affects the 
intensity of the slow particles to the greatest extent, the 
calculated distribution would then decrease less rapidly for 
low energies. This would increase the discrepancy rather 
than improve the agreement. While it is difficult to explain 
the deviation at low energies except in terms of the experi- 
mental uncertainty, the considerations of the following sec- 
tions, which are based on the form of the sea-level distribution, 
are to a large degree independent of this part of the spectrum. 

The theoretical curves have been calculated for several 
values of the lifetime and two curves are shown in Fig. 1: 
2.1 and 2.9 X 10-*sec. The normalization factors for various 
values of the lifetime differ considerably. However, the 
important point is that the shape of the distribution is not 
very sensitive to the lifetime and any rt) between the two 
values cited above would account for the observed results 
almost equally well. The insensitivity of the distribution 
function to the value of the lifetime may be regarded as 
unfortunate if it is desired to use these measurements to 
determine to. However, in view of the indirect nature of 
such a determination of the lifetime, this procedure could 
not be expected to yield a reliable value. Thus, accepting 
the calculations and measurements at their face value, a 
lifetime of as much as § X 10~-* sec. would appear to be ruled 
out. But the discrepancy between the calculated curve and 
the observations in this case is perhaps even within the 
uncertainties involved in both, and it is not possible to 
definitely exclude such a large value of to. Therefore, so far 


j 16 Pp, M. S. Blackett, Proc. Roy. Soc., A159, 1 (1937). 
17H. Jones, Rev. Mod. Phys., 11, 235 (1939). 
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as the energy distribution at sea-level is concerned, it seems 
safe to say that the observations are compatible with any 
value of the lifetime in the range reported in more direct 
experiments.” 

From the sea-level energy distribution one readily obtains 
the integral energy distribution, that is, the number of 
mesotrons with energy greater than e: 


N,(€) = | F(e, t,)de. (13) 
Here and in the following calculations the sea-level distribu- 


tion corresponding to a lifetime of 2.1 & 10~* sec. will be used. 
The integration (13) is then carried out numerically using the 


Fic. 2. The integral energy distribution at sea-level for vertically incident 
mesotrons. With the adopted normalization (see caption to Fig. 1), the total 
number of mesotrons per cm.? per sec., per unit solid angle at the zenith is 
N,(0) = 32.8. 


curve given in Fig. 1. The result is given in Fig. 2. The 
total intensity (within a constant factor) of vertically incident 
mesotrons at sea-level is thus obtained from this curve when 
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e is the minimum energy required to penetrate the walls of 
the counter or other device used as a detecting instrument. 


V. ABSORPTION BELOW SEA-LEVEL IN LIGHT MATERIALS. 


As was discussed above, the energy spectrum of the 
mesotrons in the formation layer can be determined from the 
results of absorption measurements of the particles under 
thick layers of dense materials such as rock. This conclusion 
is based on two facts. First, the energy necessary to pene- 
trate 30 meters of rock (about 100 meters water equivalent) 
or more is at least 3 X 10'°e.v. The energy at the formation 
layer is then about 10 per cent. greater. For such energies 
the attenuation due to decay in traversing the atmosphere is 
unimportant. The exponential factor (cf. Eq. 12a), whose 
departure from unity is a measure of the effect of decay, is 
0.92 for an energy at sea-level of 3 X I10!° e.v. The second 
basic fact is the following: in dense materials the effect of 
decay is negligible in all cases. Thus, only stopping by 
energy loss need be considered and the atmosphere could be 
regarded as an additional thickness of dense material whose 
thickness is given by the 10.3 meters water equivalent of the 
air layer. The absorption of the mesotrons is then directly 
determined by the primary energy spectrum since the number 
stopped in a thin layer dt at a certain depth ¢ is just the 
number of particles in the formation layer whose energy lies 
between E and E + dE, where E and E + dE are the energies 
necessary to just penetrate the layers of thickness ¢ and 
t + dt respectively. The primary spectrum is thus defined 
by the observed absorption and the given energy-range 
relationship. 

The conclusion that decay in dense materials is always 
negligible may appear surprising at first sight. It may be 
argued that the intensity of the slow particles emerging from 
a thick layer should be very much dependent on decay. 
Moreover, it may further be supposed that the energy below 
which decay is appreciable is greater the greater the thickness 
of the material. This would appear plausible since greater 
thickness means greater path length and more time spent in 
traversing the layer. Both of these arguments are misleading. 
It is true that the slower the particle the greater the attenua- 
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tion due to decay. However, when the energy is reduced 
much below the rest energy, the rate of energy loss increases 
so rapidly that the particles will be stopped before appreciable 
numbers disintegrate. Thus even for a final energy which 
vanishes, the decay factor is very little different from unity 
(cf. below).'8 Again, increasing the thickness of the material 
does very little to enhance the effect of decay since the 
increased thickness requires a greater energy at the top of the 
layer in order that the mesotron may emerge. Thus, one 
simply adds to the path a portion in which the velocity is 
increasingly larger and the added effect of decay becomes 
increasingly smaller. In the limit of infinite thickness and 
vanishing energy of emergence the decay factor for light 
materials such as rock becomes exp (— 4.4 107*/p) where p 
is the density of the material. This is virtually unity for all 
non-gaseous materials. 

Since only energy loss is of importance in dense materials, 
the energy distribution for vertically incident mesotrons at a 
depth ¢ cm. water below sea-level is 

Bie’ 

F(e, t)de = a F,(e')de, (14) 
where e’ is the energy of the mesotron at sea-level and F,(e’) 
is the energy distribution at sea-level for vertical incidence. 
The functions B refer, of course, to the dense material. For 
the comparison with the measurements to be considered in 
this section we are concerned with the energy loss in rock, 
whose heterogeneous chemical composition cannot be known 
exactly. However, the constituents of the absorbing material 
are almost entirely composed of light materials for which’ 


18 For materials of*uniform density the decay factor becomes exp (— D’), 
where 


D' = 1/(or«) [» dR/p(R). 


Here p is the density of the material and R’, R are the ranges at the top and 
bottom of the layer. The idea that decay should be important for very low 
energies amounts to the supposition that the integral D’ will diverge for vanishing 
energy (lower limit = 0). Actually the integral converges because the range 
vanishes with decreasing energy more rapidly than does the momentum. E.g., 
if for small energy we have R~ E*, then n > 1. The integrand of D’ then 
behaves as R~1/2" which gives a convergent result. 
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Z/A =~ 3, cf. Eq. 3. The dependence of B on the nature of 
the material is then contained in the average ionization 
energy J(Z) which enters only logarithmically. It is then 
permissible to assume for this ionization energy a value 
corresponding to an average Z. We take Z = 13 and 
I(Z) = 175 e.v. (Eq. 3a). The small error thereby incurred 
is of the same order as the uncertainty in the numerical 
factor in (3a). The same value of Z (and A = 2Z) is used 
in computing the radiation energy loss, Eq. 4. In this case 
the uncertainty in Z is of greater importance so far as the 
radiative rate of energy loss, B,, is concerned. On the other 
hand radiative losses play a very minor role in determining 
the absorption in light materials for energies less than 5 X 10"! 
e.v. which constitutes practically an upper energy limit for 
our considerations. 
The number of mesotrons transmitted through a thickness 
t is given by 
r nr / 
N(e, t) = { F(e, thde = mis) F,(e')de, (15) 
& 4 B(e) 7 
where ¢, is the energy required to penetrate the walls of the 
detecting instrument, counter telescope. For a given e; the 
absorption could be calculated by first determining the energy 
distribution from (14) for various thicknesses ¢ and then 
carrying out the integration (15) numerically for each thick- 
ness. However, a much simpler procedure is the following. 
First of all we note that with the usual experimental 
arrangements, €; ~ 0.1. Since the energy distribution below 
a layer of dense material will not differ greatly from that at 
sea-level so far as shape is concerned, we may conclude that 
the number of particles stopped by the counter walls, etc., is 
negligible and the number of transmitted particles may be 
written as 
; *B(s} 5 
N(t) = J, Be Fy (e')de (16) 
=-[ F(e’)de’ = N,(€0). 


0 


Here € is the energy at the top of the layer corresponding to 
zero energy at the bottom and JN,(e9) is the number of meso- 
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trons at sea-level with energy greater than €o, that is, the 
number with range equal to or greater than the thickness of 
the rock. Thus the absorption is directly obtained from the 
results of Fig. 2. 

Dividing by N,(0), the total number of particles at sea- 
level, we obtain the absorption curve normalized to unity at 
zero thickness. For all but the smaller thicknesses (¢ 5 40 m. 
water) the sea-level energy distribution for E > 10! e.v. 
becomes of importance. Since the observations do not extend 
beyond this energy, the absorption was first calculated by 
extrapolating the e?-* distribution. The result is shown in 
Fig. 3 by the curve drawn in full up to 280 m. water equivalent 
and by the dashed curve for greater thicknesses. The dis- 
agreement between this calculated curve and the experi- 
mental points due to Wilson '* indicates a larger absorption 
coefficient at greater thicknesses. This was, in fact, originally 
pointed out by Wilson. While some attempt has been made 
to explain the rather abrupt change in absorption near 280 m. 
water in terms of some additional process of absorption or by 
postulating some change in the nature of the hard component, 
it does not seem possible at this time to arrive at a consistent 
interpretation of the observations in this manner. Instead, 
one may assume that the extrapolation of the primary energy 
distribution as given by (8) is unjustified, as it certainly is, 
and the behavior of the spectrum in the formation layer for 
energies above 8 X 10!° e.v., which is necessary for the 
penetration of 280 m. water, is deduced from the observed 
absorption for the larger rock thicknesses. The full curve in 
Fig. 3 gives the calculated absorption when it is assumed that 
for e > 800, the energy distribution in the formation layer 
is proportional to e~*-*. When the multiplicative constant is 
chosen so as to make the distribution continuous, the primary 
energy spectrum is as follows: 


Fy(e)\de = 7.4 X 104e-? “de, e < 800, (17) 
2.1 X 107%e—*-4de, e > 800. 


While the observed absorption may in this way be used to 
determine the primary energy distribution up to very large 


1 V. Wilson, Phys. Rev., 53, 337 (1938). 
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energies, it must be admitted that the adopted procedure is 
decidedly ad hoc. Thus the increase in absorption may be 
only partly due to a modification in the primary energy dis- 
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Fic. 3. The absorption in light materials for vertically incident mesotrons. 
The unit of thickness is 100 g./em.2 = I meter water. The circles are the experi- 
mental points of Wilson !* and the full curve is the calculated absorption with the 
primary energy distribution (17). The dashed curve results when the distribution 
(8) is extrapolated to high energies. Both calculated and experimental intensities 
normalized to unity at sea-level. 


tribution and some as yet unspecified mechanism of additional 
absorption may still be involved. This question is best 
decided when further independent evidence concerning the 
primary spectrum is at hand. 
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VI. ANGULAR DISTRIBUTION AT SEA-LEVEL. 


The energy distribution at sea-level for vertically incident 
mesotrons and the absorption in rock essentially serve to fix 
certain parameters; mainly, the primary energy distribution 
above the energy required to reach sea-level. Therefore, 
from the comparison of calculated results and those observed 
in experiments of these two types it is hardly permissible to 
draw conclusions as to the validity of the postulated absorp- 
tion mechanisms. The validity of the assumed energy dis- 
tribution, especially at higher energies than are important for 
vertical incidence at sea-level, can be tested by a comparison 
between observed and theoretical angular distributions at sea- 
level. For the inclined rays the path through the atmosphere, 
measured both in usual geometric units as well as mass units, 
is increased by a factor sec 3, where # is the zenith angle. 
Therefore the decay and ionization loss are enhanced corre- 
sponding to the longer path. The number of mesotrons then 
depends more directly on the intensity in the formation layer 
at the higher energies. 

Before the angular distribution at sea-level may be deter- 
mined some decision in regard to the angular distribution in 
the formation layer must be made. In the absence of any 
information about the primary distribution in angle a calcu- 
lation such as the following would only serve to obtain this 
information empirically. Fortunately, such a priori know]l- 
edge of the initial angular distribution is not entirely lacking. 
Although the precise nature of the process by which mesotrons 
are produced is not known, it is fairly certain that the collisions 
responsible for such production are due to short range forces. 
It would follow then that the mesotrons are produced with an 
isotropic distribution in angle. We assume that this is the 
case. We also assume that the primary distribution in energy 
is independent of the initial direction of motion so that (17) 
applies to all mesotrons in the formation layer. 

The calculation of the angular distribution involves first 
the determination of the sea-level energy distribution for each 
angle of incidence. The procedure is exactly as described in 
Section IV (and the appendix) except that in the decay factor, 
D is multiplied by sec # and the energy loss is computed for a 
thickness (t, — to) sec 8. The total intensity for each angle 


32 M. E. Rose. (J. F. 1 


is obtained by integration over the corresponding energy 
distribution as in (13). The angular distribution, normalized 
to unity for vertical incidence, is then 


=,” (18) 


70 


where now the second argument in F, indicates the zenith 
angle. 


' rr) r) 7 4 S : 
cos § 

Fic. 4. The angular distribution for mesotrons at sea-level. The curve is 
the calculated distribution. The designation of the experimental points is as 
follows: O, Bernardini, X Cocconi and Tongiorgi, 0 Bernardini and Bocciarelli, 
A Greisen. Both calculated and observed results are normalized to unity at 
zero zenith angle. 


The computed angular distribution is shown in Fig. 4 
together with the experimental points obtained from the 
observations due to Bernardini,?° Cocconi and Tongiorgi,?! 
Bernardini and Bociarelli *? and Greisen.?* It may be seen 


20 G. Bernardini, Nat., 129, 578 (1932). 

21G. Cocconi and V. Tongiorgi, Phys. Rev., 57, 1180 (1940). 
22 G. Bernardini and D. Bocciarelli, Ric. Sct., 1, 3 (1935). 
23K. Greisen, Phys. Rev., 61, 212 (1942). 
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that the agreement between the calculated curve and the 
data is as good as the agreement between the results of the 
various observers. In the absence of greater consistency. in 
the experimental results one may only conclude that the 
theory is in accord with the measurements within the experi- 
mental error. 

It may be noted that the theoretical angular distribution 
corresponds only roughly to a cos? # variation which is the 
form usually assumed. The distribution found conforms 
quite closely to the cos? # for intermediate angles (30°-40°) 
but decreases more rapidly for larger angles where it behaves 
as cos?*? #, 

A further assumption may be verified by the above com- 
parison. It has been supposed that the effect of scattering 
may be neglected. If an appreciable amount of scattering 
took place the angular distribution as observed would sys- 
tematically deviate from the calculated one in the direction 
of greater isotropy. Since no systematic discrepancy of this 
kind is evident, the assumption of the negligible importance 
of scattering would seem to be justified. 

From the calculated angular distribution one may obtain 
the integral angular distribution which may be of interest for 
certain experimental arrangements in which a limited solid 
angle is subtended by the apparatus. Thus, 


vo 
(do) = an | UT Io) sin ddd (19) 
0 


gives the total intensity at sea-level, in units of the vertical 
intensity, as determined by a detecting instrument which sub- 
tends the solid angle 27(1 — cos 8). The result is shown 
in Fig. 5 where ¢ is given as a function of cos #). No meas- 
urements are available for comparison. 


VII. ABSORPTION IN LEAD AND OTHER MATERIALS. 


The importance of the assumption of the absorption 
mechanism in dense materials, that is, energy loss by ordinary 
electromagnetic collisions, has been emphasized in the fore- 
going. The experiments designed to furnish some confirma- 
tion of this assumption are largely concerned with absolute 
measurements of the absorption coefficient in a given absorb- 
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ing material and with the comparison of the absorption in 
different materials. 

The experiments in which we are interested are those in- 
volving the use of comparatively thin absorbers, that is, 
layers sufficiently thin so that the energy loss in the material 
is small compared to the energy at which radiative losses are 
of importance. For lead this restricts the consideration to 
thicknesses small compared to 10% meters water equivalent 
(100 meters Pb), and for absorber of atomic number 13 the 
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Fic. 5. The integral angular distribution at sea-level, see Eq. 19. 


pertinent thicknesses must be small compared to 6,000 meters 
water equivalent. These thicknesses are well above those 
involved in the experiments. As a consequence, the intensity 
of the rays emerging from the absorber is mainly due to 
mesotrons which have suffered only ionization loss. 

We consider first the comparative absorption data which 
evidently provides a measure of the dependence of ionization 
loss on atomic number. It is customary to present the data 
in the form of equivalent absorber thicknesses where the 
implication is equal absorber efficiency. More precisely, the 


: 


July, 1943-1 HARD COMPONENT OF Cosmic RADIATION. 35 


observed value of the ratio of absorption coefficients (measured 
in cm.~') in the two materials is compared with the calculated 
ratio of atomic electron densities, since according to (3) the 
rate of energy loss per cm. is, aside from the insensitive 
logarithmic factor, proportional to pZ/A, or the number of 
electrons per unit volume.** These ratios which are supposed 
to be identical within the experimental error, are identified 
with the inverse ratio of equivalent absorber thicknesses. An 
examination of the data, however, reveals what appears to 
be a discrepancy since materials of low atomic number absorb 
more strongly than is expected. 

Against this method of analyzing the data two objections 
may be raised. First, the comparison of absorption coeff- 
cients and electron densities implies that the former is equal 
to the rate of energy loss, within a factor independent of the 
material, or very nearly so. Actually the absorption coeff- 
cient depends on the energy distribution of the mesotrons 
incident on the layer and, as will be shown below, for different 
materials the pertinent portion of the spectrum will be dif- 
ferent. It must be admitted that this fact is due to the non- 
constancy of the rate of energy loss but, while the energy loss 
varies slowly, the effect on the energy distribution and on the 
absorption coefficient due to this slow variation is in this case 
fairly important. The second objection is concerned with the 
definition of equivalent absorber thicknesses. The most use- 
ful definition is one based directly on what is measured. 
Thus, we define equivalent thicknesses of two materials as 
those thicknesses which will reduce the intensity of mesotrons 
at a given place by the same amount.” This definition would 
be equivalent to that used above according to which the 
thickness ratio is the reciprocal of the ratio of absorption 
coefficients, provided these coefficients were independent of 
thickness. Actually this is not the case. 

In the following we consider mesotrons incident in the 
direction of the zenith and at first we shall consider the ab- 


*4E. G. Steinke and H. Tielsch, Zeits. f. Physik, 84, 425 (1933); H. Tielsch, 
ibid., 92, 589 (1934); A. Sittkus, zbid., 108, 421 (1938). 

25 See, e.g., M. A. Pomerantz and T. H. Johnson, Phys. Rev., 59, 143 (1941). 
In the definition given it is essential to specify the incident energy distribution 
in terms of the amount and nature of material between the absorber and sea-level. 
Obviously, the angular distribution of the particles is of importance also. 
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sorption at sea-level. A total of N, mesotrons are incident 
on the absorbing plate of thickness ¢ and the number trans- 
mitted will be denoted by N,(t,) or N2(t,) to distinguish be- 
tween two different materials. It will be convenient to 
introduce what may be called the attenuation coefficient k 


according to 


1 dN; 
k; container = S , (20) 
N, dt; 
The absorption coefficient x is then given by 
d N, 
;( = ——log N; = —&,. 21 
: * dig N, vy 


For vanishing thickness the two coefficients become identical. 
The number of transmitted mesotrons is, as before, 


Ni = [~ Fiode, (22) 


where e¢,; is the minimum energy required to penetrate the 
thickness ¢; (Eq. 5a). Then we have from (20) 


k; = B,(e,;) F.(€;)/N; (23) 


"ty —1 
kK, = k. (1 -{ eat ) (24) 


The dependence of the attenuation and absorption coeffi- 
cients on thickness is of interest. Since the energy e; is 
rougly proportional to the thickness it follows from the form 
of the sea-level distribution (Fig. 1), that at large thickness 
the attenuation coefficient k; decreases rapidly (about as ¢-7). 
At small thicknesses the factor F,(e;) increases while B(e;) 
decreases with increasing t. However, as may be seen from 
(12) and (3), F, varies more rapidly and over a wider range. 
Thus for small ¢ the attenuation coefficient increases and a 
maximum is attained for some intermediate thickness. This 
behavior is depicted in Fig. 6 where k for Pb and a light 
material (see Section V) is given as a function of thickness. 
Due to the decrease with thickness of the factor in the 
denominator of (24) the absorption coefficient x, also given in 


and from (23) 
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Fig. 6, decreases less rapidly than does k and the maximum 
is considerably flatter. Nevertheless, it is apparent that the 
assumption of constant absorption or attenuation coefficient 
is an extremely crude one. 

According to the definition given above, the equivalent 
thicknesses, ¢; and f., of two materials are related by 


ty "ta 
| kidt = i] kedt. (25) 
vo 0 
The thickness of light material ¢;, equivalent to a given thick- 


ness of Pb, fp, is at once obtained by (25) using the attenua- 
tion coefficients given in Fig. 6. The results, when ¢; is 
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Fic. 6. The attenuation coefficient k and the absorption coefficient « 
(Eqs. 20 and 21), for Pb and a light material (Z2/A ~ 3, lo << Z< 15). The 
ordinate gives the coefficients in units of 1074 cm.?/g. while the abscissa gives 
the layer thickness in units of 100 g./cm.? or meters water equivalent. 


expressed in g./cm.?, are shown in Fig. 7. It is remarkable 
that despite a non-constant attenuation coefficient the ratio 
of equivalent thicknesses is constant (within the accuracy of 
calculation and certainly within the error of experiment). 
For all practical purposes the curve given in Fig. 7 may be 
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expressed by 

ty = 0.75tp». (26) 
This may be compared with the ratio of equivalent masses 
based on the ratio of absorption coefficients as given by 
Tielsch: 24 0.78 for Fe/Pb, 0.78 for Al/Pb, 0.55 for H:O/Pb 
and 0.65 for C/Pb. 
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Fic. 7. The thickness of light material t, equivalent in stopping power 
to a thickness tp, of lead. Both ft; and tpp are given in units of 100 g./cm.?. 


The average slope of the curve is 0.75. 


Having obtained the attenuation coefficient we may com- 
pare the calculated absorption with the observations for a 
given material. The relative intensity transmitted is 

N "t 

—=][— kdt. (27) 

N 8 e700 
Fig. 8 shows the theoretical absorption curve for Pb at sea- 
level (vertical incidence) and the experimental points taken 
from the measurements of Nielsen and Morgan,?* Sittkus,”* 


26W. M. Nielsen and K. Z. Morgan, Phys. Rev., 54, 425 (1938). 
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Rossi 2? (for thicknesses up to 1 meter of Pb) and of 
Pomerantz 28 for larger thicknesses. While the experimental 
results usually indicate a linear absorption (constant &) it 
should be pointed out that over the range of observation 
(10-100 cm. Pb) the calculated curve is straight within the 
experimental uncertainty (of order 10 per cent.). However, 
a discrepancy of a somewhat more serious nature lies in the 
fact that the observed absorption coefficient is, in almost 
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Fic. 8. The absorption of mesotrons in lead for vertical incidence at sea- 
level. The curve is calculated and the designation of the experimental points 
is as follows: O Sittkus, @ Nielsen and Morgan, A Rossi, 0 Pomerantz. The 
calculated curve is normalized to unity at zero Pb thickness as are the results 
of Sittkus and of Nielsen-Morgan. The data of Rossi and of Pomerantz are 
each fitted to the calculated curve at one point. 


every case, larger than the calculated one. For the first 
meter of Pb the calculations give an average absorption 
coefficient of 3.9 X 10-* cm.-! Pb. The data of Nielsen and 
Morgan give 4.5 X 10-* cm.~!, that of Sittkus 4.9 X 107% 
cm.~' and Rossi’s data when fitted at the larger thickness 
give 5.0 X 107% cm.~! Pb. These deviations would seem to 
lie outside the limit of error of the measurements. At larger 
thicknesses the data of Pomerantz, when fitted at the point 
indicated in the figure, gives an absorption coefficient of 


27 B. Rossi, Zeits. f. Physik, 82, 151 (1933). 
*8 M. A. Pomerantz, Phys. Rev., 57, 3 (1940). 
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2.5 X 10-* cm.~! Pb compared to the calculated coefficient, 
averaged over the range of thicknesses concerned, equal to 
2.7 X 10-* cm.-!. The comparison in this case is more 
favorable and it is also satisfying that the theory accounts for 
the diminution of absorption with increasing thickness. 

The absorption coefficient in Pb has also been measured 
at various depths below sea-level. Thus, for mesotrons 
incident in the zenith, the Pb absorption coefficient at 20 
meters water below sea-level is 3.0 X 10-* cm.~! (measure- 
ments from 10-100 cm. Pb).2*° At 60 meters water below 
sea-level the absorption coefficient is 2.3 X 107% cm.~! accord- 
ing to the measurements of Nielsen and Morgan.?* This 
decrease of the absorption coefficient with depth can also be 
understood from the theory. 

Suppose the depth of rock or other material at which the 
Pb absorption measurements are made is ¢; and the incident 
intensity is V,(t;). The thickness of Pb is ¢, and the intensity 
at this thickness is N2(t:). Let e. be the minimum energy 
required to penetrate the Pb alone and e¢,; the corresponding 
energy at sea-level; that is, the energy required to just pene- 
trate the rock and Pb together. Then if B, and B, denote 
the rate of energy loss in rock and Pb respectively, we find 


Bi (e:) F,(€1) ; 
B,(€2) Ni(t:) 


If the thickness ¢; is much greater than f2, as is approximately 
the case, we can write 


k2(t:) = B2(e2) (28) 


d N, 
ky = xo(t;) ~ — 0.8—log—, 28a) 
2 o( 1) dt, BN, ( 
where we have used B.(e2)/B,(e,) ~ 0.8. Using the absorp- 
tion measurements of Wilson !* which can be expressed in the 
form 
N,i/N, = (1 + t/t) -’, (28b) 
where y’ ~ 1.8 (t; = 280 meters water) and fy is the atmos- 
pheric thickness, we find 


Se “0x 
de 2M -) 


29 V. Wilson, Phys. Rev., 53, 908 (1938). 
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The decrease of Pb absorption coefficient with depth is thus 
directly related to the decrease of the absorption coefficient 
in rock. The result (28c) is too rough to allow a significant 
quantitative comparison with the observed values. 

The results obtained in this section are of interest in still 
another connection. Because of the prime importance of 
obtaining direct evidence for decay in the form of cloud 
chamber photographs of the disintegrating mesotron, it is 
essential to examine the way in which the number of slow 
particles depends on the thickness of absorber placed above 
the chamber. As is evident from the definition (20), the 
number of slow mesotrons is proportional to the product of 
the attenuation coefficient & and the range of these mesotrons. 
Since the range depends only on the definition of a ‘‘slow”’ 
mesotron (e.g. ionization equal to twice the normal value 
associated with fast particles), the dependence of the in- 
tensity on thickness is given by the attenuation factor. The 
number of slow mesotrons beneath a given thickness of lead is 
therefore proportional to kp,.*° From Fig. 6 it is seen that 
the addition of Pb up to 50 cm results in a rapid increase in 
the number of slow mesotrons. However, for larger thick- 
nesses the number of such mesotrons decreases so that the 
thickness given is an optimum value. While the existence of 
an optimum thickness is, perhaps, not obvious it is to be 
understood very simply: The number of slow mesotrons is 
essentially the number in a narrow energy band at the top 
of the absorber which have just enough energy to penetrate 
the layer. For increasing thickness this band is selected from 
a portion of the sea-level distribution which lies at ever- 
increasing energy and the number of mesotrons in this band 
(which is of practically constant width) varies with thickness 
essentially as the sea-level distribution varies with energy. 
The behavior here predicted is in qualitative agreement with 
preliminary observations obtained in this laboratory.*! 

30 Tf account is taken of the fact that there is, in addition to the lead, usually 
some amount of light material (wood, concrete, etc.) above the chamber, the 
number of slow particles becomes proportional to the attenuation coefficient as 
given by (28). Since for large lead thicknesses this is proportional to the value 
of k for lead alone, the conclusions are qualitatively the same in this case as 


those given below. 
3} R. P. Shutt, unpublished. 
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Vill. ABSORPTION AT ELEVATIONS ABOVE SEA-LEVEL. 


The procedure utilized to calculate the energy distribution 
and intensity of mesotrons at sea-level may be applied to 
obtain the intensity in the atmosphere as a function of eleva- 
tion above sea-level. Retaining the assumption of a produc- 
tion layer at 1.0 meter water equivalent below the top of the 
homogeneous atmosphere, the results of the calculation should 
show how serious this assumption is. At the same time an 
indication of the importance of mesotron production at lower 
elevations may be obtained. 

The energy distribution for vertically incident mesotrons 
has been calculated for the elevation of Mt. Evans, 6.2 meters 
water equivalent below the top of the homogeneous atmos- 
phere. The procedure used is as given in Section IV and the 
appendix. From the energy distribution the intensity relative 
to that at sea-level is readily determined. With a lifetime of 
2.1 X 10-*sec. we find for the ratio of intensities at Mt. Evans 
to that at sea-level the value 1.98. The observed value * 
is 2.25. The difference is significant and cannot be explained 
by errors of calculation or experiment. 

It is essential to note that the calculated ratio is based 
on the extrapolation of the e-?-* energy distribution in the 
formation layer down to ¢ = 12.8 which is the minimum 
energy required for transmission through the 5.2 meters water 
from the assumed position of the formation layer to the Mt. 
Evans elevation. Of course, no evidence can be cited in 
support of this extrapolation and it is moreover rendered 
uncertain by the necessity of some modification of the energy 
distribution due to latitude effect. This arises from the 
exclusion of slow primary particles by the action of the mag- 
netic field of the earth.** If anything, it is more plausible to 
assume that the extrapolation of the distribution (8) over- 
estimates the number of slower mesotrons and therefore gives 
too large a value for the relative intensities at elevations 
above sea-level. Hence the source of the discrepancy noted 
above is to be sought elsewhere. 


8B. Rossi, N. Hilberry and J. B. Hoag, Phys. Rev., 57, 461 (1940). See 
also reference 5. 


38 T, H. Johnson, reference 2. 
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By far the most plausible reason for the low value of the 
calculated intensity ratio is the possibility that mesotron 
production occurs much below the assumed position of the 
formation layer. It is readily seen that if this is the case 
the increase of intensity for the elevation of Mt. Evans, for 
instance, would be greater than that for sea-level. Conse- 
quently, the calculated ratio would be increased as the data 
demands. If the discrepancy between calculated and ob- 
served intensity ratios is to be entirely attributed to pro- 
duction of mesotrons at lower elevations, then the results 
presented above indicate appreciable mesotron production at 
elevations comparable with the 4,300 meters at Mt. Evans. 
Essentially the same arguments have been used by Greisen ** 
to arrive at the same conclusion. Finally, this conclusion is 
in accord with the observations of Regener and Regener and 
Rossi. *4 

APPENDIX. 

The decay integral, upon which the energy distribution 

depends, can be written in the form (cf. Eq. 12a) 


ko dR I 

D{R.) = J. 7 — ROR)’ paca 
where » = R +t is a constant for the integration (Eq. 12b). 
Obviously, the lower limit R, may refer to the range at other 
elevations than sea-level. The integral can be evaluated 
numerically for each energy by some procedure as the method 
of trapezoids. However, since we are interested in the values 
of D over a wide energy range (from ¢ = 0.1 to about 500) it 
is necessary to compute the integral for many points in the 
energy scale. Again, the integrand varies rapidly, especially 
in the region where it has its largest values, and it is necessary 
to resort to a rather fine division of the range scale in order 
to obtain a reasonably accurate value for D. Finally, the 
results of a straightforward integration would apply for only 
one elevation and it would be necessary to repeat the process 
for other elevations. A method which avoids many of these 
objectionable features has been devised and is described 
below. 


4 V. H. Regener, Phys. Rev., 61, 105 (1942); V. H. Regener and B. Rossi, 
Phys. Rev., §8, 837 (1940). 
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It is clearly advantageous to carry out the integration in 
an analytical way so far as this is possible, and for this it is 
necessary to find an analytical representation of the range- 
momentum relationship. It is also essential that this repre- 
sentation be valid especially at low energies and that it be 
such as to permit the integration to be carried out in a not- 
overcomplicated manner. For this purpose it was found most 
suitable to introduce functions which vary slowly. First, we 
note that the momentum is roughly proportional to the range 
and therefore it is convenient to introduce R/p in (1A). This 
also suggests the introduction of log R instead of R as a 
variable of integration. This has the advantage that R/p is 
then a slowly varying function of the integration variable 
over the entire range of integration. We may therefore 
replace the integrand, or R/p. by a series of straight line 
segments. 


Rip =a,+6,logR, Riris=RER, (2A) 


where the R, are selected points chosen so that the straight 
line segments give a good representation of the function R/p. 
In many cases only one or two such straight lines are neces- 
sary. It will be noted that once the representation (2A) is 
fixed it applies for all elevations. 

We now introduce the function G(R) according to 


dR! 
ead Rp(R) : 


where R,, is smaller than the minimum value of the range for 
which the decay integral will be evaluated: R,, = 1 g./cm.? is 
a suitable value. Evidently, we have 


D(R,) = G(Ro) — G(R.), (4A) 


so that it is sufficient to evaluate G(R). In view of (2A) 
this function will be expressed as a series of terms of the form 


oe Sexi... See : 
GAR) = J any + log R), (5A) 


G(R) = (3A) 


where the first R, is R, and the last is Ry. The integration 
(5A) is readily carried out. We introduce y as an integration 
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variable where 
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J(y) = si *® tog 2, (7A) 
7 . 


J(y) cannot be evaluated in closed form but is easily evaluated 
in terms of rapidly convergent series. We have 
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The method of evaluating the decay integral D now con- 
sists in first computing R/p as a function of log R. Then, 
having chosen a suitable representation of the form (2A) the 
values of a,, 6, are fixed and for a given point of observation 
(t), the values of y, are prescribed. From the curve of J(y) 
plotted against y (Eq. 7’A) the values of J(y,) can be read 
and G,(R) computed from (6A). In this way for each value 
of ¢ it is not difficult to find the function G(R) and hence the 
decay integral D at sea-level and at other elevations. 
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Bomb Shelters in Hudson Palisades Proposed.—( Compressed 
Air Magazine, Vol. 48, No. 2.) Mr. George C. Atwell, a prominent 
New York contractor, has made a suggestion for utilizing the 
Palisades, bordering the Hudson River directly across from New 
York City, for bomb shelters. The scheme embodies measures 
similar to those taken at Chungking, Dover, the Island of Malta, 
and Gibraltar, where underground galleries offer the citizenry and 
even some essential industries refuge from bombings. The topog- 
raphy and geological features of the Palisades are highly favorable 
to the plan. The cliffs rise sheer from the west side of the river to 
a height of 200 feet and more, extend for several miles along the 
stream and have a horizontal thickness of approximately 13 miles 
throughout that section. They are composed of doleritic rocks 
which were intruded from subterranean lava reservoirs and spread 
out laterally between the horizontal beds of sedimentary deposits 
that formed the surficial covering. During the intervening cen- 
turies—geologists estimate the age of the Palisades at 150,000,000 
years—the softer sedimentaries were worn away, exposing the 
igneous rocks. Along their eastern edge a deep gorge was cut by 
the Hudson River, and this was followed by a period of general 
subsidence, accompanied by filling of the stream bed with silt to 
its present level. Mr. Atwell, who says the plan was suggested to 
him by John C. Evans, chief engineer of the Port of New York 
Authority, proposes to drive a tunnel, 200 feet wide and 100 feet 
high, through the Palisades at a point just above the high water 
level of the river and emerging in the meadows of New Jersey to 
the westward. Within the depths of the sheltering rock he would 
enlarge this passageway to provide immense, pillar-supported rooms 
to serve as airplane hangars, depositories for national art treasures, 
sites for vital manufacturing plants, and havens during air raids 
for a large number of people living in the environs. All this, 
Mr. Atwell contends, could be accomplished with equipment and 
materials that are plentiful. As the rock, which is already pene- 
trated by several transportation tunnels, is self supporting, no 


buttressing steel would be required. 
KR. Hi. 0. 


A QUANTITATIVE EVALUATION OF THE INFLUENCE OF 
THE LITHOSPHERE ON THE ANOMALIES OF GRAVITY.* 


BY 
ROSS GUNN, 


Technical Adviser, U. S. Naval Research Laboratory, Washington, D. C 


ABSTRACT. 

A quantitative investigation is made of the influence of a strong elastic 
lithosphere on departures from isostasy. Isostasy is shown to be a special case 
of a more general principle that may be called the principle of isobaric equilibrium. 
This principle which expresses the condition for vertical equilibrium is formulated 
quantitatively both in terms of the resultant vertical stresses in the lithosphere 
and in terms of measurable gravity anomalies. Employing the principle, a stress 
function may be determined that permits an estimate of the fraction of the 
uniform superposed load carried by the lithosphere. By applying this method, 
it is found that for the central region of a block to be in substantial isostatic 
adjustment (96 per cent.), the linear dimension of the load must approximate 
330 Km., a value more than six times that often assumed. The basic constant 
characterizing the elastic deformations of the lithosphere is determined by appeal 
to geological data from deltas and regions near a thrust fault. The values so 
obtained suggest a value of 50 + 15 Km. for the effective thickness of the litho- 
sphere. The deformed figures of the lithosphere for certain special types of load 
have been worked out together with the stress function that measures both the 
gravity anomaly and the resultant vertical stress. 


Modern ideas concerning the equilibrium of the lithosphere 
were introduced about 1855 when J. H. Pratt prepared, as a 
result of his studies in India, an extensive report on the figure 
of the earth. Measurements were made between two places 
in northern India separated by an angle of 5° 23’ 42.294” as 
determined by a careful triangulation. Observations at the 
same stations by astronomical means led to a value which 
differed from that just mentioned by 5.236 seconds of arc. 
It was at first suggested that this discrepancy was due to the 
gravitational attraction of the Himalaya Mountains some 
60 miles north of the point of observation, but Pratt showed 
that if the Himalayas had a density of 2.75, the discrepancy 

* The opinions or assertions contained herein are the private ones of the 


writer and are not to be construed as official or reflecting the views of the Navy 
Department or the naval service at large. 
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should have been three times that observed. A year or so 
later George B. Airy discussed the discrepancy and advanced 
the idea that an unusually high mountain was analogous to a 
‘“‘log whose upper surface floats much higher than the upper 
surfaces of others.’”’ The original ideas proposed by Airy 
and Pratt have been examined in the light of numerous 
observational facts by many scientists and these have estab- 
lished the principle of isostasy, but only as applied to very 
large sections of the earth’s crust. Wm. Bowie! writes, 
‘According to the principle of isostasy, if the earth’s crust 
were cut into prisms of the same cross-section by imaginary 
vertical planes, the prisms would have the same mass if the 
isostatic conditions were perfect.’’ The principle of isostasy 
as developed within the last 30 years has thus been built upon 
the hypothesis of constant mass per unit area in these prisms 
which extend downward to some selected level where the 
strength of the crust in shear becomes quite negligible. 

Asa general rule, the difference between the calculated and 
the observed values of the acceleration due to gravity or the 
isostatic anomalies in undeformed and level country or over 
most of the ocean are satisfactorily low. However, in recent 
years as methods for correcting for elevation, density, and 
topography have become more accurate, and detailed data 
have become available, it has become increasingly clear that 
many regions of the earth are definitely not in isostatic ad- 
justment. Misleading statements have frequently appeared 
in the literature of geology suggesting that these isostatic 
anomalies imply an unstable geological structure even though 
observational evidence is available showing that in such re- 
gions no actual geological activity is in progress. The basic 
concepts of isostasy, therefore, need review and modification 
so that known departures can be interpreted quantitatively 
in terms of fundamental physical laws. R.A. Daly ? has well 
summarized the difficulties with present ideas concerning the 
term ‘‘isostasy”’ by writing ‘‘. . . its acceptability and its real 
value demand a certain vagueness in its definition.’’ ‘‘Al- 
though phrased in terms of lithosphere and asthenosphere, an 


1Wm. Bowie, Bull. No. 78, Nat. Res. Council, p. 103. 
2R. A. Daly, “Strength and Structure of the Earth.” Prentice-Hall, 1940, 


p. 340. 
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exact statement of thickness for either earth-shell is now out 
of the question.”” In this paper we will attempt to remove the 
‘‘demanded vagueness” and answer, in part at least, funda- 
mental questions regarding the strength and thickness of the 
lithosphere. It will be shown that in addition to isostatic 
forces, the vertical superposed stress in the lithosphere at any 
selected point must be considered and specified before any 
statement can be made as to the stability of any section of a 
continent or ocean bottom. 


IDEAL ISOSTASY. 


Consider a new born earth made up of rock material at 
such a temperature that its strength in shear is negligible. 
After circulation and mixing have acted for a few million 
years, such an earth will, except for a thin surface shell, 
approach the hydrostatic state characterized by: (a) a 
steadily increasing density toward the center; (0) the coinci- 
dence of surfaces of equal density, equal pressure and equal 
gravitational potential; and (c) a thin and probably floating 
surface layer of light rock. 

Suppose that great prisms of continental rock are ulti- 
mately formed at the surface by selective crystallization or 
some other means and that these aggregates are systematically 
less dense than the underlying weak molten rock. Obviously 
they will be supported and, according to Archimedes’ principle, 
will float in such a way that the weight of each prism will be 
exactly equal to the weight of the displaced underlying 
magma. Consider a number of prisms supported in the above 
manner, all of different but comparable density. Assume 
that each prism is self supported and no friction whatever 
exists between the lateral faces. Then if a reference prism is 
selected and if the acceleration due to gravity is measured on 
it, making suitable corrections for the distance of the ob- 
serving station above some reference level, and for the figure 
of the earth, it will be found that the measured values on all 
prisms will not appreciably differ. This is a consequence of 
the laws of gravitational attraction and potential theory which 
show that to a sufficient approximation the vertical accelera- 
tion due to a horizontal plane plate of considerable extent is 
given by 


Ag = 277M, (1) 
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where Ag is the change in acceleration or the gravitational de- 
fect due to the plate, M is the mass per unit area of the plate, 
and y¥ is the gravitational constant. If the plate is not exten- 
sive, corrections must be applied, but Eq. (1) will represent the 
principal term. Since the mass per unit area of the floating 
block is on the average exactly equal to the mass per unit area 
of the displaced liquid, and since by Eq. (1) the gravitational 
defect depends on this quantity, it is clear that different thick- 
nesses or densities of a floating block will not produce impor- 
tant changes in the measured value of the acceleration due to 
gravity. Because such a system of freely floating blocks will 
show a constant gravitational acceleration at some selected 
level, we may describe the state as one of ideal isostasy. 


DEPARTURES FROM IDEAL ISOSTASY. 


Suppose the separate prisms considered in describing ideal 
isostasy are brought into contact. Frictional and cohesive 
forces between the vertical faces of the prisms now prevent 
the individual adjustment of each. For example, select a 
central prism and assume that sediments are piled on it alone. 
If this prism were supported only by the weak asthenosphere, 
it would sink in the substratum and displace a weight of 
material exactly equal to that added. Due, however, to the 
assumed forces of friction acting on the prism, this cannot 
take place. The applied load attempts to force the prism 
down but now, because it is mechanically associated with the 
neighboring prisms, they are forced down while the original 
prism is supported by shearing stresses to exactly the same 
extent. Thus under the selected prism not enough mass is 
forced out laterally to support it from the bottom, and part 
of the support is supplied by vertical shearing stresses ap- 
pearing on the vertical faces of the prism. Conversely, as a 
result of such stresses, considerable material is forced out 
from under the adjacent prisms even though no load of 
sediments is applied to them. Of course, considering the 
aggregate of prisms, just enough material is displaced to 
carry the superposed load of sediments on the original prism 
and a condition of regional isostasy is achieved. 

The deficiency or excess of mass displaced from under any 
prism is readily calculated in terms of the applied force or the 
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vertical stress. For simplicity, consider a prism of unit 
cross-section supported ideally by the substratum. Suppose 
this prism is pushed down by an added load of sediments, by 
some arbitrary vertical force acting on unit area, or by 
shearing stress due to the finite strength of the lithosphere. 
This force per unit area is the vertical resultant stress SS. 
This stress, acting on a freely moving prism, will displace 
substratum material of density d until a restoring stress equal 
to the weight of the displaced liquid restores equilibrium, or 


S = M’g, (2) 


where M’ is the mass per unit area of the displaced material. 
It is clear, therefore, that a downward or positive resultant 
stress from any source whatever results in a deficiency of 
mass under the prism, while an upward stress results in an 
excess of mass. If this change in mass is assumed to be a 
flat sheet, the principal term of the accompanying gravita- 
tional defect Ag by Eqs. (1) and (2) is 


ghg = — 2ry75S, (3) 


where y is the gravitational constant and g is the acceleration 
due to gravity in the region where the displaced substratum 
is forced out. Attention is directed to the fact that, due to 
the frictional forces or cohesive forces acting on the vertical 
faces of the prisms, these prisms are just as much in equi- 
librium as if they were supported ideally from the bottom. 
It follows from the above discussion that if the resultant 
vertical stress distribution can be determined for any arbitrary 
load on the lithosphere either by measurement or by calcu- 
lation, then it is possible to evaluate the gravity anomalies 
associated therewith. Conversely, a knowledge of the space 
distribution of gravity anomaly and mass permits an estimate 
of the vertical superposed stresses and the resulting deforma- 
tions of the lithosphere. Subsequent paragraphs consider 
two interesting cases in which the distribution of the vertical 
stress or the stress function can be determined. 


ISOBARY, A CRITERION FOR THE STABILITY OF THE LITHOSPHERE. 


Because the lithosphere is everywhere supported on a 
substratum having negligible strength, it is now possible to 


2 Ross GuNN. [J. F. I. 
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establish the criterion for vertical equilibrium in any selected 
prism. In the asthenosphere below the lowest section of the 
lithosphere, we adopt some surface of equal pressure or equal 
gravitational potential as a reference. If an applied vertical 
load on the lithosphere should increase the pressure at this 
selected level, transient horizontal pressure gradients will be 
set up that force the weak rock outward more or less hori- 
zontally until the constancy of pressure over the selected 
equipotential surface is reéstablished. Calculations suggest 
that such transient states last only about 25,000 years and 
therefore are not of much geological interest. It is a property 
of the hydrostatic state that when constancy of pressure is 
reéstablished, the system is in equilibrium. Now the pressure 
at this selected level is made up of three components, namely: 
(a) the weight per unit area of the column of weak rock 
extending from the reference level up to the lithosphere; 
(6) the weight per unit area of the column that extends from 
the bottom of lithosphere to the surface; and (c) the vertical 
stress imposed on the column by the finite strength of the 
lithosphere or by some superposed load. The total pressure 
arising from (a) and (6) may be conveniently written M9 
where M is the average mass per unit area lying above the 
reference equipotential surface and g is the average gravita- 
tional acceleration in any selected prism. 
The criterion for vertical stability therefore is 


P = Mg+ S = Constant, (4) 


where P is the pressure at the selected equipotential or 
isobaric level and S is the superposed vertical stress. The 
relation is generally applicable even though each prism may 
be rigidly attached to its neighbors. Sis taken to be positive 
when it acts in the direction of g or downward. This exact 
result, discussed by the author in an earlier paper,? may be 
stated in words as follows: For sections of the lithosphere to 
be in stable or isobaric equilibrium it is necessary and suff- 
cient; in prisms reaching from a selected surface of equal 


*R. Gunn, “Quantitative Study of Mountain Building,’”’ Jour. FRANKLIN 
INST., 224, 19, 1937. 
T. A. Elkins, Geophysics, 7, 45, 1942. 
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pressure in the asthenosphere to the actual surface, that the 
the average externally applied vertical stress added to the prod- 
uct of the mean mass per unit area and the average acceleration 
due to gravity, shall be a constant. 

When the lithosphere is adjusted to this condition, one may 
conveniently say it is in isobaric equilibrium, since the word 
isobaric implies equal pressure. Isobary * is, therefore, the 
state of vertical mechanical equilibrium even though it per- 
mits real anomalies of gravity to exist. Obviously isobary re- 
duces to isostasy as a special case in level regions where ver- 
tical stresses are negligible. It is convenient practically to 
pick some standard column asa reference. A suitable column 
might be like Heiskanen’s ‘ isostatic standard column except 
in most cases it would include a portion of the asthenosphere. 
It is important, however, to select the standard column near 
the center of a large undeformed region. At such a place one 
may be assured that the superposed vertical stresses usually 
due to deformation are small and the regional isostatic 
adjustment is nearly ideal, irrespective of the finite strength 
of the lithosphere. Now if gravitational measurements at 
other places are compared with those of such a reference 
point, it is clear that gravity anomalies will be an indication 
of vertical stress. Eq. (3) may be combined with Eq. (4) 
and the principle of isobaric equilibrium written in the 
approximate but useful form 


1 
ee 


P = Mg- Bag _ Constant. ( 
27Y 

This relation may be made exact if the geometrical distribution 
of stress, mass and topography is available in any given case. 
The principle of isobaric equilibrium as here expressed, to- 
gether with observed numerical values for the gravity anom- 
alies, provide a valuable tool for the study of lithospheric equi- 
librium, particularly special geological structures that are 
definite exceptions to the principle of isostasy as usually 
defined. 


* C. E. Dutton considered the use of the term isobary to describe isostasy in 
1889 but rejected it for reasons that now appear trivial. 


‘W. Heiskanen, Gerlands Beitrige fiir Geophysik, 36, 201, 1932. 
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DEFORMATIONS OF A HYDROSTATICALLY SUPPORTED ELASTIC LITHOSPHERE. 


The lithosphere of the earth is the outer shell that is 
capable of maintaining relatively large shearing stresses for 
long periods of time. This outer layer is supported every- 
where by the underlying asthenosphere that is assumed to 
have no such strength. The lithosphere, below about 10 Km. 
is subject to a vertical compressive stress that exceeds 3.10° 
dynes/cm.?. Such a stress can initiate plastic flow in any 
ordinary rock so that voids due to fracture probably cannot 
exist. Further Bridgeman’s experiments® strongly suggest 
that fractures occurring under such circumstances will be self- 
healing. For this reason it seems logical to suppose that, 
except for the upper 5 or 10 Km. that may be traversed by 
many fractures, the lithosphere is usually a homogeneous 
sheet covering the earth, and one may discover its properties 
by considering the deformations produced by known loads 
superposed thereon. The mathematical complications im- 
posed by three dimensional solutions make such investigations 
most laborious and impractical. However, several important 
two dimensional problems may be solved and these, because 
of the essentially linear nature of many earth features, are 
capable of giving an adequate description of the properties of 
the lithosphere. In the following discussions we consider the 
deformations of the lithosphere in sections taken at right 
angles to long linear blocks of sediments, long narrow oceanic 
deeps, and continental boundaries. On comparing these 
idealized sections with earth geology, a good quantitative 
idea of the fundamental physics involved in the production 
and maintenance of such surface features may be derived. 
It will be shown in this and succeeding papers that a litho- 
sphere having uniform elastic properties and finite strength 
will explain in some detail most of the observed anomalous 
departures from isostasy. Moreover, a proper application of 
the principles of isobary as outlined above shows that many 
regions heretofore classed as questionable are, in fact, in 
definite mechanical equilibrium even though large gravity 
anomalies exist. 

Exclude the upper 5 kilometers of the earth’s crust and 
consider a section of a uniform lithosphere resting on magma 


" 5 P, W. Bridgeman, Jour. Geology, 44, 666, 1936. 
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that supports it hydrostatically. For convenience assume 
that the section is of unit width, effective thickness 7, its 
uniform elastic modulus is £, and its section moment of 
inertia is J. Reference to any good text on strength of 
materials will show that as such a section or beam is deformed, 
the radius of curvature r is given by 


+ EI’ (9) 
where M, is the magnitude of the bending moment at the 
point where 7 is measured. The product EJ is called the 
flexural rigidity and is a fundamental constant of the litho- 
sphere. Now if we ignore the square of small quantities in 
the factors which determine 7 and neglect the slight curvature 
of the earth, both of which are permissible, it follows that 
d’y vv 7 
aa (7) 
where y refers to the vertical deflection of the neutral axis of 
the section, and x is the codrdinate measured along its length. 
As a general thing, M, is an unknown factor in the case of the 
lithosphere so that Eq. (7) is differentiated twice, yielding 
the useful form 


El 


oP i dV 
EI — = —-— = - 0,, (8) 

dx* dx 
where V is the shear per unit width and ¢, is the resultant 
restoring stress applied to the section as a result of its defor- 
mation. If the neutral axis of the original undeformed section 
coincides with the x axis and if vertical displacements are pro- 
duced by superposed loads or even by horizontal compression,’ 
these vertical displacements will force out the underlying mag- 

ma and produce a counter-stress 


0; = dogy, (9) 


where dp is the density of the underlying magma and y is the 
downward displacement. It is clear that the deformation at 
the surface is substantially the same as that calculated for 
the neutral axis. Further, if the y displacement results from 
the loading of the section by sediments of density d washing 
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into the region, the resultant restoring stress is evidently 


a. = (dy — d)gy, (10 
whence combining Eqs. (8) and (10) 
~@ 
EI‘ ies a = (dy — d)gy (11) 
dx' : 


or, since J the sectional moment of inertia per unit width is 
7* 12, one may write the basic equation that describes the 
deformation of the lithosphere as 


ET*d'y 


gee oe + (dy — d)gy = 0. (12) 


This expression may be compared with reference 3, Eq. (10), 
which contains another term involving the tangential com- 
pressional stress. 

The solution of Eq. (12) is of the form 


y = € "(A cos bx + B sin bx) 
+ €*(C cos bx + Dsin bx). (13) 


It is reasonable to suppose that the deformations of the crust 
are zero at great distances from any system of forces near the 
origin, so that one may set C = D=0. The values of A 
and B depend on the boundary conditions of the particular 
problem under investigation. 

Consider first the simplest case of a continuous section in 
which a downward force F is applied at the origin, and because 
the section is horizontal there, it is easily shown that A = B. 
The vertical deflection, therefore, is 


3Fe= ; - ; ; 
Y = Siepapa (Cos bx + sin by), (14) 
where 
ps = 3640 — DE (15) 


ET? 


Another useful solution of Eq. (13) corresponds to the case 
in which the lithosphere is vertically fractured. Select the 
origin to coincide with the fracture and suppose that a 
vertical force F and simultaneously a bending moment JJ» 


IO 
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are there applied. The force F may be due to the overriding 
of an adjacent section, or it may be due to some externally 
applied force. The solution for this case is readily shown to be 


6 e br 
y= 


4 OE [ F cos bx — bMo cos bx + bMy sin bx]. (16) 


Upon comparison with Eq. (14) it may be seen that the 
maximum deflection of a broken beam for a given vertical 
load is four times that of a continuous one. Thus, the most 
unusual geological characteristics of the lithosphere will ap- 
pear in such regions of failure. The deflection of the litho- 
sphere in all these cases is evidently wavelike, in character, 
the wave-length \ being 27/5. For convenience and reference, 
Eqs. (14) and (16) have been plotted in Fig. 1, employing 
identical constants and taking My) = 0. 
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Because most loads on the earth are distributed over a 
considerable area, it is important to calculate with greater 
accuracy the deformation of the lithosphere due to any 
arbitrary linear distribution of load. Suppose that the 
section is continuous, uniform, of unit width, and that a 
load W per unit length is applied. The deflection dy due to 
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an element of this load Wdx is, by Eq. (14), obviously 


ty = SUS (cos bx + cin bx) - 
dy = RT (cos bx + sin bx)dx, (17) 
where x is the distance from the point at which the elementary 
load Wdx is applied to the point where the beam is deflected 
a distance dy. Eq. (17) has been integrated for certain 
important load distributions but we here consider only 
uniform loading. Select the origin as the point at which the 
total deflection y is to be evaluated. Then if the superposed 
uniform load per unit area Aodg = W extends from x = — L 
to x = N, Eq. (17) may be evaluated. The total deflection 
at the origin is found to be 


hod 

ve ce 

4d, 

where d is the density of the superposed block, do that of the 
substratum, and hy is the height of the block. In this relation 
y may be expressed more conveniently as a function of a new 
variable x measured from the center of the block. Evidently 
then N = m+ xgand L = m — x, where m = (N + L)/2, 
so that 
hod 
ita do 


4 — 26°" cos DN — 2€7°* cos dL], (18) 


fr — de" cos b(m + x) 


— 56°» cos b(m — x) }. (19) 


Now if m+ x> 1/b, that is, if the block has reasonably 
large dimensions compared to 1/6, then Eq. (19) may be 
further reduced to a form most useful near the boundary of 
the block or 
__ hed 
7 dy 


where g is measured inward from the end of the block. 
Because the lithosphere just under the end of the block is 
a point of inflection and M,; = 0 there, it is easy to show that 
for distances measured away from the edge of the block 


[1 — 3€°* cos bq ], (20) 


hod _, 
__ —Op - 
4 i cos bp. (21) 
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It should be clear that in all of these relations the ampli- 
tude of any deformation decreases from its maximum value 
to zero in distances varying from x = 2/2) to x = 2/b so that 
b is evidently a most important lithospheric constant. Its 
value and significance will be discussed in a later section. 


FRACTION OF SUPERPOSED LOAD CARRIED BY LITHOSPHERE AS A 
FUNCTION OF LOAD SIZE. 


One of the important problems of isostasy is the determi- 
nation of the size of the smallest superposed load that will be 
substantially in isostatic adjustment. This size is especially 
important in evaluating regional isostatic anomalies. The 
critical size may readily be calculated from the above relations 
and the principle of isobaric equilibrium. Let Z be the 
distance from the top of the originally undeformed lithosphere 
to the selected uniform pressure level in the asthenosphere 
and JT be its thickness. Then Z — T is the distance from 
the bottom of the lithosphere to the reference level. Further, 
let d be the density of a superposed block of sediments of 
height Ao, D be the density of the lithosphere, and dy be the 
density of the underlying magma. Without loss of generality, 
one may suppose that the acceleration due to gravity g is 
constant. At the assumed reference level the pressure outside 
the deformed region is 


bo = [Z — T ]dog + TD (22) 
while the pressure at the same reference level under the 
deformation is 

bu = hdg + TDg + (Z — T — y)dog + S». (23) 


But by the principle of isobaric equilibrium and Eq. (4) these 
pressures are identical and may be equated to each other. 
Thus it is found that 


Sy + hodg — ydog = 0. (24) 


This may be generalized and the issobaric stress function 


written as 
S, = ydog — Q, (25) 


where Q is the applied load per unit area at a selected point, 
y is the downward deformation of the lithosphere and S is 
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the vertical stress provided by the finite strength of the 
lithosphere at the same point. As before, a downward stress 
is taken to be positive. It will be noted that Eq. (25) is 
independent of the density of the lithosphere, and reduces to 
a statement describing isostasy if. S = 0. This useful relation 
will be called the stress function and is invaluable in all 
discussions involving gravitational anomalies. 
Employing Eqs. (19) and (25) one may write 


— S, = hodg — ydog = (Visos — V)dog 


hidg 
e Ss [e b(m+rz) cos b(m + x) + € b(m—r) cos b(m aA x) | (26) 
~ 
where Yisos = Nod/do. 


Now because erosion and sedimentation will round off 
any practical block of the assumed type and transfer material 
to the adjacent depressions, Eq. (26) can probably not be 
accurately applied to regions near its boundary. It is 
possible, however, through this relation to express in regions 
near the center of the block the ratio of the expected gravita- 
tional anomaly to that which would be observed were the 
load carried by a perfectly rigid lithosphere. This same ratio 
also is the fraction of the load carried by the lithosphere and 
will therefore be zero for a region in perfect isostatic adjust- 
ment. Employing Eq. (3) and taking x = 0 because regions 
near the edges may not be significant, it is found here that 


iM Ag 
( - = { —— = €°™ cos bm. (27) 
hodg r=0 2ryhod z=0 


This ratio is plotted in Fig. 2 as a function of bm. It is 
interesting to note that small negative anomalies are possible 
even at the center of a large positive load of sediments, 
provided bm lies between 7/2 and $r. This detail, however, 
may well be illusitory in geological structures if the erosion of 
the block edges is considerable. In any case, it is evident 
from Eq. (27) and Fig. 2 that a lithosphere having a finite 
strength in shear does not permit a superposed block load to 
be even in approximate (96 per cent.) isostatic adjustment 
unless its angular semidimension bm exceeds 1.4. The per- 
fection of isostatic adjustment is measured, of course, by the 
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fraction of the superposed load carried by the immediately un- 
derlying magma, and this in turn is unity minus the fraction 
of the load carried by the lithosphere. For convenience, the 
average fractional load carried by the lithosphere for blocks of 
various semidimensions )m is also plotted in Fig. 2. 
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Now in order to calculate m the semi-linear dimension of 
the superposed load for any selected degree of approach to 
isostasy, it is necessary to know the value of b. Approximate 
values of 6 may be calculated immediately from Eq. (15) by 
assuming a value for the effective thickness of the lithosphere 
T. If in Eq. (15) one takes E = 10” dynes/cm.?; dy — d 
= 3 gms./cm.’ (i.e., assume d = 0), T = 10’ cm. and g = 10° 
cm./sec.”, it is found that b = 5.4-107* cm.~', or if T = 5.10° 
cm., 6 is 8.4-10-§ cm.~!.. Adopting this latter value as the 
most probable because of the values deduced in a later 
paragraph, it follows that a block of at least a semidimension 
m = 167 Km. will be necessary before even its center will 
approach 96 per cent. compensation or within 4 per cent. of 
perfect isostatic adjustment. A block of semidimension 
m = 25 Km., or total width of 50 Km., a value sometimes 
given as the extent of the block necessary to achieve isostatic 
adjustment, will only be 19 per cent. compensated, and 
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81 per cent. of the load will be carried by the lithosphere. 
This lithospheric load is passed on to the asthenosphere in 
regions outside the block for a distance approximating 1/b cm. 
It seems clear, therefore, that unless the values of the im- 
portant lithospheric constants employed above are seriously 
in error, the areas necessary for the establishment of isostatic 
equilibrium must be at least three hundred kilometers square 
and the methods currently employed to estimate gravity anom- 
alies are not valid. In a later section, the important constant 
b is calculated by two other independent methods not related 
to isostasy and a value approximating that adopted above 
is found. 


GRAVITY ANOMALY AND DEFORMATION NEAR A THRUST FAULT. 


We next investigate the distribution of the gravity 
anomalies that may be expected to result from a thrust fault 
if isobaric equilibrium is maintained. Assume that the 
horizontal compressive stress S; in the lithosphere has been 
great enough to produce a thrust fault and that the plane of 
the fault makes an angle ¢ with the horizontal. Then the 
downward force at the origin due to the overriding section is 
S,T cos ¢ sin ¢ where T is the thickness of the section. The 
stress function acting on the depressed section from Eqs. (16) 
and (25) is evidently 
old, hey ee ee ee 

: EBT? 
whence, from Eq. (3), the gravity anomaly in the depressed 
region approximates 
— 121r(dy — d)yS;, cos ¢ sin ge” cos bx 
: ct leat osama ce cecum IE (29) 
Eb®T? 


Gb 
re 
II 


A nearly similar expression may be obtained for the positive 
anomalies in the overriding region by making allowance for 
the erosion after the original equilibrium is attained. Upon 
substitution of b = 7-107§ cm.~!, dy — d = 2 gm./cm.* (be- 
cause the depression is at first filled with water), T = 5-10° 
cm. and S = 10° dynes/cm.’ it is found that 


Ag = — 0.280e~* cos bx. (30) 
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The actual tangential compressional stress acting on the 
thrust fault cannot be directly evaluated but because thrust 
fracture and the formation of arcuate mountains cannot take 
place without sufficient stress to induce considerable plastic 
deformation, the assumed stress seems reasonable. A value 
of 2 X 10° dynes/cm.? may be equally probable so that one 
may only say that Eq. (30) is consistent with the known 
observational facts and that the maximum anomalies approxi- 
mate from 200 to 500 milligals. Further, although Eq. (30) 
implies a sharp discontinuity in gravity anomaly coincident 
with the fracture, it is clear that erosion and sedimentation 
will always act to erase such rapid changes. 


THE VALUE OF THE LITHOSPHERIC CONSTANT b. 


Because the value of the lithospheric constant 6 is of first 
rate importance in all discussions of the stress and deformation 
of the lithosphere, it is of interest to examine various known 
geological formations to see how } changes under widely 
varying conditions. In general, one selects simple structures 
in which the acting forces are both known in position and 
localized in nature. For example, the downward force acting 
on a thrust fault is applied at the fault. Gravity anomalies 
near such a fault are given by Eq. (29) and in many cases 
may be observed. From Eq. (29) the width of the principal 
negative anomaly strip and the adjacent principal positive 
strip are both given by bx = 7/2 whereas the secondary 
positive strip outside the principal negative region has a 
width bx = x. Upon comparison of these values with actual 
faults, 6 may readily be determined. There is considerable 
geological evidence showing that arcuate island or mountain 
chains result from great systems of tangential stresses in the 
lithosphere and that this opportunity for overthrusting 
frequently produces long linear geological features. One of 
the best examples is the East Indian chain, particularly 
south of Java, where the geometry of the deformation is 
remarkably simple and is suitable for comparison with the 
present calculations. It is assumed that this anomaly strip 
is a direct result of a long thrust fault. 

The gravitational characteristics of the sea south of Java. 
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have been well studied by V. Meinesz ® and many facts are 
available that have been summarized in a well known map 
by van Bemmelen. Scaling from this map it is found that 
the negative anomaly strip varies in width from 150 to 185 
Km. Adopting an average value of 170 Km. and setting this 
equal to 2/20 it is found that 6 = 9.2 X 107§cm.~!._In the 
same manner the principal and secondary positive strips 
yield values for 6 of 9 X 107-8 cm.-! and 7.8 X 1078 cm.~! 
respectively. 

As another case of importance, consider the geological 
implications of Eq. (14) that represents the downward 
displacement of a continuous section of the lithosphere due 
to a localized load. The sediments initially deposited at the 
single mouth of a great river might represent such a great 
concentrated load, and instead of these sinking at the point 
of deposition establishing isostatic equilibrium, Eqs. (14) and 
(27) show that the first sediments deposited will be almost 
entirely supported by the lithosphere. Therefore, it will be 
bent down in a wide basin of half width bx = 37/4. This 
wide basin of downwardly deformed ocean bottom encourages 
the deposition of sediments at a considerable distance from 
the original mouth, and a great symmetrical delta will be 
built up whose semidimension at right angles to the axis of 
its development is 37/4). In general, sedimentation will 
proceed laterally until the downward deflection is filled, at 
least thinly, to the edge of the original depression. When 
such a delta is well developed, as in the case of the Nile and 
Niger Rivers,’ the fact that no large gravity anomalies are 
observed guarantees that the piled up sediments on the ocean 
floor reflect the figure of the downward deformation of the 
lithosphere. Barrell’s data indicate the semidimension 37/4 
for both of the above-mentioned deltas to be about 300 Km. 
Thus the characteristic constant for the lithosphere 0 is found 
to be 6 = 7.85 X 107° cm.~! or the characteristic distance 
1/b = 127 Km. The average of the four estimates above is 
8.4-10°° cm.~1. Perhaps real variations of thickness or 
modulus of elasticity exist, and one is therefore justified only 


®6V. Meinesz, ‘‘Gravity Expeditions at Sea,’’ Vol. 2. 
7]. Barrell, Jour. Geology, 22, 1914. 
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in adopting 
b = 8.4 + 2.0 X 107° cm... 


Thus the characteristic relaxation distance of the lithosphere 
is 1/6 or 120 + 25 Km. 

Further by aid of Eq. (15), one may estimate the effective 
thickness of the lithosphere by assuming E = 10” dynes/cm?’. 
and dy —d = 2 gms./cm.?. From this it is found that 
T = 50 + 15 Km.,avalue that happens to be in fair agreement 
with that inferred from isostasy using the Airey hypothesis. 


CONCLUDING REMARKS. 


The foregoing analysis appears to show that the presumed 
approach of the earth’s crust to exact isostatic adjustment has 
been somewhat exaggerated in the past and that some of the 
corrections applied to geodetic data have been carried to an 
accuracy beyond those values actually justified. Large sec- 
tions of the earth do tend toward a precise regional isostatic 
equilibrium, but notable departures are necessary in special 
phases of geological evolution. The formation of great 
geosynclines, mountains, ocean deeps, and overthrusts at 
some time in their geological history were all characterized by 
departures from isostasy even though isobary was approxi- 
mately preserved. Generally speaking, isostasy describes the 
equilibrium of the lithosphere to a sufficient approximation in 
level or undeformed regions but near mountains, continental 
boundaries, oceanic deeps, or fracture zones, the more general 
principle of isobary must be employed to account accurately 
for all the observed phenomena. 

In a future paper the methods outlined above will be 
applied to the study of the Hawaiian mountain range and to 
continental boundaries. 
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Reclaiming Worn Files.—(/Jron Age, Vol. 151, No. 6.) The 
useful life of an average file can be extended approximately three 
times its normal life by a reclaiming process used by the R.C.A. 
Mfg. Co.’s plant at Indianapolis. This process was suggested by 
Jason R. Harris, a plant employee, who received a W.P.B. certificate 
of individual production merit for the suggestion. Mr. Harris’ 
suggestion was thoroughly investigated by the plant’s engineering 
department and it was found to be very effective and applicable 
to files used on both non-ferrous and ferrous metal. At present 
all files used at the Indianapolis plant are being reclaimed by this 
method. The reclaiming process is as follows: (1) Clean the file 
with a wire brush wheel to remove deeply imbedded particles. 
(2) Clean and rinse file thoroughly, using ordinary degreasing 
solution. (3) Dip in solution made up of 800 c.c. water, 150 c.c. 
nitric acid, 100 c.c. sulfuric acid. (4) Cold water rinse. (5) Hot 
water rinse. (6) Cyanide dip to neutralize. (7) Cold water rinse. 
(8) Hot water rinse. (9) Dry thoroughly. The time of immersion 
in the solution (item number three) is dependent upon the condition 
of the files. Experiments will indicate the best time periods. 
The minimum dipping time is approximately 1 min. To save 
time, presorting of the files is recommended. Groups of files 
showing approximately the same wear should be dipped together. 
The files should be kept in an upright position while immersed so 
that the solution can act on both sides. R.C.A. engineers report 
that some of the rinse steps listed may be eliminated and fairly satis- 
factory results obtained. However, the steps listed here were found 
to give the best results. Using this procedure, the total handling 
time for 25 files is about 30 min. 


R. H. O. 


THE HORIZONTAL PATH ALONG THE EARTH’S SURFACE 
OF A PROJECTILE OR PLANE UNDER CON- 
STANT TANGENTIAL ACCELERATION. 


BY 
W. S. KIMBALL, Ph.D., 


Lt. Comdr. U.S.N.R., U. S. Naval Academy, Annapolis, Maryland. 
OUTLINE. 


The additional term in Coriolis’ normal force due to the 
earth’s rotation previously derived ! by vector methods, and 
resulting from the increased rate of direction change that 
arises from tangential path acceleration along the surface of 
the rotating earth, is here derived again by the wedge method 
of showing normal path displacements. It is further shown 
that every tangential path acceleration along the surface of 
the earth rotating at constant angular velocity w, = w sin @ 

P P as da ZW 
at latitude ¢ involves additional path curvature, —- = ~—, 
eo 29 
where v is the velocity, and due to the earth’s rotation; and 
when this curvature is multiplied by the velocity, v, of the 
acceleration field, we obtain the extra rate of direction change, 
9 
me me due to the earth’s rotation where tangential 
acceleration occurs. Thus, Coriolis’ force using the entire 
rate of direction change due to the earth’s rotation becomes 
20'v = 2UV (w. -+- nA nt 
oe 
rived, and the second term is now seen to arise from the extra 
path curvature that is due to and occurs with tangential 
acceleration along the surface of the constantly rotating earth, 
as well as to the effect itself of the earth’s rotation. 

A special study is made of the path of a body all of whose 
energy of motion is due to constant tangential acceleration, 
such as an ideal airplane of constantly increasing speed. 
The path is found to be the involute of a circle, and the rate 
of change of direction is constant and about a degree per 


) = 2uw, + 3wiv, as previously de- 
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minute in middle latitudes. If 100 miles per hour per hour 
is taken as the constant tangential acceleration, the plane will 
completely reverse its direction in about three hours at a per- 
pendicular distance of about 285 miles from the starting point 
of rest, and at a final speed of about 300 miles per hour. 
These results are shown to be in accord with the speed and 
distance requirements of elementary mechanics using only 
the familiar Coriolis’ force term C = 2w ,v for constant gradi- 
ent velocity across a uniform parallel force field on the earth 
as previously ! treated. 
NO. 1. THE NORMAL FORCE DUE TO TANGENTIAL ACCELERATION FOR 
HORIZONTAL MOTION ON THE ROTATING EARTH. 

We, here, treat the case of horizontal motion on the earth 
rotating at constant angular velocity w, = w sin ¢ at latitude 
@. The acceleration with respect to fixed codrdinates is given 
by the well-known vector experssion: 

d’s_ dv ; 

— = — + 2w, Xvu+orz X (wr X 1), (1) 

er 
where v is the velocity relative to the moving earth, and 7 is 
the radial distance to an origin on the earth’s surface about 
which a tangent plane rotates, and the Gibbs’ notation for 
vector products is used. The extra term due to rate of change 
of w, is absent since w, = constant. Since we are dealing 
with horizontal motion v and w,; are perpendicular and so (1) 


reduces to: 
d’s = dv 4 = 
cmeeee «ieee 20:19 — 217 
ar 
Since we here consider motion along the earth’s surface rather 
than out on the tangent plane away from it, we set r = 0 and 
consider the acceleration at the origin where the tangent plane 
touches the earth: 
d’s dv " (2) 
~ = —U, + 2w Vue 2 
dt? dt 


where the unit vectors “4, and u» are taken along and normal, 
respectively, to the path, and the second term will be recog- 
nized as Coriolis’ normal acceleration for the constant rate 
of direction change w,. 
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The present paper is concerned with departures from (2) 
arising from the fact that w, will not give the rate of direction 
change due to the earth's rotation where tangential accelera- 
tion occurs on a curved path. For example, the rate of di- 
rection change for constant speed on a circular path is not 
the same as when tangential acceleration occurs on the same 
path. 

It was! shown by vector methods, such as above, that 
tangential acceleration along a horizontal path entails an ex- 


° - . . ‘ U; 
tra rate of change of direction Aw, = 5 OL due to the earth’s 
v 


rotation and expressed in terms of the ratio of the velocity z, 
of the acceleration field to the complete velocity v.  Pre- 
viously, v, was ¥ v? — 2/1 where // is the energy of the har- 
monic field which is the straight line projection of energy 
stored at constant speed on the potential energy circle, where 
it is balanced by the effect of the earth’s rotation, and hence 
contributes nothing to the tangential acceleration field energy 
with its velocity v,. We treat now the more general case 
where the acceleration field energy, necessarily a part of the 
total energy of motion (relative to the earth’s surface) is 
arbitrary and no longer restricted as above. This gives rise 
to a corresponding term of the Coriolis’ force K = 3w:?, 
which, where tangential accelerations exist, is essential to 
the complete treatment of the problem: 


20'v = 2w,v + 30,7, = C+ K, 


, 3 vy a 
AO =a,t+ ane § (3) 
-~ c 


where 6’ is the complete time rate of change of direction due 
to the earth’s rotation and 26’v is the complete Coriolis force, 
given by replacing wz of (2) by @’ of (3). 

For the case of unit mass going from rest to the steady 
state of constant gradient velocity across a uniform field, it 
was found that the work of path deflection was equally di- 
vided between the two normal forces C and K, showing that 
K = 3w 2, was of the same order of magnitude as the familiar 
C= 2W VU. 

Its importance would thus seem to justify a second deri- 
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vation by the wedge method of showing normal displacements 
due to constant angular velocity w, in case of a given path 
velocity and acceleration according to Fig. 1. These dis- 
placements are shown for constant velocity v, and constant 
acceleration ‘‘a,’’ where it is understood, as usual in calculus, 
that these are instantaneous velocity and acceleration at a 
given place and time, the velocity being the distance that 
would be covered in unit time if the rate of change of distance 
were kept constant over unit time interval, and similarly with 


7 3 
Le S=29stawt 


S=A9« Vay ? 


mir 


Fic. 1. The normal displacements due to the earth’s rotation for 
constant velocity and constant acceleration along the path. 


the acceleration. The normal displacement s = 7@ is ex- 
pressed in terms of the time ¢, using the constant angular 
velocity relation 6 = w,t, and radial distances along the path 
are given by r = vt andr = af? asin Fig. 1. Thus: 


Il 


s = r0 = ww_f* for velocity, 


II 


s = 70 = }aw,f* for acceleration, 
and the second time derivative of these is: 
d*s 
dt? 


= 200, = C, 


d?s 


dt? 


= 3atw 1 = 3001, = K, 


FS Ee Re OO 


See Se 
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which are seen to be the same as expression (3) previously 
derived for the Coriolis’ normal acceleration terms due to 
velocity and acceleration, and we note how v, comes into the 
formula for K as the velocity due to tangential path accelera- 
tion. In the previous vector treatment the additive nature 
of C and K was demonstrated * by showing each to be an 
additive constant of integration obtained by the partial in- 
tegration of (11)*, yielding K as the constant of integration 
for (10)* and C as the constant of integration for (13)*. In 
the present treatment, their additive nature is justified by 
the elementary fact that displacements due to acceleration 
may be separated from and superposed upon those due to 
velocity, to give the resultant effect. Thus for constant 
initial velocity vo, and constant negative acceleration —a, we 
have r = vot — at, and the corresponding normal displace- 
ments, shown in Fig. 1, and their time derivatives above, will 
clearly be of opposite sign and additive (or subtractive) rates 
of change, and, in general, the resultant normal acceleration 
becomes: 


dy = 2001, + 30,1, = 200, ( + 3 =) =C+K, (3) 


which is the same as (3) above. Here v is the actual velocity 
of the particle, and v, is the part of the velocity due exclusively 
to the tangential acceleration field. 

One or two special cases of interest may be noted. Sup- 
pose a particle has a constant initial velocity vp and negative 
tangential acceleration —a, then, we have for the time 
t = 2v0/5a; 

C = 2w1V = 2wz(Vo — at) = 6/5w Vo, 
K = 3010, = — 3wrat = — 6/5w 1. 


This particular time is chosen because it shows a resultant 
zero value for the normal acceleration (3) since C and K are 
of opposite sign but equal magnitude. 

Another simple case of interest is where a particle starts 
from rest and the actual velocity is entirely due to the ac- 
celeration along the field so that v = v,; then (3) becomes: 


Gn = 2010 + 301d = 5wr.v. (4) 
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NO. 2. THE PATH CURVATURE DUE TO TANGENTIAL ACCELERATION 
ON THE ROTATING EARTH. 


ra : 3 Ur f A 
Che correction term Aw; = . OL to the rate of direction 
change for horizontal motion involving tangential accelera- 
tion on the rotating earth may be shown to be a direct con- 
sequence of the path curvature which always results where 
such motion exists in a tangent plane rotating at constant an- 
gular velocity w,;. Thus, take constant radial acceleration 
‘‘a.”’ Then, the distance traversed and the angular dis- 
placement will be: 
r = jal’, 7 
= @® ut. \S) 
Elimination of ¢ between these gives the path equation in 
polar coérdinates: 


ae? 
r= -, 
20 1° 
(6) 
dr _ as 
dé w,?- 


Divide the first of (6) by the second to obtain the tangent 
of the angle between curve and radius vector: 


ne rdé 6 (7) 
any =- =—, 
dr 2 ‘ 
Then the angle of inclination of the path will be: 
i] 
a=60+y = 6+ arctan =f 
Hence the curvature is: 
1d 
' f] 2 ds i] I 
da ‘ d ae ds ae d oe (8) 
ds ds #? ds rw” 
Lj ani 
4 2a 


At the origin where the horizontal tangent plane shown in 
Fig. 2 touches the earth, r = 0 and (8) reduces to 


Pe aul (9) 


ds 2 ds 2v 


Re I Re eat 
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where on the right we have introduced d@ = w,dt and ds = vdt 
whereby the curvature is restricted to the case in hand of 


Y 


t? 


x 


Fic. 2. The curved path where acceleration is along radii that are 
turning with the earth’s constant rotation. 


motion with velocity v in a frame rotating with angular 
velocity w,, being measurable physical quantities. 

When the curvature (9) is multiplied by v, = at, the ve- 
locity of the acceleration field, we have the corresponding 
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rate of direction change: 


Aw, = v’—- = Of (10) 


If the XY plane of Fig. 2 were fixed, any particle of 
velocity v, pursuing the path there shown will experience a 
rate of direction change given by (10). If in addition the 
XY plane rotates with angular velocity w,, its total rate of 
direction change will be given by the sum of these two rates: 


=o. t+ 301. (11) 


The above results (9) to (11) are not restricted to the case 
of constant acceleration ‘‘a,’’ but hold for any acceleration 
where the distance is given by a series instead of (5): 


r=ta+be+c*..., (12) 


with first and zero power terms absent since we treat only 
acceleration effect with zero initial velocity and displacement. 
The path is given as before by removing ¢ by means of 0 = wzt. 
The series expression for the path when divided by its de- 
rivative gives a series expression for tan y, just as (7) is ob- 
tained from (6), and when this is used in (8), the curvature 
becomes: 

da dé, dé ( I bé ) 

ds ds + S32 a, 
for terms below second and higher powers of @. At the origin 
where ry = 0 and @ = Oo, this curvature becomes the same as 
(9), which is thus seen to apply to any acceleration law given 
by (12), provided only it occurs on the constantly rotating 
earth where @ = wzt. 

We shall now apply the new expression (4) where constant 
horizontal tangential acceleration is completely responsible 
for all the motion with v, = v, to obtain the corresponding 
path equation along the surface of the earth. 
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NO. 3. THE ROUTE OF AN AIRPLANE TRAVELLING AT CONSTANT 
HORIZONTAL ACCELERATION ALONG ITS PATH IS THE 
INVOLUTE OF A CIRCLE OF RADIUS ¢ = «/25w7?. 

It is well known that unit mass actuated by constant 
velocity on the earth’s surface will follow the path of the 
inertia circle, in absence of outside forces, such as pressure 
gradient, etc., with the Coriolis force balancing the centrifugal 
force: 

C = 20,0 = v/r;7r = 0/201, (13) 


where r, the radius of the inertia circle, is shown in terms of 
the velocity and the earth’s constant rotation w; = w sin ¢ 
at latitude ¢. 

In this section we derive the path which will be followed 
by unit mass actuated by constant tangential acceleration 
along the earth’s surface and in complete absence of all other 
outside forces, so that as with (13) the Coriolis force (4) due 
to the earth’s rotation will balance the centrifugal force. No 
airplane or projectile is, of course, ever actuated by such 
exclusive acceleration, but if it were, then this is the path it 
would follow. We take zero initial velocity, so that the 
velocity of the tangential acceleration field is the same as 2, 
the complete actual velocity, and all of the energy of the ac- 
celeration field goes into energy of motion along the path, 
none of it being lost in the form of energy stored on a poten- 
tial energy circle, since there is no outside force field to 
balance it. 

Thus instead of (13), we have, using (4): 


C+ K = So. = v/r = — vda/ds, (14) 


where on the right the curvature is introduced using a@ for 
the angle of inclination of the path and s for the distance 
along the path. Since we have constant tangential accelera- 
tion, a, the energy and velocity are given by integration 
along the path: 
8 y 
f ads = as = —;v* = 24s. (15) 


2 


Hence from (14): 


da/ds = — 5w,/V 2as, (16) 
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which integrates to: 


T 2S ‘ 
—-a=B = 50o1.V— (17) 


a 


E T i ei : 
taking a = >, ats = 0, asin Fig. 3, and @ is the total angular 


~ 


range of the path. Then: 


; , | 2s : 
sin (7/2 — a) = cosa = dx/ds = sin { 5a, V— ) =sin B, 
Wr 


a 
; 2s 
~= sin 5 — ds 
a 


= f sin \ ds = cf asin 6BdB, (18) 


y 2s 
where c = @/25w 7”, B = J - and ds = cBdB. 
c 


Upon integrating again, we have: 
x = c(sin 8 — B cos B) (19) 


Again taking the cosine of 7/2 — @ as given by (17): 


; | 2s | 
cos (7/2 — a) = sina = dy/ds = COS 5W1 \ — = cosf, (20) 
a 


and with the same substitution as in (18), we obtain by in- 
tegration: 


y +c =c(8sinB + cos £8), (21) 


where the constant of integration, c, provides that y = 0 at 
the start where 6 = 0 according to Fig. 3. Note that the 
maxima and minima of y occur at x = +c where (20) is 
zero, and the maxima and minima of x occur at y = 0 or 
—2c, and that x = + $c at y = 0, as well as at its maxima 
on the line y = — 2c. Equations (19) and (21) are well 
known as the parametric equations of the involute of the 
circle of radius c shown in Fig. 3, which is the locus of the 
centers of curvature of the involute. 

Elementary mechanics affords a check on the path shown 
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in Fig. 3, as may be seen by comparing it with the path fol- 
lowed in a uniform force field on the rotating earth by unit 
mass starting from rest, as shown in Fig. 2 of the previous ! 
article. In that case, the path extended distance 6 = a/8w,? 
along the direction of the force field before attaining gradient 


bd 


Care, 0) (8.9) 
, X 


e4 19 
(7 <,-a¢) 


Fic. 3. The involute of the circle of radius ¢ = a/25w,? which is the path 
of a plane or projectile actuated by constant tangential acceleration ‘‘a” along 
the surface of the constantly rotating earth. 


rectilinear velocity across the field after a 90° deflection. 
Compare this with the first maximum of x in Fig. 3 where the 
derivative of x given by (18) is zero, and hence 6, the angular 


ie ae a 
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path range is m and by (19) xX» = mc = ma/25w,? after a 180° 
path deflection. 

The notable thing about these two perpendicular dis- 
tances, 6 and x,,, from the starting point of rest is that they 
are approximately equal (compare (22) and (25) below) with, 
for b, a 90° deflection and for x, a 180° deflection, using the 
same acceleration, a. This shows that the deflecting force 
of the present case of Fig. 3 averages about twice as great as 
for the previous problem, since the angle is doubled for about 
the same normal distance to the starting point, a result which 
checks the mechanics of the situation. For, while the present 
problem involves no field acceleration to resist deflection, the 
previous case of a uniform parallel field, supplied a force 
squarely against the y component of the deflecting force due 
to just such a field resistance and extending over half the 
angular range in making a 90° turn. Thus for the uniform 
parallel field half the angular deflection is to be expected in 
about the same distance from the starting point. Moreover 
this distance, 6, was seen to be independent of the path and 
entirely determined by the classical Coriolis term C = 2w 1, 
so that when the deflection angle is doubled in the absence 
of a resisting field of acceleration as in the present case of 
Fig. 3, the result is seen to be consistent with and thus checked 
by this classical mechanics. 

We may compute this first maximum of x given by (19) 
when its derivative (18) is zero. It occurs where the angle of 
path deflection 8 = 7m (see Fig. 3) and we may take for middle 
latitudes ¢ = 53°; sin @ = .8, so that w,; = w sin @ becomes 
4/5 the rate of the earth’s rotation about its polar axis, i.e. 


7 wa ‘ 
w, = 12°/hr = —rad./hr. Thus using (19), we have 
15 
Ta 2254 900 : 
Xm = 1C = ——, = —— = — = 286 miles. (22) 
25 1~ 257 T 


If we assume a small * but constant acceleration of 100 miles 


* 100 miles/hr.2 = 4-4 inches/sec.2 = about } inch/sec.2 = about 1} thou- 


sandths of the acceleration of gravity. For large accelerations leading to dis- 
tances comparable with changes in latitude, w, would not be constant and a 
corresponding modification of the present treatment would be required. 
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per hour *, then (22) becomes about 286 miles as shown on 
the right. 

From the distance formula for constantly accelerated mo- 
tion and the relations between a, 8, c and s given by (18) 
we have: 

B 2 
s = fal? = hcp? = fa (£ ) ort = B/5wx. (23) 
OWL 
Hence, we note that when we use a = 100, this first maxi- 


mum, X» occurs three hours after starting from rest, and at 
a velocity of 300 miles an hour: 


t = B/5w1, = (m/5)(15/4) = 3 hours, 
v = at = aB/5w, = 300 miles per hr., 
dB/dt = 5, = 2/3 radians/hr. = 60°/hr. = 1°/min. (24) 


Thus an airplane actuated by this constant tangential 
acceleration of 100 miles per hour ? will completely reverse its 
direction (8 = 7m) in about 3 hours and at a perpendicular 
distance of almost 300 miles from its starting point of rest. 

On the other hand for the uniform parallel field of accelera- 
tion previously discussed the gradient velocity at right angles 
to the field where the Coriolis force is (13) and equals the 
field we have, using the above data: 


@ = 2w1¥;v = a/2w, = 100/2 X 15/m = 240 miles/hr. 
and the perpendicular distance from the starting point is 
b = a/8w,? = 100(225)/82? = 285 miles (25) 
for 8 = r/2 = 90° or half the angular deflection of (22). 
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Glass Fiber.—Parachute flare shades woven of glass fibers are 
being used to shield Army bombing planes and the eyes of bom- 
bardiers from the glare of flares dropped to illuminate the target. 
Field telephone cable used by the Signal Corps is now made with a 
tension core of glass yarn to replace the metal formerly used for this 
purpose. Navy deck motors employ glass fiber insulation because 
of its resistance to dampness. These and other Army and Navy 
uses of glass fiber materials were described by Major Hubert D. 
Keiser, Services of Supply Resources and Production Division, 
United States Army, in an address at the 158th meeting of the 
American Institute of Mining and Metallurgical Engineers at the 
Engineering Societies Building, New York City. The glass cloth 
parachute flare shades are used in connection with brilliant mag- 
nesium flares of about one million candle power which are employed 
to illuminate the objective of a bombing attack with almost the 
brilliance of daylight. If the light of the flare were not shielded 
overhead the bombardiers could not see the ground because of the 
intense glare, and the attacking planes would be revealed to enemy 
anti-aircraft positions. A light flare shade is therefore required, 
just above the incandescent element. Glass cloth, darkened with a 
colloidal graphite finish, is able to meet the requirements because 
it is light in weight and can withstand the intense heat of the flare. 
Describing Navy uses of glass fiber materials, Major Keiser said, 
“The most important single property required of materials used in 
Navy ships is incombustibility or fire safety. Reliability under 
extreme conditions of service is the primary requisite of equipment, 
and reduction of weight is always desired to increase fire power, 
cruising radius, or speed. For these reasons fibrous glass is ex- 
tensively used in Navy cables, motors, generators and fire controls.”’ 
Aircraft uses of glass fiber materials, according to Major Keiser, 
include glass fiber blankets for heat and sound insulation in hull 
and cabin, according to type of ship. 


R. H. O. 
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CONGRESS OF PHYSICS AT PUEBLA, MEXICO. 


The First National Congress of Physics was held under the 
auspices of the Governor of the State of Puebla, the Rector 
of the University of Puebla, and the Director of the National 
Astrophysical Observatory of Tonanzintla, at Puebla City, 
Mexico, on May 3 to 8, inclusive. 

Approximately 50 members were present, including the 
following from the United States who received special invita- 
tions from the President of Mexico to attend the Congress: 
Dr. J. G. Baker, Harvard University; Dr. J. W. Beams, Uni- 
versity of Virginia; Dr. P. W. Bridgman, Harvard University; 
Dr. L. J. Briggs, National Bureau of Standards; Dr. S. 
Chandrasekhar, University of Chicago; Watson Davis, Sci- 
ence Service; Dr. W. V. Houston, University of California; 
Dr. L. D. Leet, Harvard University; Dr. I. 1. Rabi, Columbia 
University; Dr. Harlow Shapley, Harvard University; Dr. W. 
F. G. Swann, Bartol Foundation of the Franklin Institute; 
and Dr. O. R. Wulf, United States Weather Bureau. 

The regular sessions were held at the University, with 
certain special lectures at the University and at two of the 
city’s theatres. At the opening session, the delegates were 
welcomed by Dr. C. Gonzalo Bautista, Governor of the State 
of Puebla. After the president of the Congress, Dr. Manuel 
Vallarta (formerly professor at Massachusetts Institute of 
Technology), had taken the chair, Governor Bautista pre- 
sented the Cross of Honor of the State of Puebla to Dr. 
Harlow Shapley in recognition of his achievements in galactic 
astronomy. 

In addition to five special evening conferences, seven 
regular sessions were held at which 45 papers were presented. 
Many of the visiting American scientists reviewed fields of 
activity with which their names are closely associated. They 
were later given a special audience by the President of 
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Mexico, and many other outstanding courtesies were ex- 
tended. 

A most interesting feature of the Congress was a visit to 
the new Mexican astrophysical observatory located near 
Puebla with Popocatepetl towering in the distance. Work in 
this observatory is devoted largely to the study of galaxies, 
and the observatory is equipped with a very fine Schmidt 
reflector for this purpose. Here also is located the new 
seismograph recently purchased by the Mexican Government. 
At present the seismological stations in our eastern States 
are able to get a good north-south range on seismological 
disturbances in the Caribbean, but an east-west range for this 
region has been lacking. The new seismological station will 
provide the much-needed intersecting range and thus fix the 
center of the disturbances. 

Great interest was shown by the Mexican scientists in the 
progress of research in the United States. In his address on 
the relation between the Government and scientific investi- 
gation, Dr. Briggs referred to the fact that the International 
Bureau of Weights and Measures at Paris is at present iso- 
lated from the republics of the world. In recognition of this 
condition, the Secretary of Commerce had requested him to 
convey the message that ‘‘the National Bureau of Standards 
will be most happy to carry out for the Mexican Government, 
and for other American republics, without cost, any com- 
parison of your standards with our own that you may deem 
desirable.”” This invitation was the subject of much favorable 
comment, and a special article on the work of the National 
Bureau of Standards was published in one of the leading 
newspapers. 


EFFECT OF AERATION ON SOILS. 


A number of investigators have studied the change in pH 
value of soils in passing from the moist field condition to the 
air-dried state. It is evident from these studies that drying 
of soils is occasionally accompanied by a slight change in 
reaction, generally in the direction of increased acidity. 
Since the air-dried soil offers a more uniform basis, the soil 
in this condition is generally used for determining the pH 
value. 
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Samples of a poorly-drained soil showed large differences 
in -hydrogen-ion concentration between the fresh soil and the 
soil in the air-dried condition. By varying the degree of 
aeration, the soil was passed through a cycle—from neutrality, 
as found in the natural field condition, to high acidity, 
which occurs in the air-dried state. As a result, it has been 
concluded that pH measurements of very poorly-aerated soils 
should be made only on samples maintained in the field condi- 
tion. The change in pH value observed has been associated 
with the activity of sulfate-reducing bacteria. 

The conclusions reached in this study, which will be re- 
ported in detail by Melvin Romanoff in a forthcoming issue 
of the magazine Gas, have an important bearing on the loca- 
tion of soil areas which are corrosive because of the activity 
of sulfate-reducing organisms. As these bacteria function 
only within a narrow pH range, the reaction of the soil is 
considered to be an important criterion in identifying corrosive 
areas where these bacteria are active. 


ULTRAVIOLET SOLAR RADIATION. 


Sun bathing and the resulting burns from overexposure of 
the untanned skin attest to the importance of reliable infor- 
mation on the amount of biologically effective ultraviolet 
radiation from the sun and from the whole sky, under various 
weather conditions. 

For some years W. W. Coblentz and R. Stair have investi- 
gated the instruments and methods of evaluating the bio- 
logically effective component of ultraviolet solar and sky 
radiation useful in heliotherapy; including instruments for 
measuring ultraviolet radiation incident on a surface exposed 
normal to the sun, and instruments exposed horizontally to 
receive the direct solar rays and also the rays incident from 
the whole sky. The latter instruments are useful to meteor- 
ologists in their studies of ultraviolet climatology. 

The instruments (photoelectric cells) now used in making 
ultraviolet radiation measurements are highly selective in 
their spectral response. Moreover, thus far it has been im- 
possible to reproduce photocells having exactly the same 
response, which should, preferably, cover only the biologically 
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effective ultraviolet of wave-lengths shorter than about 3200A. 
It is this radiation that produces vitamin D and is effective in 
preventing and curing rickets; also in curing lupus and other 
skin diseases 

In the June Journal of Research (RP1542), measurements 
of ultraviolet solar radiation are given, as obtained with 
photocells differing widely in spectral response. If and when 
manufacturers succeed in making photoelectric cells in which 
the response range is identical, and is preferably one that 
coincides with the band of wave-lengths that heals lupus, 
surgical tuberculosis and rickets, that produces vitamin D, 
etc., the whole procedure used in evaluating biologically 
effective ultraviolet solar and sky radiation will be greatly 
simplified. 

The data in the present paper show marked variations in 
the spectral quality and total intensity of the biologically 
effective component of ultraviolet solar radiation caused by 
slight changes in turbidity of the atmosphere and by the well- 
known cyclonic and seasonal changes in ozone concentration 
in the stratosphere—ozone being very opaque to ultraviolet 
solar radiation. It is further shown that at a sea-level station 
in mid-latitude (Washington, D. C.) on the clearest days, the 
ultraviolet intensity at the noon-hour in midwinter is only 
about one-sixth that of midsummer. 


A TUNGSTEN-IN-QUARTZ LAMP, AND ITS APPLICATIONS 
IN PHOTOELECTRIC RADIOMETRY. 


A paper by R. Stair and W. O. Smith of the Bureau of 
Plant Industry (RP1543) in the June Journal of Research 
gives constructional data of a tungsten-filament-in-quartz 
lamp, and discusses its applications to precision filter radi- 
ometry. The lamp is adapted for use in calibrating photo- 
tubes for spectral and total radiation sensitivity. In its 
applications, it supplements the spectroradiometer for cali- 
brating phototubes, while employing the tubes with the 
associated filters used in measuring the radiation under in- 
vestigation. 

The spectral output-radiation from such a lamp depends 
upon the temperature of the filament, the emissivity of 
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tungsten, and transmission of the fused quartz envelope. 
A table of relative black body radiation intensities has been 
calculated for the range of temperature from 2500° to 2900° K. 
in the wave-length interval 2,300 to 3,500 Angstroms. 


FLASH POINTS OF HEAVY-BODIED PAINTS. 


The Tagliabue (Tag) Closed Cup Method for determining 
flash points is prescribed in certain specifications for paints. 
The use of this method in tests of heavy-bodied paints, such 
as the camouflage paint covered by Specification T-1279 of 
the Corps of Engineers, U. S. Army, is open to question be- 
cause no stirring device is used to maintain an even distribu- 
tion of temperature, as is done in the Pensky- Martens tester. 

In an article which has been submitted for publication in 
American Paint Journal and in the Bulletin of the American 
Society for Testing Materials, David Busker reports that 
during tests of heavy paints at the Bureau, wide differences 
were noted in the results obtained with the two types of tester. 
In the Pensky-Martens tester, the heat was uniformly dis- 
tributed throughout the paint, whereas in the Tag Closed Cup 
tester, the temperatures were higher at the periphery than 
at the central point where the thermometer is inserted, and 
corresponded quite closely with the flash-point temperatures 
indicated in the Pensky-Martens tester. 


COMMERCIAL STANDARD FOR MINERAL WOOL. 


Commercial Standard CS105-43 for mineral wool, loose, 
granulated or felted, for use in low-temperature installations, 
is now available in printed form. It is designed to improve 
the quality of mineral wool and to set up installation require- 
ments to be followed in insulating refrigerated compartments 
or cold areas. 

Mineral wool is a substance composed of very fine inter- 
laced mineral fiber, having the appearance of loose wool or 
cotton. It is a fibrous, glass-like material composed prin- 
cipally of silicates of calcium and aluminum and is manu- 
factured from rock, slag, or glass in numerous localities 
throughout the United States. 
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This voluntary standard was proposed in April 1942 by 
the Industrial Mineral Wool Institute. Following several 
preliminary meetings a proposed commercial standard was 
prepared and submitted to manufacturers and to technical, 
distributor, and consumer organizations for their comment. 
The standard was then adjusted according to the composite 
recommendations of those concerned and circulated to the 
trade for written acceptance. It covers minimum physical 
and chemical requirements for mineral wool, as well as for 
auxiliary materials used in insulating cold areas. Detailed 
descriptions and diagrams are included in the pamphlet to 
illustrate the proper way to insulate various types of con- 
struction, as well as the minimum thickness of insulation 
required for various operating temperatures. 

Printed copies of CS105-43 are obtainable from the Super- 
intendent of Documents, Government Printing Office, Wash- 
ington, D. C., at 5 cents each. 
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THE FRANKLIN INSTITUTE. 


NEW MEMBERS OF THE FRANKLIN INSTITUTE 
ELECTED, DECEMBER 1, 1942, TO MAY 1, 1943. 


HONORARY MEMBERS. 


Carl Van Doren Dr. George Washington Dr. Harold C. Urey 
Pierce 
ACTIVE LIFE MEMBER (NON-RESIDENT). 

Mr. Sidney R. Lindboe 

SUSTAINING MEMBERS. 
Charles M. Alrich Mr. Frank Koons Mr. George C. Tatem 
Harry Alfred Mr. Harry M. Pflager Mr. William F. 

Eberhardt, Jr. Mr. Joseph P. Rainey Wagner, Sr. 

Lionel L. Jacobs Mr. John S. Smith Mr. C. J. Waterman 

Mr. Stanley Switlick 

ACTIVE MEMBERS WITH FAMILY PRIVILEGES. 
Edward A. Allison Mr. Walter C. Eberhardt Mr. Joseph P. Hoffner 
Samuel B. Mr. Abraham Ehrlich Mr. Frank J. Hoguet 
Amsterdam Mr. French McC. Mr. Oskar Huber 

J. Paul Austin, Jr. Emmons Mr. John F. Igoe 
Louis Back Mr. Harold N. Ewertz Mr. Fred Jessar 
Albert Bailey Mr. Samuel S. Feldgoise Mr. Nelson Kelley 
W. Heston Bailey Mr. Edward J. Fisher Mr. W. O. Kellogg 
Somerville N. Baker Mr. Wilmot Fleming Mr. John E. Kieffer 
Aaron Batchelor Mr. William N. Force Mr. William K. Kimbel 
Joseph Earl Borton = Mr. Frederick W. Fox Mr. Stanley C. Kirn 
Louis L. Borucki Mr. Philip Frank Mr. W. A. Kofke 
David Brein Mr. Alex A. Gettlin Mr. Charles La Scala 
Bennett K. Burris Mr. Joseph M. Goldberg Mr. Samuel Lashof 
C. Jack Cohan Mr. N. R. Guilbert, Jr. Mr. Horace M. Laurence 
Edward Corner Mr. Walter F. Graham Mr. John W. Lewis 
Paul Chrstos, Jr. Mr. W. Howard Green Mr. Albert E. Lipp, Jr. 
Ralph Cornell Mr. Theodore C. Hammel Mr. William Lu 
Jesse G. Croll Mr. Robert S. Mr. George A. Lyons 
Lloyd H. Daub Harrison, Jr. Mr. Patrick L. Mackey 
Michael Donnelly Mr. Jesse G. Haydock Mr. Charles F. May 
William H. Dorrance Mr. J. Albert Herz Mr. John J. McCarthy 
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Mr. J. A. McPeak 

Mr. Charles A. Miller 
Mr. William Moennig 
Mr. J. G. Morris 

Mr. Herbert S. Murphy 
Mr. George F. Nelson 
Mr. Carl F. Norberg 
Mr. Robert Odlen 

Mr. Harry K. Pape 
Mr. Joevert G. Parks 
Mr. Ernest Peacock 
Mr. Henry J. Perazzoli 
Mr. G. M. Pfundt 

Mr. George K. Porter 
Mr. Robert W. Radcliff 
Mr. Albert L. Reidl 
Mr. Philip F. Ries 


Miss Berte Abramson 

Mr. G. N. Ackley 

Mr. Robert Goss Adams 

Mr. Grover C. Albright 

Mr. Erwin H. Aldinger 

Mr. J. F. Alexander 

Mr. Frank B. Allen 

Mr. Paul Angstadt 

Mr. B. F. Antrim, Jr. 

Mr. E. Hulings Antrim 

Mr. J. Milton Armstrong 

Mr. William Arthur 

Miss Marguerite 
Aspinwall 

Mr. Edward Austin 

Mr. Joseph Axelrod 

Mr. Julius Axelrod 

Mr. Harry Bailis 

Mr. Philip Bailis 

Mr. Charles Bain 

Mr. Hermann F. Ball 

Mr. William H. Ballman 

Mr. William H. Barnes 

Mr. Harry Batt 

Mr. J. J. Baylson 

Mr. Charles B. Beck 

Mr. Lloyd Bergeson 

Mr. Joseph Berlin 


New MEMBERS. 


Mr. D. Rosenberg 

Mr. J. K. Rothman 

Mr. H. W. Rudoler 

Mr. John B. Rulon 

Mr. William W. 
Rutherford 

Mr. W. H. Sackman 

Mr. Albert Schnyder 

Mr. J. F. Scott 

Mrs. Lucille Seelig 

Mr. Carl J. Seydel 

Mr. Herbert L. Shallcross 

Mr. Franklin Sheble II 

Mr. William M. Sheehan 

Mr. Barton P. Sloane, Jr. 

Mr. Benjamin Small 

Mr. Dudley S. Smith 


ACTIVE MEMBERS. 


Mr. Joseph Berliner 
Mr. Benjamin P. Berman 
Mr. George W. Berry 
Mr. George E. Billings 
Mr. William Black 
Mr. Louis Bless 
Mr. L. J. Bloche 
Mr. Nathan Bloom 
Mr. Harry Bodenschatz 
Mr. Thomas J. Boerner 
Mr. Ernest Boghosian 
Mr. Harry H. Bonner 
Mr. Elmer H. Booz 
Mr. Ralph C. Borthwick 
Mr. Howard A. 

Bradley, Jr. 
Charles S. Brangan 
A. Balfour Brehman 
Mr. Charles W. Brelle 
Mr. Morris Brenner 
Miss Alice Brothers 
Mr. Harry F. Brown 
Mr. George M. Brunner 
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Mr. Stephen R. Buck 
Mr. George B. Bush 
Mr. William J. Bush 
Miss Louise Calloway 
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Victor C. Smith 
Felix Spatola, Jr. 
Samuel Spector 
Benjamin J. 

Spotts, Sr. 
R. C. Stearnes, Jr. 
William Swan 
David G. Swanson 
Raymond E. Tarpley 
William O. Tatcher 
Clarence Tolan, Jr. 
David M. Trilling 
F. Addison Urie 
Frank J. Valetti 
Ludwig Weber 
Albert W. Zackey 
Francis G. Zeiser, Sr. 


Charles Campbell 
John J. Campbell 
Joseph A. Carlin, Jr. 
Thomas A. Cashin 
Philip Cass 
F. C. Castelli 
Morris Chaikens 
Walter W. 
Chamberlain 
Rowland J. Cheney 
Daniel R. Cianciulli 
Ray S. Clarke 
Lloyd G. Clift 
John P. Coghlan 
Raymon Coghlan 
Samuel A. Cohan 
Benjamin Cohen 
Samuel Cohen 
William Cohen 
S. Z. Cole 
Robert C. Collins, Jr. 
Watson R. Collins 
Robert R. Colsher 
Joseph M. Comerford 
Ira L, Conkling 
T. A. Connelly 
Frank C. Connolly 
Michael Cortizas 
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Mr. James F. Couch 
Mr. William F. Covolus 
Mr. H. L. Cox 
Mr. Maurice A. Crews 
Mr. John P. Croasdale 
Mr. William J. Cronin 
Mr. Horace Cullen 
Mr. James E. Currie 
Mr. G. M. Curry 
Mr. William B. Curry 
Mr. Herman Z. Cutler 
Mr. Clarence R. 
Daugherty 
Charles M. Davis 
R. N. Davis 
William M. Davis 
Mr. W. H. S. Deckard 
Mr. John E. Deininger 
Mr. J. Carroll Deisher 
Mr. F. M. Delrose 
Mr. Americo A. 
DeMartinis 
Guido D’Emilio 
John C. Dettra 
Alexandro Di Bona 
Mr. W. J. Diederichs 
Mr. Harry A. Dieffenbach 
Mr. John E. Diem 
Mr. Joseph A. Dietrick 
Mr. David Di Prespi 
Mr. Stan J. Dirvin 
Mr. H. E. Doughty 
Mr. Orney E. Dunnum 
Mr. Rudolph H. Durnell 
Mr. Leon T. Eckell 
Mr. Alfred K. Edwards 
Mr. K. F. Ehmann 
Mr. Elmer T. Erb 
Mr. Attilio Esposito 
Mr. George M. Ewing 
Mr. James D. Faires 
Mr. Harry M. Farber 
Mr. James F. Feehery 
Mr. John Thorpe Feidt 
Mr. Samuel P. Felix 
Mr. Harry J. Ferguson 
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Mr. Joseph Ferraro 
Mr. Epentus L. Fetteroff 


New MEMBERS. 


Mr. Frank Feuse 

Mr. Thomas J. Finnerty 

Mr. Arthur M. Fischman 

Mr. Adam Fisher 

Mr. Leo Fisher 

Mr. Irvin J. Fishman 

Mr. Charles Flachbarth 

Mr. M. L. Flashman 

Mr. Nathaniel R. 
Fleischer 

Mr. David A. Fleming 

Mr. Norton D. Fleu 

Mr. H. G. Fliess 

Mr. Oliver K. Follweiler 

Mrs. Irene W. Ford 

Mr. Aug. C. Frank 

Mr. Edward Frank 

Mr. H. B. Frazer 

Mr. Frank J. Frei 

Mr. George H. Frew, Jr. 

Mr. G. William Friedrich 

Mr. B. H. Frishmuth 

Mr. Robert S. Fulton 

Mr. Samuel E. Fulton 

Mr. Joel Gahuse 

Mr. Charles I. Gallagher 

Mr. Hugh F. Gannon 

Mr. Fred H. Ganther 

Mr. Andrew S. Gardner 

Mr. L. R. Garretson 

Mr. William C. Gaskill 

Mr. B. E. Gerhring 

Mr. Charles Gerlach 

Mr. Richard G. Gerlach 

Mr. Elbridge F. Gerry 

Mr. George D. Gideon 

Mr. Julius D. Gilbert 

Mr. Harry Gilmore 

Mr. Joseph F. Glaser 

Mr. Thomas R. Glynn 

Mr. Samuel Goldberg 

Mr. Solomon Goldman 

Mr. Morris Goldstein 

Mr. G. F. Goodman 

Mr. Harry Gordon 

Mr. William R. Govett 
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W. F. Grau 
Rowland N. Grauch 
R. A. Gravdahl 
Alfred M. Green 
Louis B. Greenstein 
John Greenwood 
George H. Greiss 
Charles H. Greyer 
Elwood F. Groom 
Samuel Gross 
Samuel P. Gross 
David Guither 

F, Eldridge Guither 
Joseph Guss 

Albert M. Haas 
Harry L. Haeberle 
Samuel Haftel 
Daniel J. Haley 
Edward C. Hall 
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George Purdy Hall 
Charles D. Hamburg 
Albert J. Hamlin 
Elisha Hammonds 
J. Walter Hammonds 
B. J. Hansner 
Harry K. Hardcastle 
Josiah P. Harmer 
Charles S. Harper 
Robert L. Harper 
A. M. Harrington 
Edson S. Harris 
John H. Harrison 
Frederick E. Hartweg 
Herman Hassinger 
C. E. Hastings 
George K. Heebner 
C. Oscar Hegardt 
Alfred G. Heinemann 
Ralph E. Hellmer 
William M. 
Hendrickson 
Harry T. Henning 
John J. Henry 
H. C. Henssler, Jr. 
Robert S. Herbert 
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. Arthur S. Graham 
. Mario Grasso 


Mr. Wilmer J. Herrick 
Mr. Joseph F. Hess 
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Mr. Michael Hibner 
Mr. Newton Howard 
Holden 
Mr. Charles A. Holm 
Mr. Robert M. Holtzman 
Mr. Frank P. Homer 
Mr. Albert M. Hooven 
Mr. Charles A. Hooven 
Mr. Russell Hopkins 
Mr. Clifford F. Howard 
Mr. James Huey 
Mr. Ernest N. Hyde 
Mr. Theodore Irvin 
Mr. Albert Jackalous 
Mr. Leroy S. James 
Mr. Casper G. Jay 
Mr. John E. Johnson 
Mr. J. Carroll Johnston 
Mr. C. Walker Jones 
Mr. Charles D. Jones 
Miss D. E. Jones 
Mr. Edward Worthington 
Jones 
Mr. Whitney W. Jones 
Mr. Yern Jorgenson 
Mr. Alexander Kamens 
Mr. Joseph Katz 
Mr. William H. Kay 
Mr. Alfred R. Kennedy 
Mr. L. M. Kennedy 
Mr. Matthew Kenney 
Mr. Harry H. Kepley 
Mr. Max Kersun 
Mr. G. Walker Ketcham 
Mr. Joseph D. Kevorkian 
Mr. Carl Kirsch 
Mr. Walter Kiwi 
Mr. Charles Klause 
Mr. William 
Kleimenhagen 
Mr. J. Francis Kleinz 
Mr. Maurice Kline 
Mr. Harry N. Knauer 
Mr. Charles Knox, Jr. 
Mr. B. H. Brewster Koons 
Mr. Lewis A. Kopp 
Mr. Louis Korr 
Mr. J. Kosoy 


New MEMBERS. 


Mr. Abraham L. 
Krakovitz 
Mr. William C. Kroberger 
Mr. Alfred Kronschnable 
Mr. Henry J. Kunzig 
Mr. John M. Kupp 
Mr. Claude Lachman 
Mr. Charles Laessle, Jr. 
Mr. G. Duke Lambert 
Mr. Frederick William 
Lang 
Mr. Percy L. Lanning, Jr. 
Mr. Frederick B. Larsen 
Mr. George F. Lawrence 
Mr. Robert T. Lawson 
Mr. Sandor Lederer 
Mr. Stanley Lehmberg 
Mr. Robert Paul Leopold 
Mr. Robert T. Leopold 
Mr. Aaron Levine 
Mr. Benjamin Levy 
Mr. Martin M. Levy 
Mr. John Locke Lewis 
Mr. Martin Lewis 
Mr. Harry C. Lieberman 
Mr. Arthur H. Lockhart 
Mr. Isaac Lodge 
Mr. Joseph Loeffler 
Mr. Henry W. Long 
Mr. Fred Lorber 
Mrs. Ruth P. Loveland 
Mr. William A. 
MacAvoy, Jr. 
Mr. Robert MacMullen 
Mr. David Maimin 
Mr. Paul Maleson 
Mr. Thomas D. Marino 
Mr. Joseph Marziani 
Mr. Earl L. Mason 
Mr. Peter Masucci 
Mr. Joseph Mayer 
Mr. Robert McCurdy 
Mr. James F. McCormick 
Mr. Ralph D. Megee 
Mr. Walter E. Michaelis 
Mr. Samuel C. Milbourne 
Mr. J. A. Mitchell 
Mr. James H. Mitchell 
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Mr. Frank R. Moore 

Mr. J.Frank Motson 

Mr. J. N. Moyer 

Mr. Frank A. 
Muenzenberger 

Mr. A. P. Mynders 

Mr. Charles R. Nalle 

Mr. Martin Nebel 

Mr. Palmer Nelsen 

Mr. Walter A. Neuman 

Mr. W. Herbert Newbold 

Mr. Richard E. Newman 

Mr. Russell E. Nittinger 

Mr. William L. Noe 

Mr. Richard T. Orth 

Mr. F. C. Osterhout 

Mr. Charles A. Paoli 

Mr. Joseph W. Parkin, Jr. 

Mr. R. A. Parkinson 

Mrs. Robert W. 

Parkinson 

Mr. George H. Pechin 

Mr. H. F. Peirce 

Mr. Nathan Perilstein 

Mr. L. Franklin Peters 

Mr. W. J. Pickwell 

Mr. Arthur Platt 

Mr. Harold W. Polin 

Mr. Lawrence W. Post 

Mr. Lynn F. Powell 

Mr. Thomas J. Power, Jr. 

Mr. Dana Pratt 

Mr. F. V. Rais 

Mr. Willard M. Rambo 

Mr. Charles A. R. Rath 

Mr. R. Stewart Rawlings 

Mr. D. H. Redmond 

Mr. D. E. Regan 

Mr. Charles Reitzes 

Mr. J. Walter Rex 

Mr. Horace S. Rhoads 

Mr. J. Cecil Rhodes 

Mr. John H. Rice, 3rd 

Mr. F. W. Riess 

Mr. John P. Rife 

Mr. Albert Riggall 

Mr. John H. Ringe, 3rd 

Mr. Joseph Roberts 
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Mr. Archie M. Robertson 
Mr. David T. Robinson 
Mr. Francis Churchill 
Rogers 
Mr. Henry A. 
Romberger, Jr. 
Commander W. P. 
Roop, U.S.N. 
Mr. Benjamin Rose 
Mr. Morton M. Rose 
Mr. Charles Rosenthal 
Mr. Edward J. Rosser 
Mr. Howard T. Rossler 
Mr. Louis Roth 
Mr. Guy F. Rowlett 
Mr. Karl Rugart 
Mr. Warren E. Ruttle 
Mr. Reginald Ryberg 
Mr. Irvin Sachs 
Mr. George Sachsenmaier 
Mr. Ludwig Saha 
Mr. Harold B. Saler 
Mr. Samuel Saler 
Mr. Julius Saltzman 
Mr. Edgar M. 
Scattergood 
Mr. Fred C. Schadrack 
Mr. A. E. Schaffer 
Mr. H. V. Schaffer 
Mr. Irving Schaffer 
Mr. Morris Scherb 
Mr. John Schmidt 
Mr. D. F. Schmit 
Mr. Fred J. Schmitt 
Mr. Joseph Schmitz, Jr. 
Mr. Charles F. Schneider 
Mr. E. Schneider 
Mr. William L. Schneider 
Mr. Charles Schoenfeld 
Mr. W. Richison Schofield 
Mr. Herman Schulty 
Mr. David Schwartz 
Mr. Frank S. Schwarz 
Mr. Nicholas Sclufer 
Mr. David Scott, Jr. 
Mr. Donald C. Scott 
Mr. G. Frank H. 
Senderling 
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Mr. Paul J. Seng 

Mr. Milton C. Seib 

Mr. C. J. Seibel 

Mr. George Senn 

Mr. Bernard Shapiro 

Mr. M. M. Sharf 

Mr. Edgar Y. Sharlip 

Mr. F. Dickman Shaw 

Mr. Nathan Sherman 

Mr. L. H. Shingle 

Mr. W. L. Shirley 

Mrs. Stella D. Shober 

Mr. D. F. Showell 

Mr. Morris I. Shoyer 

Mr. George Siegel 

Mr. Charles Silver 

Mr. Maurice A. Singer 

Mr. William H. Simkins 

Mr. James A. Simmons 

Mr. Omer Simon 

Mr. E. A. Skowbo 

Mr. Albert J. Slap 

Mr. H. Royer Smith 

Mr. James J. Smith 

Mr. Joseph A. Smith 

Mr. Nicholas A. Smith 

Mr. Parker M. Smith 

Mr. Rudolph Smith 

Mr. Joseph J. Smyth 

Mr. Willard H. Snell 

Mr. A. D. Snyder 

Mr. Baer Soffe 

Mr. Norwood R. Solly 

Mr. K. Theodore 

Sorensen 

Mr. Morris J. Sostmann 

Mr. Meyer Sperber 

Mr. Leo J. Spillane 

Mr. William Spring 

Mr. Cameron F. Stains 

Mr. Howard B. ‘ 
Steigelman 

Mr. Michael Stein 

Mr. Seymour E. 
Steinbrook 

Mr. G. H. Stevens 

Mr. William H. Stilz 

Mr. John H. Stortz 
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. Joseph A. Stott 

. Albert H. Stover 

Charles W. Strong 

A. Fred Sturgis 

Jacob Sucknow 

W. E. Swanger 

A. W. Swartz 

L. H. Swind 

Joseph H. Thomas 

Eric E. Tinney 

J. B. Tomlinson 

J. Edward Topper 

Lew Trenner 

Julius J. Trumper 

Peter Tursi 

Paul S. Tyson 

Alfred Paul Urban 

Roland G. E. Ullman 

Daniel M. Vail 

J. T. Vallbrecht 

C. K. Van Rensselaer 

Edmond F. Venzie 

Frederick C. Vieser 

Joseph S. Vila 

Morris Vinskin 

Harry C. Viti 

Thomas A. Voelker 

John C. Vogtle 

H. G. von Cullin 

Conrad Wacker 

Charles A. Wagner 

Alfred S. Wall 

Curtis C. Wallace, Jr. 

William J. Walsh 

Joseph H. Walter 

Edgar S. Ward 

Martin L. 
Wasserman 

J. Hampton Webb 

Franklin Wedge 

David M. Weiner 

Morris Weiner 

Irvin I. Weinreich 

Philip H. Weinreich 

Morton B. Weinstein 

F. W. Weir 

Samuel Weisfield 

Harry Weissinger 
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Mr. Frank A. Wellman 
Mr. H. O. West 

Mr. S. Rogers Wetherill 
Mr. Clarence B. White 
Mr. James J. White 
Mr. R. R. Wiggins 

Mr. G. Philip Wild 

Mr. I. W. Wilenchik 
Mr. Charles Wilhelm 
Mr. James S. Wilkes 
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Mr. T. Carroll Willson, Jr. 
Mr. T. Carroll Willson, Sr. 
Mr. Alfred J. Williams 
Mr. J. R. Williams 

Mr. Charles D. Witman 
Mr. Clarence W. Wood 
Mr. Howard Woodward 
Mr. Earl Wunderlich 

Mr. Charles A. Yeager 
Mr. Alexander Younger 
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Mr. Aram Zakian 
Mr. Henry G. Zeiher 
Mr. Edward E. 
Zimmerman 
Mr. J. Milton 
Zimmerman 
Mr. John 
Zimmermann, Jr. 
Mr. George Zorn, Jr. 
Mr. Nathan Zucker 


ACTIVE NON-RESIDENT MEMBERS. 


Mr. Donald P. Campbell 

Mr. Randolph Church 

Mr. Thomas Harley 
Davis 

Mr. Frank S. Fawcett 

Mr. A. A. Green 

Mr. James G. Goodwillie 


Active Life Members. 


Mr. Charles Fearon 
Mr. Robert Moran 


Associate Life Members. 


Mr. G. H. Blakely 
Mr. William J. 

McCahan, Jr. 
Mr. George W. Morris 
Mrs. Arthur V. Morton 
Mr. William L. Nevin 
Mr. Frank A. Wills 
Mr. Wilson S. Yerger 


Mr. Isaac E. Johnson, 3rd 
Mr. Arthur H. Jordon 
Mr. Albert H. Klein 

Mr. Richard Kravits 

Mr. Joseph E. Mertle 
Mr. Carl N. Martin 

Mr. Wilbur E. Meserve 


NECROLOGY. 
Active Members. 


Mr. William H. Balls 
Mr. George E. Deming 
Mr. Thomas W. Evers 
Mrs. Jennie M. Fels 
Mr. Henry Reed Hatfield 
(Also Associate Life 
Member) 
Mr. Edwards Kneass 
Mr. Charles McDowell 
Mr. Charles A. Schranz 
Mr. Joseph F. Stockwell 
Mr. Walter A. Sweeting 
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Mr. Allan W. Morton 
Dr. Archie D. Power 
Mr. Arthur L. Stern 
Mr. L. E. Swedlund 
Mr. Harry L. Winston 
Mr. Howard Wright 


Mr. John P. Sykes 
Rear Admiral J. G. 
Tawresey 
Mr. Paul Thompson (also 
Associate Life Member) 
Mr. Charles W. Werst 


Active Non-Resident Members. 


Mr. Frederick M. Becket 
Dr. Alice G. Bryant 

Dr. Charles E. Caspari 
Mr. H. H. Vaughan 


The Committee on Library desires to add to the collections any technical 


works that members would wish to contribute. 
acknowledged and placed in the library. 


to other libraries as gifts of the donor. 
Photostat prints of any material in the collections can be supplied on request. 
The average cost for a print 9 X 14 inches is thirty-five cents. 
The library and the reading room are open on Mondays, Tuesdays, Wednes- 
days and Fridays from nine o'clock A.M. until five o’clock p.m., Thursdays from 
two until ten o'clock p.m. and Saturdays from nine until twelve o’clock noon. 


Contributions will be gratefully 
Duplicates received will be transferred 
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RECENT ADDITIONS. 
AERONAUTICS. 


Apams, Harotp W. Aircraft Hydraulics. 1943. 

AuMENT, HuGH C., Jr. Airplane Hydraulic Systems. 1943. 

GLAUERT, H. The Elements of Aerofoil and Airscrew Theory. 1943. 
Jorcu, WILLIAM F. Aircraft Electrical Systems. 1943. 

KINERT, REED. America’s Fighting Planes in Action. 1943. 

MILBANK, JEREMIAH, JR. The First Century of Flight in America. 1943. 


AGRICULTURE. 


Brown, Harry BATEs. Cotton. Second Edition. 1938. 
ASTRONOMY. 
GOLDBERG, LEO, AND LAWRENCE H. ALLER. Atoms, Stars and Nebule. 1943. 
BIOCHEMISTRY. 


HARRIS, ROBERT S., AND KENNETH VY. THIMANN, Editors. Vitamins and Hor- 
mones. Volume lr. 1943. 

PETERSON, WILLIAM H., JOHN T. SKINNER AND FRANK M. Stronc. Elements 
of Food Biochemistry. 1943. 


BIOGRAPHY. 


Hott, RACKHAM. George Washington Carver. 1943. 
KELLY, FRep C. The Wright Brothers. 1943. 
CHEMISTRY AND CHEMICAL TECHNOLOGY. 

American Institute of Chemical Engineers. Transactions. Volume 38, 1942. 
1943. 

BENNETT, H., Editor-in-Chief. The Chemical Formulary. Volume 6. 1943. 

BENNETT, H. Practical Emulsions. 1943. 

Chemical Rubber Company. Handbook of Chemistry and Physics. Twenty- 
seventh Edition. 1943. 

KALICHEVSKY, V. A. The Amazing Petroleum Industry. 1943. 

KILLIAN, J. Crystals. 1941. 

MIDLETON, ARTHUR R., AND JOHN W. WILLARD. Semimicro Qualitative Analy- 
sis. I940. 

OtneEy, Louis A. Textile Chemistry and Dyeing. Part 1: Chemical Tech- 
nology of the Fibers. Eighth Edition. 1942. 

Organic Syntheses. Volumes 22 and 23. Collective Volume 2. 1942-1943. 

SmiTH, PAuL I. Glue and Gelatine. 1943. 


DICTIONARIES AND ENCYCLOPEDIAS. 


Americana Annual. 1943. 
CRISPIN, FREDERIC SwING. Dictionary of Technical Terms. Revised Edition. 


1942. 
ELECTRIC ENGINEERING. 


FIELD, SAMUEL, AND A. DupLEY WEILL. Electroplating. Fourth Edition. 
1943. 
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Harvey, A. F. High Frequency Thermionic Tubes. 1943. 

Lincotn, E. S., Lincoln’s Industrial-Commercial Electrical Reference. First 
Edition. 1941. 

MacLaren, Matcotm. The Rise of the Electrical Industry during the Nine- 
teenth Century. 1943. 

ScHELKUNOFF, S. A. Electromagnetic Waves. 1943. 

ScHIEPE, EpwARD M. The Inductance Authority. 1942. 


ENGINEERING. 


BILLIG, Kurt. Pre-stressed Reinforced Concrete. Part 1. 1943. 
Hoo., GEORGE A., AND W. S. KINNE, Editors. Structural Members and Con- 
nections. Second Edition. 1943. 


PARKER, HARRY. Simplified Design of Reinforced Concrete. 1943. ‘ 
GEODESY. 

SHARP, H. OAKLEY. Geodetic Control Surveys. Second Edition. 1943. : 
GEOLOGY. ‘ 


GROVER, NATHAN CLIFFORD, AND ARTHUR WILLIAM HARRINGTON. Steam Flow 
Measurements, Records and their Uses. 1943. 
SVERDRUP, H. U., MARTIN W. JOHNSON, AND RICHARD H. FLEMING. The Oceans. 
1942. 
INDUSTRIAL MANAGEMENT. 


Atrorp, L. P. Cost and Production Handbook. 1942. 


MANUFACTURES. 


Hunt, Antony. Textile Design. No date. 
Ler, RutH WEBB. Sandwich Glass. Second Edition. 1939. 
MONAGHAN, JAMESF. Processing and Finishing Cottons. Volumes 1-2. 1935. 


MATHEMATICS. 


American Society for Testing Materials. A.S.T.M. Manual on Presentation of 
Data. 1943. F 
CraiG, HOMER VINCENT. Vector and Tensor Analysis. First Edition. 1943. ‘ 
MARGENAU, HENRY, AND GEORGE MosELEY Murpuy. The Mathematics of 
Physics and Chemistry. 1943. 
TircHMAkRSH, E. C. The Theory of Functions. Second Edition. 1939. 


MECHANICAL ENGINEERING. 
BILLINGS, J. HARLAND. Applied Kinematics. Second Edition. 1943. 
New Jersey Zinc Company. Die Casting for Engineers. 1942. 
METALLURGY. 


Haywarp, CARLE R. An Outline of Metallurgical Practice. Second Edition. 


1940. 

PAINTON, EpGAR T. The Working of Aluminium. 1927. 

PORTEVIN, ALBERT. Introduction to the Study of Heat Treatment of Metal- 
lurgical Products. 1939. 
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NAVAL ARCHITECTURE. 


CARMICHAEL, A. W. Practical Ship Production. Second Edition. 1941. 
EDDINGTON, WALTER J. Glossary of Shipbuilding and Outfitting Terms. 1943. 
HARMAN, F. WARD. Ship Models Illustrated. 1941. 

RIMINGTON, CRITCHELL. Fighting Fleets. 1943 Edition. 


PHYSICS. 


DINGLE, HERBERT. Sub-Atomic Physics. 1943. 

Jameson, ALEX. H. An Introduction to Fluid Mechanics. 1942. 

LeMON, HARVEY BRACE, AND MICHAEL FERENCE, JR. Analytical Experimental 
Physics. 1943. 


RAILROAD ENGINEERING. 
Railway Engineering and Maintenance Cyclopedia. Fifth Edition. 1942. 
VAN NAME, WILLARD G. Collection of Photographs of Old Locomotives and 
Cars, Taken before the Year, 1898. 


SANITARY ENGINEERING. 


American Foundrymen’s Association. Tentative Recommended Good Practice 
Code and Handbook on the Fundamentals of Design, Construction, Opera- 
tion and Maintenance of Exhaust Systems. 1938. 


NOTES FROM THE BARTOL RESEARCH FOUNDATION. 


NOTE ON THE NATURE OF THE GAS MIXTURE IN 
SELF-QUENCHING GEIGER-MUELLER TUBES.* 


BY 
PAUL WEISZ 


From studies on the breakdown of spark gaps, and from 
theoretical and experimental work on the discharge process 
of Geiger-Mueller tubes itself, we conclude that the intrinsic 
difference between self-quenching and _ non-self-quenching 
G-M (Geiger-Mueller) tubes lies in the particular nature of 
the action of ultraviolet photons in each case. 

In the G-M tube which is not self-quenching, a great 
number of photons are emitted during the discharge process. 
In so-called self-quenching G-M tubes, no evidence has been 
found for the presence of photons at appreciable distances 
from the wire. Previous experimental evidence,! and a find- 
ing to be reported, seem to suggest that ultraviolet light 
quanta are still produced, but are quickly absorbed by the 
‘quenching gas,’’ probably within a distance of a few mean- 
free paths in the gas. 

Thus it is the purpose of the quenching gas to absorb all 
ultraviolet quanta produced, or ultimately to make the prob- 
ability of production of negative ions by way of photons and 
photoelectrons as nearly as possible equal to zero. Since 
this probability will depend on the density of the quenching 
gas, an investigation was made of the quenching properties 
as a function of the partial pressure of the quenching gas. 
This was accomplished by studying the width of the “ pla- 
teau,’’ for mixtures of argon as the ‘‘multiplication gas’”’ and 
ethyl-ether as the quenching gas. The plateau is defined as 


* Reprinted from an article published in full in The Physical Review, 62, 


477 (1942). 
1W.E. Ramsey, Phys. Rev., 61, 96-97 (1942); H. G. Stever, Phys. Rev., 59, 


765 (1941), and Phys. Rev., 61, 38-52 (1942). 
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that region of applied electrode potentials which lies between 
two limits A and B, such that beginning at the potential A 
all pulses are of a convenient pulse size (three volts in this 
case), and that the tube starts to show multiple discharges 
when the potential is increased beyond B, and the tube is 
operated with a normally small external resistance of I meg- 
ohm. The appearance of such multiple pulses can be in- 
terpreted in the light of the picture of the mechanism as 
being evidence for the production of photoelectrons, their 
subsequent capture and release, and thereby the repetition 
of the discharge process. 
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Fic. 1. The plateau region for the argon/ethyl-ether counter. The counter 
is }’’X3"’. The central wire is 0.003’, -he total pressure, 6 cm of Hg. 


Figure I shows the results of the plateau measurements 
for a tube containing argon/ethyl-ether at 6 cm. total pres- 
sure; the partial ether pressure is experssed in percentage of 
total pressure. It can be seen that for this total pressure, an 
ether content of about 20 to 25 per cent. yielded the optimum 
plateau. At relatively high ether contents, the pulses were 
no longer of equal amplitudes but showed appreciable dis- 
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tribution in size up to potentials indicated by curve C, al- 
though the average pulse size was still approximately 3 volts 
at A. In the light of the picture outlined the following ex- 
planation suggests itself: At great densities of the quenching 
gas the absorption of photons becomes so intense that even 
the propagation of the mechanism along the wire cannot 
always take place along its entire length. The pulse size, 
being proportional to the total amount of charge produced, 
may then be nearly proportional to the length along which 
each discharge is able to propagate before it fails in producing 
at least one more photoelectron. The pulse size may thus 
be subject to appreciable fluctuations. 

Similar measurements concerning optimum plateau con- 
ditions were made at different total pressures, and the fol- 


lowing data obtained: 
Approx. % ether 


Total pressure for opt. plateau 


6 cm. 20 to 25% 

od /¢ 
15 cm. 7 to 10% 
40 cm. approx. 3% 


It appears that the absolute pressure or density of the quench- 
ing gas alone determines the quenching characteristics re- 
gardless of what the total gas pressure may be, the partial 
ether pressure for optimum conditions being about 1.5 cm. 
of mercury. It seems noteworthy that this finding seems to 
favor the concept of photons being absorbed by the presence 
of a sufficient amount of the quenching gas, in preference to 
a picture of the mechanism involving the assumption of some 
action by which this gas prevents the creation of quanta 
during the ordinary discharge process. 


NOTES FROM THE BIOCHEMICAL RESEARCH 
FOUNDATION. 


Effect of Propionic Acid and of Na Propionate on Para- 
mecia.—HuGu E. Potts. Paramecium has been used as a 
test animal for various chemicals (Ferguson, Frederick F., 
Holmes, Jane R., and Lavor, Edward, Journal of the Elisha 
Mitchell Scientific Society, 58: 53, 1942; Feiler, Marie, Archiv 
fuer Protistenkunde, 59: 562, 1927; Mel’nikov, N. N., Avetes- 
yan, A. M., and Rokitskaya, M. S., Comptes rendus de 
l’académie des sciences de l’U. R. S. S., 31: 123, 1941). It has 
the advantages of being large enough for observation of 
changes in size, shape, contractile vacuole activity, and of 
being easy to cultivate in large numbers. 

Since propionic acid in minute traces (Miller, F. W., Jr., 
Agricultural News Letter, Public Relations Department, E. I. 
du Pont de Nemours and Co., 10: 62, 1942) is capable of in- 
hibiting growth of mold on bread, it was thought that it 
might prove to be of some value as a therapeutic agent. 
Therefore it was tested in different concentrations on Para- 
mectum caudatum. Small watch glasses were used for the 
experiments. These measure about 20 mm. in diameter and 
can conveniently hold 1 cc. or a little more of fluid. The 
bottom of each glass is grooved and the upper rims are bevelled 
and ground so that they can be stacked one upon another, 
thus avoiding any great amount of evaporation. 

Paramecia were cultured in small glass jars. The pH 
was kept at about 7.0 by addition of NaOH or NaHCQs. 
An original inoculation of Serratia marcescens was made to 
the paramecia in aquarium water, and the organisms were 
fed on a few flakes of dried rolled oats per jar twice a week. 
Subcultures were made about once a month by putting sev- 
eral cc. from old cultures into clean jars and adding distilled 
water, a few cc. per day until the jars were filled. 

For testing purposes paramecia were pipetted into the 
watch glasses and the desired amount of propionic acid added, 
proportions being reached by counting the number of drops 
of culture and of chemical, the same pipette being used for 
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both. Pipettes were rinsed thoroughly with distilled water 
whenever necessary. 

The propionic acid used was acid propionic pure, pur- 
chased from The Coleman and Bell Co., Norwood, Ohio. 
This was titrated in our laboratory and found to be 13.13 
normal. It was diluted into a series of stock solutions, each 
1/10 the concentration of the previous one, down to 13.13 
x 107° normal. From each stock solution 10, 8, 6, 4, 2 and 
o drops of the acid were added to 10, 12, 14, 16, 18 and 20 
drops respectively of paramecium culture. After addition of 
the acid each watch glass was tilted and agitated to mix the 
acid evenly throughout the culture. 

The paramecia were then examined under a magnification 
of 40X, through a dissecting microscope at half-hourly, hourly 
or greater intervals, as seemed advisable. The pH was ap- 
proximated at the beginning or shortly after the beginning of 
each experiment by use of bromthymol blue. 

All concentrations down to 2.63 X 107° N killed all para- 
mecia in 3 hours. The pH on all these concentrations of 
propionic acid was 6.0 or less, the 2.63 X 107° N and 1.31 
x 107° N being approximately 6.0. With dilutions beyond 
1.31 X 10°* N the pH rose gradually, being approximately 
7.0 at 2.63 X 194 N, and none of these had any appreciable 
effect. 

In order to determine whether or not it was merely the 
acid effect which killed the paramecia, the propionic acid was 
~eutralized with normal NaOH, 2 cc. of the alkali being added 
i 10.5 cc. of acid. The resulting solution had a molarity of 
4. 3 X 107%. 10, 8, 6, 4, 2 and o drops of the neutralized 
aci. were added respectively to 10, 12, 14, 16, 18 and 20 drops 
of paramecia culture. The pH of all concentrations was 
about 7.1. 5.51 X 10°? M killed all paramecia in 24 hours. 
4.41 and 3.31 X 10° M killed a few in 24 hours. Greater 
dilutions had no noticeable effect. 

These results indicate that both the acid and its sodium 
salt have an appreciable action on paramecia. The acid, 
however, is many times as active as its sodium salt. It is 
probable that the sudden drop in pH, on addition of pro- 
pionic acid, is in itself capable of causing the death of 
paramecia. 
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BOOK REVIEWS. 


ATOMS, STARS AND NEBUL#&, by Leo Goldberg and Lawrence H. Aller. (The 
Harvard Books on Astronomy.) 323 pages, 150 illustrations. Philadelphia, 
The Blakiston Company, 1943. Price $2.50. 

The fifth of the Harvard Books on Astronomy has now appeared, and will 
soon find its place on the shelf of most-used books of astronomers and physicists. 
For while the book is intended for the serious non-professional reader, it has so 
summarized the present status of astronomical spectroscopy that astrophysicists 
will find it constantly useful. 

The narrative style is maintained in excellent fashion through the whole 
work, and as in some of the other books of the series the photographs of con- 
temporary contributors to advances in the field add a worthwhile touch. Tables 
summarizing the data on star temperatures, masses, densities and other charac- 
teristics are scattered liberally, and every other page contains an illuminating 
photograph of a spectrum, a graph, or a diagram. 

Supplemented with readings in the current literature, this book fills the need 
for a modern introductory textbook in astrophysics, yet it can be read with 
pleasure by those who are merely inquisitive about the methods and findings of 


modern astronomical spectroscopy. 
Roy K. MARSHALL. 


An INTRODUCTION TO FLUID MEcHANICS, by Alexander H. Jameson, M.Sc., 
M. Inst. C. E., F. K.C. 245 pages, illustrations, 14 X 22 cms. New York, 
Longmans, Green and Company. 1942. Price $3.40/ 

Although the principles of fluid mechanics have been applied in various 
fields for a good many years, it has only been comparatively recently that applica- 
tions have been multiplied to a great extent. Perhaps the most common peace- 
time use today is the stream-lined automobile and only a short time before the: 
appeared, was there such broad application of fluid mechanics. The Unive - y 
of London decided that an elementary knowledge of this fundamental s ct 
shall be required of all candidates in the Part I, B.Sc. (Engineering) Exam it. a, 
and this book is intended for such second year students. 

The treatment is, of course, of a rigorous mathematical nature. The method 
of presentation, the logical progression, and the somewhat unusual approach in 
the way of demonstrating the solution of examples in the text and the means 
indicated in the exercises at the end of chapters all make for a lightening and mor 
impressive text. The book is divided into eight chapters. They begin with 
an explanation of liquids and gases, and follow through with pressures, th_ per- 
fect fluid, motion of a perfect fluid and motion of a cylinder and sphere through a 
fluid, flow through orifices, flow in pipes, channels and notches, and the measure- 
ment of velocity. Appendixes contain Bernoulli's Theorem for Compressible 
Fluids, and alternative treatment in the method of dimensions, and the derivation 


of exponential flow formulae. A subject index is in the back. 
R. H. OPPERMANN. 
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THe ELEMENTS OF AEROFOIL AND AIRSCREW THEORY, by H. Glauert. 228 pages, 
illustrations, 14 X 22 cms. New York, The Macmillan Company, 1943. 
Price $3.50. 

The actions of an aerofoil are not new to those thousands who have been so 
fortunate as to have observed them in the museum of The Franklin Institute. 
That they are of importance is unquestioned as is the theory behind them. The 
main object of aerofoil theory is to explain and to predict the lift and drag ex- 
perienced by an aerofoil. Much of this has been done, as is shown in the exhibi- 
tion and in everyday practice. A satisfactory theory has been developed in 
recent years for the ordinary range below the critical angle. The determination 
of the maximum lift of an aerofoil and of the critical angle at which it occurs is 
not yet possible, although some insight has been obtained into the cause of the 
phenomenon. 

This is an American edition of the British book. It aims to give an account 
of this mathematical subject in terms that may be understood by those not having 
previous knowledge of hydrodynamics and by minimizing complex mathematical 
analysis. After a brief introduction, Bernoulli’s equation is discussed followed 
by an analysis in stream function which led to the determination of the theoretical 
flow past a cylinder in a uniform stream. Another type of flow in which the 
fluid circulates round the cylinder is next taken up and finally the velocity po- 
tential and the potential function. All this constitutes an introduction to those 
aspects of hydrodynamics which are required for the development of aerofoil 
theory. The following chapters deal successively with the lift of an aerofoil in 
two dimensional motion and an examination of the fundamental basis of this 
theory of flow governed by Joukowski’s hypothesis, the aerofoil in three dimen- 
sions, the monoplane aerofoil, and biplane aerofoils. The last three chapters 
deal with the development of the airscrew theory. This treatment begins with 
the momentum theory and goes through an investigation of the forces experi- 
enced by airscrew blades regarding each blade as an aerofoil element moving in 
its appropriate manner, and then wind tunnel interference. 

The book gives a rather complete treatment within its limitations. It 
forms a workable basis of knowledge of the subject. 

R. H. OPPERMANN. 


THe First CENTURY OF FLIGHT IN AMERICA, by Jeremiah Milbank, Jr. 248 
pages, illustrations, 15 X 22 cms. New Jersey, Princeton University Press, 
1943. Price $2.75. 

New industrial developments or inventions when presented for use or con- 
sumption by the public are seldom recognized as being the result of evolution 
over a long period. Without such a process there would be few of the many 
things available today that make up the tempo and ease of modern life. Aviation 
is certainly one of these and the book at hand is a means of bringing this realiza- 
tion. Mr. Milbank here reveals much that is generally unknown about flight in 
America for the nineteenth century. His account gives details of experiences of 
great inventors and aeronauts who provided most spectacular feats for early 
America, and at the same time laid the foundations for one of the greatest wonders 


of our present generation. 
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The first authenticated ascent in America was made June 24, 1784, by a boy 
of only 13 years, Edward Warren, in a balloon. Almost ten years later Jean 
Pierre Blanchard made the first actual aerial voyage in this country. George 
Washington was Blanchard’s patron, and the voyage lasted 46 minutes from 
Philadelphia to Gloucester County, N. J. These and other events are vividly 
described by the author in more or less chronological order. Special attention is 
given to the events occasioned by outstanding men such as John Wise and others. 
An entire section is devoted to T. S. C. Lowe and the Civil War and it is interest- 
ing to note here that in 1859 Lowe had decided to attempt flight in a balloon over 
the Atlantic Ocean. The voyage was scheduled to start from New York but the 
New York Gas Company was unable to deliver hydrogen fast enough to effect 
full inflation of the balloon. Dr. John C. Cresson, President of the Franklin 
Institute and of the Point Breeze Gas Works, of Philadelphia, then invited Lowe 
to his city where he guaranteed sufficient speed of inflation. Weather condi- 
tions postponed the trip, however, temporarily and later postponement was made 
final due to the Civil War. 

The Civil War offered the first demonstration of the military value of the 
balloon. Lowe organized a balloon division in the Union Army and did valuable 
work which was recognized by Federals and Confederates alike. Major General 
A. W. Greely, Chief of Army Signal Corps, stated ‘‘it may be safely claimed that 
the Union Army was saved from destruction at the battle of Fair Oaks, May 31 
to June 1, 1862, by the frequent and accurate reports of Professor Lowe.” Not 
only did the Confederates try to eliminate Lowe’s aeronautical corps but they 
even tried to organize their own balloon outfit. 

The book is a very interesting account that furnishes a background for those 
in aviation today and a fund of information for the layman. 

R. H. OPPERMANN. 


THE AMAZING PETROLEUM INbDusTRY, by V. A. Kalichevsky. 234 pages, dia- 
grams, 13 X 19 cms. New York, Reinhold Publishing Corporation, 1943. 
Price $2.25. 

There is a definite place in the literature for books of an informative char- 
acter for general reading. Whenever a field of human endeavor becomes promi- 
nent enough to affect the lives of large numbers of people, there is a duty, educa- 
tional if not otherwise, to make available to the people a means of understanding 
the processes involved. The petroleum industry is such a field. Its position as 
one of the foremost branches of chemistry is undisputed. Measured in tonnage 
it is probably the country’s largest supply of chemical raw materials. Measured 
in achievement, it ranks second to no other branch of chemistry as a supporter 
of the war effort, both in its role as a motor fuel and lubricant source, and as the 
actual savior of military transport through its basic contribution to the synthetic 
rubber program. Dr. Kalichevsky’s book, therefore, is timely. 

The book embraces the entire industry including the production and nature 
of petroleum, its transportation and storage, and the processing of petroleum. 
Special sections are devoted to petroleum emulsions and natural gas, gasoline 
and crude oil distillations, cracking, hydrogenation and aviation gasolines, chemi- 
cal treatment of light petroleum products, and chemicals from petroleum. The 
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book is a non-technical resumé—understandable by anyone. In the back there 
is a glossary of terms used in the industry. 

Readers of all kinds can profit by reading this book. Coming from a man 
closely connected with the industry for many years it is an authoritative explana- 
tion. He rightly states that it ‘‘should be of assistance to college students, to 
young engineers beginning their careers in the petroleum industry, and to persons 
already connected with it but not directly engaged in production or refining of oil. 

R. H. OPPERMANN. 


VECTOR AND TENSOR ANALYsIs, by Homer Vincent Craig. 434 pages, diagrams, 
16 X 23cms. New York, McGraw-Hill Book Company, 1943. Price $3.50. 


In the preface of this book the author refers to the economy of time in the 
training of engineers and physicists. While of course, this has been an important 
item before educators for years, it is refreshing to see it brought out so boldly 
when considering a subject for study, especially so fundamental a subject as con- 
cerns this book. Vector analysis is undoubtedly the best introduction to certain 
related branches of mathematics. Every student of theoretical physics, en- 
gineering, and differential geometry should have at least a reading knowledge of 
this subject. As to tensor analysis, this is a somewhat newer study and at present 
it can be said to be in a state of rapid expansion. The state of development is 
such that a broad and voluminous content is available. Nevertheless the book 
at hand undertakes to give an introduction to these subjects, sufficient to indicate 
their value as tools for the solution of practical problems and their flexibility of 
application. 

The book is divided into three parts the first of which is devoted to a mathe- 
matical background necessary for the main and later divisions. Here reference is 
made to the calculus of functions of two variables, parameterized arcs, deter- 
minants, the e-systems, and certain transformation equations. This takes up 
roughly, one-quarter of the content of the book and prepares the way for the study 
of elementary vector analysis which is the second division. It begins with a 
treatment on the algebra of vectors containing material basic for the examination 
of the simple N-dimensional space. Subsequently the differential calculus of 
vectors and integral transformations are taken up, followed by Stoke's theorem, 
proving for certain parameterized surfaces a modified form of the second of the 
two fundamental integral transformations. 

The third division of the book is devoted to the introduction of tensors and 
the explanation of extensions and the differential calculus of tensors. Division 
number four provides some material for practice reading in vector notation and 
covers classical dynamics based on the notion of absolute space and time, special 
relativity characterized by a three-dimensional Euclidean space and other con- 
cepts, and concluding with a formulation of the most important assumptions ol! 
general relativity. To this division is appended a bibliography. 

The construction of the book is one of logical progression with careful selec- 
tion of topics. Within limitations it is thorough and provides a course of study 
for practical application. 


R. H. OprpeERMANN. 
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ANALYTICAL EXPERIMENTAL Puysics, by Harvey Brace Lemon and Michael 
Ference, Jr. 584 pages, plates and diagrams, 23 X 30 cms. Chicago, The 
University of Chicago Press. 1943. Price $7.00. Second Printing $5.75. 


In view of the many textbooks on physics and laboratory manuals which 
are continually appearing, it seems strange that little mention is made of con- 
centrated study of the subject of laboratory practice in a broad light. For in- 
stance, most of the laboratory procedures are in effect standardized, a routine 
performance without much attention directed to the training received compared 
to hours expended, and methods for improving this. Also, aside from the dia- 
grams and tabulated results, there is little upon which can be based reflection of 
actions and functions during the experiment. A review and subsequent rearrange- 
ment of these as well as other details, incorporating any new devices, and co- 
ordination with an adequate text, will provide a course in physics which has out- 
standing merit. This is the thought behind the book at hand. Needless to say, 
it has many unusual features. 

The presentation leans toward the narrative type wherever possible. The 
sequence of topics is mechanics, heat, electricity and magnetism, wave motion, 
sound and light. There is a mathematical appendix in the back which contains 
all the formulas needed. The story of each topic is so told that the experiment is 
woven into it, making it an integral part. The authors use the novel method of 
enlarged reproductions of strips of motion picture film of many of the experi- 
ments in action and under observation. These are photographed demonstration 
experiments. By running over the successive frames of the motion picture repro- 
ductions the student sees what is happening; he sees the measuring instruments as 
they are in position, making quantitative records. The readings of these instru- 
ments on the original films when too small for satisfactory reproduction by the 
half tone processes are set in type alongside the pictures, together with the neces- 
sary mathematical operations required to reveal the law that is dictated by the 
phenomenon. Insofar as possible, the experiments are set up in such a manner 
as to provide two methods of approach to the conclusion, and thus independent 
check of the results is available. Of course not all the work is susceptible to the 
photographic method and, too, all may not agree on the value in the method that 
shouldbe weighed carefully by educators. 

The text has been for nearly eight years at the University of Chicago in the 
sophomore sequence courses. It is stated that about twenty different instructors 
of a great variety of training and talent have been free to criticize, suggest altera- 
tions, and discuss results; and rewritings have been made to meet these suggestions. 

R. H. OPPERMANN. 
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PUBLICATIONS RECEIVED. 


Basic Physics for Pilots and Flight Crews, by E. ]. Knapp, Ph.D. 118 pages, 
charts and diagrams, 13 X 19 cms. New York, Prentice-Hall, Inc., 1943. 
Price $1.25. 

College Physics, Abridged, Revised Edition, by Henry A. Perkins, Sc.D. 
593 pages, diagrams, 15 X 23 cms. New York, Prentice-Hall, Inc., 1943. Price 
$4.00. 

The Rise of the Electrical Industry During the Nineteenth Century, by Malcolm 
MacLaren. 225 pages, illustrations, 16 X 24cms. New Jersey, Princeton Uni- 
versity Press, 1943. Price $3.75. 

Basic Electricity for Communications, by William H. Timble. 603 pages, 
illustrations and diagrams, 15 X 22cms. New York, John Wiley and Sons, Inc., 
1943. Price $3.50. 

The Chemistry of the Aliphatic Orthoesters, by Howard W. Post. 188 pages, 
tables, 16 X 23 cms. New York, Reinhold Publishing Corporation, 1943. 
Price $4.00. 

Radio Troubleshooter’s Handbook, by Alfred A. Ghirardi. 744 pages, il- 
lustrations, 22 X 29 cms. New York, Radio & Technical Publishing Co., 1943. 


Price $5.00. 


CURRENT TOPICS. 


Termite Research.—Research for the last 30 years—particularly 
since 1939—has enabled the Department of Agriculture’s Research 
Administration to assist in developing and applying measures to 
protect buildings in military cantonments and in Federal housing 
projects against termites—among the worst insects hampering this 
country’s war effort. Suggested improvements in plans for these 
projects include provision for adequate ventilation and clearance of 
wood debris beneath buildings. These measures have largely super- 
seded the metal shields formerly used to protect buildings from 
termite damage. At the request of military authorities, Depart- 
ment entomologists also are seeking measures to prevent termite 
damage to military buildings and supplies in tropical and subtropical 
lands. Under way is a study to determine which of the many 
tropical woods are most resistant to termite attack. As facilities 
are not available for treating these woods with preservative chemi- 
cals in the areas where they are needed for construction, it is most 
important that the wood to be cut and used locally without treat- 
ment should have as much natural durability as possible. Research 
has shown that the life of most untreated woods used in contact 
with the soil, such as in temporary buildings and bridges, can be 
prolonged by poisoning the adjacent soil. Tests are in progress to 
determine which chemicals do this job best under tropical conditions. 
Creosote, pentachlorophenol, and some of the arsenical compounds 
are effective for the purpose in the United States. Department 
entomologists believe that tents, tarpaulins, clothing, and other 
materials that must be stored or used where they are subject to 
damage from termites, especially in the tropics, can be protected by 
chemical treatments at the time of their manufacture or use. Re- 
sults of limited tests—the only ones yet made—indicate that some 
of the chemicals that prevent mildewing and decay of many fabrics 
may also prevent termite damage. 


ae Me 


Largest Steel Chain Ever Made in Canada Now in Service.— 
(Mining and Metallurgy, Vol. 24, No. 437.) This was made by 
Sorel Steel Foundries Ltd. for use on a marine railway in a new ship 
building plant. This plant, situated on a comparatively narrow 
river, builds 10,000 ton ships which are winched into place on the 
launching platform of the railroad; the railroad is then slowly 
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lowered into the water, thus floating the ship’s hull. The chain, 
which was to be used as the restraining medium in lowering the 
platform with the 10,000 ton hull upon it, was of course an important 
factor in the plant. The job required 1} miles of 2} in. chain, each 
link weighing 38 Ib., and cast into 90 ft. lengths. One of the many 
difficulties encountered had to do with the proper testing of the 
chain to ensure that each length would properly meet the drag it 
would be subjected to when letting the 10,000 ton hull and launching 
platform down the inclined marine railway. An ingenious solution 
was devised by Sorel engineers, who had eight lugs welded around 
the sides of a 180 ton scow. Hawsers were run from each lug to the 
center of the scow and attached in turn to each individual go ft. 
length of chain, the other end of which was hooked to a shear leg. 
The shear leg was then operated, lifting against the scow, a total 
dead weight of 125 tons, and the scow was raised a minimum of 27 
in. on the draft line, which exerted a strain on the chain beyond the 
maximum called for in the specifications. Out of the 8400 links 
tested, only nine were rejected or less than 75 of one per cent. of 
the total production, an outstanding record of foundry control under 


the circumstances. 
R. H. O. 


Dry Ice Kills Coal Pile Fires.—(Heating and Ventilating, Vol. 40, 
No. 4.) A new and ingenious use for dry ice as an extinguisher of 
coal fires has been developed recently. An estimated 80,000,000 
tons of bituminous coal is now in storage throughout the country, 
a figure double the usual season’s reserve. This likewise doubles the 
risk of loss through spontaneous combustion, notwithstanding the 
protective effect of modern methods of storage. The dry ice is 
applied directly at the seat of the fire, through 15 ft., 3 in., driven 
pipes fitted with perforated tips. Charges of 50 lbs. of the solid 
CO, are placed in the pipes. Heat liberates the gas to blanket the 
hot zone and smother the smoldering coal, while the chilling action 
of the dry ice reduces the local temperature, this effectively elimi- 
nating the hot spot. To locate such danger points in the pile a hot 
spot indicator, consisting of an iron pipe containing a spring loaded 
signal, held back by a chain with fusible links at one foot intervals 
has been made available by the same manufacturer. Only the hot 
test link lets go, permitting the signal to rise, and this occurs only 
when temperature rises to 150 deg. F. Thus the location and depth 
of the danger point are advertised 72 hours before combustion can 
take place, at the normal rate of temperature rise for the average 
coal composition. 
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Two Traveling Tower Derricks Erect Timber Blimp Hangar. 
—(Engineering News-Record, Vol. 130, No. 11.) Steel Erector’s 
tower-mounted derricks are used to raise pre-assembled timber 
hangar sections for blimp housing at the Lakehurst, N. J., Naval Air 
Station. The hangars have a span of 246 ft. and a height above the 
ground of more than 170 ft. The timber arch trusses are set 20 ft. 
C. to C., 51 of them being used to include a length of 1000 ft. Over- 
all measurements of the hangar are 1058 ft. by 296.5 ft., the added 
length being taken up by door assemblies at both ends of the struc- 
ture. Timbers are preassembled on the ground into four sections 
for each arch truss. Two similar sections are framed into a “‘brace- 
bay”’ 20 ft. wide that forms the unit to be handled. The bottom 
section of the truss is eight panels long and the top portion is the 
remaining ten panels to the center line of the structure. The 
bottom sections of the arches are set in place in braced pairs by the 
individual derricks. The top quarters are then joined on the ground 
by the two derricks working together and raised to a connection 
with the first erected part to form the complete ring. Timber 
members for the arches are accurately precut to template in mills on 
the West Coast and are given a pressure-cylinder impregnation with 
chemicals to fire-proof and preserve them. The fire-proofed wood 
is required to stand a ? hour exposure to a standard time-temperature 
test in which the wood is exposed to gas flames of 1300 deg. at } 
hour, 1550 deg. at 3 hour, and 1658 deg. at ? hour. The wood will 
char at these temperatures but will not support combustion. Prog- 
ress of the work of erection when weather conditions have not been 
too adverse has been a complete 40 ft. advance in each working day. 
This requires that a braced-bay of the lower } of the arch be set, 
followed by connection and erection of the two top quarters. Inter- 
mediate bracing then is immediately placed to stabilize the section 
just set. Winds with a 5 minute average velocity of 67 knots 
(80 m.p.h.) have been successfully withstood. 

R. H. O. 


Invisible ‘* Raincoat ’”» Waterproofs Many Kinds of Materials.— 
An invisible ‘‘raincoat’’ which can be formed on cloth, paper and 
many other materials by exposing them to chemical vapors, thereby 
making them water-repellent, has been developed in the General 
Electric Research Laboratory, Dr. William D. Coolidge, G. E. Vice 
President and Director of Research, has announced recently. Many 
possible uses of this new method of waterproofing are now being 
studied. One of the most important so far is the treatment of 
ceramic insulators for radio equipment used by the armed forces of 
the United States. Such insulators lose much of their electrical 
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resistance when exposed to high humidity. This causes formation 
of a thin film of moisture on their surfaces and allows the currents 
to leak away. It has been the practice to coat them with wax, but 
the new chemical treatment is about nine times as effective, and its 
results are permanent. The process was discovered by Dr. Winton 
I. Patnode, G. E. research chemist. Various chemicals are used, one 
class being known as methyl chlor silanes. In a closed cabinet, the 
articles to be treated are exposed for a few minutes to vapors from 
such a compound. Then they are taken out and, if necessary, are 
exposed to ammonia vapor. This is to neutralize corrosive acids 
which may collect during treatment. Dr. Patnode is not able to 
explain exactly what happens in the process, but the result is that 
an extremely thin film, which resists wetting, is formed on the sur- 
face. This ‘‘raincoat”’ is so thin that its structure cannot be de- 
termined by chemical analysis. It cannot be seen under a high 
powered microscope. But, whatever its nature, it prevents water 
from spreading to form’a continuous film. If moisture does collect, 
it is in the form of small isolated drops. In the case of ceramic 
insulators, these do not allow electricity to pass because they are 
separated from each other. Another use is a laboratory one. The 
top of water in laboratory glassware, such as measuring cylinders 
and hydrometers, is ordinarily curved. Such a curved surface or 
‘‘meniscus”’ is prevented if the inside of the apparatus is treated 
with the waterproofing vapors. Then the water surface is flat and 
its height may be read more easily 


R. H. O. 


